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Preface 


For  one  unfamiliar  with  radar,  it  is  often  very  diffi¬ 
cult  to  gain  a  complete  and  thorough  understanding  of  a  parti¬ 
cular  aspect  of  radar  engineering  because  of  the  multitude  of 
disciplines  and  variety  of  backgrounds  and  viewpoints  of  the 
individuals  explaining  radar  engineering.  A  large  percentage 
of  radar  literature  deals  with  the  analysis  of  specific 
aspects  of  radar  engineering  as  opposed  to  synthesis.  I  have 
found  this  to  be  especially  true  in  understanding  the  signal 
processing  that  is  needed  to  cancel  the  clutter  present  in 
airborne  or  spaceborne  radar  returns  so  as  to  obtain  target 
range  and  velocity  information.  There  appears  to  be  two 
reasons  for  this  inadequacy:  1 )  to  adaquately  model  and  des¬ 
cribe  the  clutter  problem  a  varied  and  extensive  mathematical 
background  is  needed  and  2)  much  of  the  literature  that 
relates  to  the  clutter  cancellation  problem  either  broadly 
treats  the  subject  or  concentrates  on  a  specific  feature  or 
aspect  and  causes  one  to  lose  perspective  of  the  interrela¬ 
tionship  between  the  specific  topic  and  the  clutter  cancella¬ 
tion  problem.  It  is  thus  my  objective  to  present  a  concise 
and  succinct  system  explanation  of  the  signal  processing 
needed  to  perform  clutter  cancellation  by  blending  mathemati¬ 
cal  theory  and  modeling  of  signal  analysis  techniques  with  a 
heuristic  description  of  the  physical  problem. 

To  understand  the  clutter  cancellation  problem  it  is 

iii 


important  that  one  maintain  a  proper  perspective  between  the 
reality  of  the  physical  situation  and  the  occassional  abstrac¬ 
tion  of  mathematical  modeling.  Mathematical  models  allow  one 
to  describe  and  handle  the  clutter  problem  manageabably. 
These  models  are  idealizations  or  approximations  of  the  actual 
physical  phenomena  involved;  yet  to  be  useful  they  must  convey 
the  essential  features  of  the  situation.  The  balancing  of  the 
complexity  of  reality  against  the  simplicity  of  the  math 
models  of  that  reality  is  one  of  the  central  issues  that  I 
wish  to  explore.  While  it  is  important  to  understand  the 
theory  and  proof  of  various  mathematical  principles,  because 
many  insights  into  the  natural  order  of  events  can  be  brought 
to  light,  I  do  not  wish  to  get  bogged  down  in  elaborate  ab¬ 
straction  and  precision  of  various  particular  branches  of 
mathematics.  I  hope  to  build  upon  the  mathematical  proofs  and 
combine  physical  reality  to  give  one  a  full  appreciation  and 
grasp  of  the  problem. 

Based  upon  these  somewhat  philosophical  thoughts,  I  will 
proceed  in  the  next  several  chapters  explaining  the  math 
fundamentals  that  1 )  justify  signal  representation  and 
analysis  and  2)  provide  a  way  of  modeling  clutter.  These 
mathematical  ideas  will  then  be  applied  in  developing  the 
signal  processing  techniques  that  minimize  and  even  eliminate 
clutter  from  a  radar  return.  Hopefully,  after  reading  this 
text,  one  will  have  a  fairly  substantial  understanding  of  the 
clutter  cancellation  problem  and  can  then  effectively  channel 
one's  efforts  into  specialized  research  of  a  particular  aspect 


of  the  clutter  cancellation  problem 


Joseph  DeVenuto 
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Abstract 

The  ability  of  a  space-based  radar  surveillance  system 
to  obtain  target  information  from  a  clutter  corrupted  radar 
echo  signal  is  dependent  upon  the  clutter  cancellation 
technique  employed.  To  thoroughly  understand  the  nature  of 
clutter  so  as  to  efficiently  and  effectively  design  and 
develop  clutter  cancellation  techniques,  fundamental  mathema¬ 
tical  relationships  dealing  with  radar  signal  representation, 
probability  theory,  detection  and  estimation  theory  and  signal 
comparison  techniques  were  established  and  defined.  Since  the 
implementation  of  the  mathematical  expressions  which  describe 
clutter  and  clutter  cancellation  techniques  is  accomplished 
digitally,  digital  signal  representation  and  digital  filter 
characterization  via  FFT's  were  also  established  and  defined. 
Drawing  upon  the  established  mathematical  principals,  as  well 
as  temporal  probability  models,  spectral  models  and  average 
backscatter  coefficient  models,  clutter  was  described. 
Capitalizing  on  the  different  spectral  characteristics  between 
clutter  and  targets,  three  representative  doppler  cancellation 
schemes  were  developed.  The  three  schemes  were  delay  line 
cancelers,  doppler  filter  banks  and  multiple  look  systems. 

Using  the  models  developed  for  clutter  and  clutter 
cancellation  schemes,  an  interactive  simulation  was  developed, 
and  used  to  analyze  the  clutter  cancellation  ability  of 
various  cancellation  methods  against  different  backgrounds. 


Joint  United  States  -  Canadian  Air  Defense  Policy  states 

that  Canada  and  the  U.S.  will  jointly  maintain  an  air  defense 

system  that  can  1 )  provide  tactical  warning  and  attack 

characterization  so  as  to  limit  damage  to  strategic  and 

3 

retaliatory  forces  and  command,  control  and  communication  (C  ) 
nodes  and  2)  exercise  control  of  peace-time  access  to 
continental  airspace  (Ref  4:1).  A  viable  air  defense  system 
is  predicated  on  a  surveillance  system  that  will  provide 
timely,  reliable  and  unambiguous  warning  at  a  level  sufficient 
to  support  survival  of  strategic  forces  located  in  North 
America. 

In  October  1981,  in  support  of  the  North  American 
Tactical  Warning  (ATW)  Mission,  Air  Force  Space  Division 
proposed  a  constellation  of  space-based  radar  sensor  systems 
to  provide  indication  and  warning  of  an  impending  air  attack 
on  North  America.  Performance  of  the  surveillance  system  is 
based  upon  the  sensor's  ability  to  detect,  track  and  provide 
intercept  vectors  of  potential  hostile  aircraft.  The  radar 
system  aboard  each  space  platform  is  judged  by  its  ability  to 
extract  range  and  velocity  information  describing  a  target 
(aircraft)  from  a  returning  radar  signal. 

The  extraction  of  target  range  and  velocity  information 
would  be  fairly  straightforward  if  it  were  not  for  the  fact 
that  the  return  signal  is  corrupted  by  energy  returns  from 


unwanted  scatters.  This  unwanted  energy  limits  the  radar 
system  performance.  The  limitation  to  target  detection  and 
target  estimation  by  unwanted  energy  returns  can  be  so  severe 


that  at  times  the  target  is  completely  masked.  To  eliminate 
the  masking  of  targets,  techniques  must  be  employed  which 
reject  or  cancel  out  the  unwanted  energy  returns. 

Problem  Statement 

The  return  from  unwanted  targets  is  defined  as  clutter. 
Clutter  is  the  return,  or  group  of  returns  of  reflected  radar 
energy  which  results  from  scattering  points  of  background 
targets  such  as  buildings,  birds,  insects,  weather,  land 
and/or  sea.  The  magnitude  of  the  clutter  returns  is  a 
function  of  not  only  the  orientation  of  the  radar  antenna 
energy  pattern  with  respect  to  the  background  targets,  but 
also  the  amount  of  surface  area  illuminated.  The  radar 
antenna  energy  pattern  describes  the  magnitude  of  energy  the 
antenna  radiates  as  a  function  of  direction.  The  pattern  is 
characterized  by  a  main  lobe,  or  main  beam,  and  sidelobes. 
Generally  the  greatest  amount  of  clutter  energy  returns  result 
from  the  illumination  of  background  targets  with  the  mainlobe, 
but  since  the  surface  area  which  the  sidelobes  illuminate  is 
so  large,  the  sidelobe  clutter  return  cannot  be  ignored. 

Because  a  radar  based  on  a  space  platform  illuminates  a 
large  surface  area,  the  return  signal,  or  echo  signal  is 
dominated  by  ground,  sea,  and  weather  clutter.  The 
distribution  of  the  power  magnitude  versus  frequency  is  called 


a  power  spectral  density  (PSD) .  The  echo  signal  which  results 
from  clutter  is  generally  described  by  a  PSD  curve.  Current 
clutter  cancellation  techniques  employed  on  mobile  radar 
platforms,  such  as  aircraft,  ships,  or  trucks,  use  filtering 
schemes  or  clutter  cancellation  techniques  to  filter  out  the 
power  spectrum  which  corresponds  to  stationary  and  unwanted 
background  returns.  This  allows  for  detection  of  targets  of 
interest  such  as  aircraft  or  missiles. 


Scope 

This  study  analyzes  the  clutter  cancellation  capability 
of  various  techniques  available  that  could  be  used  in  Air 
Force  Space  Division's  space-based  radar  (SBR)  system  to 
cancel  ground  and  sea  clutter.  This  study  compares  the 
various  cancellation  techniques  presently  employed  on  mobile 
radar  systems  in  a  space-based  radar  system,  and  is  not  a 
study  intended  to  be  an  exhaustive  effort  to  collect  all 
methods  available  to  cancel  clutter.  Clutter  cancellation 
techniques  can  be  broadly  divided  into  doppler  and  nondoppler 
processing  techniques.  This  study  concentrates  on  the  doppler 
processing  techinques  that  are  described  in  Chapter  III. 

Assumptions 

Although  this  effort  was  motivated  by  Air  Force  Space 
Division's  space-based  radar  design  development,  the  satellite 
and  radar  parameters  are  generic.  The  radar  aboard  the  space¬ 
craft  is  a  pulse  doppler  phased  array  radar.  The  radar 


antenna  size,  radiated  power,  orientation,  and  transmitted 
waveform  are  described  in  Appendix  A. 

There  is  still  some  disagreement  as  to  the  exact 
statistics  of  earth  and  sea  clutter,  therefore  the  clutter 
model  defined  in  Space  Division's  Space-based  Radar  System 
Specification  was  used.  This  model  is  detailed  in  Appendix  A. 

Evalution  Standards 

The  clutter  attenuation  factor,  CA,  will  be  used  to 

evaluate  the  performance  of  the  various  clutter  cancellation 

techniques  that  could  be  employed  on  the  space-based  radar 

system.  The  clutter  attenuation  factor  is  the  ratio  of  the 

received  clutter  power  normalized  by  the  output  clutter  power 

of  a  clutter  cancellation  technique,  or  CA  =  C  /  C 

in  out 

This  criterion  was  selected  because  it  clearly  shows  the 
capability  of  a  cancellation  technique  to  remove  clutter 
energy.  (Clutter  attenuation  factor  has  the  added  advantage 
of  having  a  universal  definition  throughout  literature.)  The 
detectability  of  signal  energy  is  the  end  product  we  wish  to 
achieve,  thus  for  this  study  the  cancellation  technique  which 
yields  the  largest  CA  is  considered  as  the  optimal 
cancellation  technique.  Obviously,  other  factors  such  as 
cost,  weight,  space  qualification  and  complexity  could  have 
been  considered  in  choosing  an  optimal  technique,  but  due  to  a 
wide  variety  of  components  available  to  implement  a  specific 
technique,  cost  and  complexity  as  well  as  the  other  factors 
were  not  discussed  so  as  to  limit  the  size  of  the  study. 


Approach 

Several  steps  were  necessary  to  analyze  the  various 
clutter  cancellation  techniques: 

1 •  Various  fundamental  mathematical  relationships  dealing 
with  radar  signal  representation,  probability  theory, 
signal  comparison  techniques,  i.e.  correlation, 
spectral  analysis,  etc.,  detection  and  estimation 
theory,  and  digital  design  needed  to  all  be 
established. 

2.  Since  typically  clutter  cancelers  filter  or  mask  out 
the  power  spectral  components  which  correspond  to 
clutter,  the  design  and  implementation  of  filters  had 
to  be  understood  from  both  the  time  and  frequency 
domain  viewpoint.  The  frequency  and  time  domain 
viewpoints  allow  one  to  comprehend  exactly  what  the 
filter  is  doing  to  the  voltage  signal  that  represents 
clutter. 

3.  A  complete  and  thorough  understanding  of  the 

characteristics,  parameters  and  models  which  describe 
land,  sea,  and  atmospheric  clutter  as  seen  from  a 
space -based  radar,  had  to  be  gained.  Clutter  power 
intensity  as  a  function  of  background  had  to  be 
understood.  The  clutter  power  spectral  density  had  to 
be  thoroughly  comprehended.  An  understanding  of 
clutter  was  essential  because  clutter  cancellation 
methods  could  not  be  effectively  developed  to  cancel 
an  unknown  quantity. 
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4.  Once  clutter  and  filter  were  understood,  clutter 
cancellation  algorithms  and  techniques  could  be 


synthesized. 

5.  The  clutter  and  the  filters  which  make  up  the  clutter 
cancellation  techniques,  were  computer  modeled  to  more 
fully  understand  the  cancelling  problem  and  to  test 
the  feasibility  of  various  cancellation  techniques. 

6.  The  attenuation  factor  for  the  selected  cancellation 
techniques  were  then  computed  and  compared  to 
determine  the  optimal  technique. 

These  six  steps  will  be  used  to  answer  the  question  of 
what  is  the  optimal  clutter  cancellation  technique  to  be 
employed  on  a  space-based  radar  system. 


Summary  of  Current  Knowledge 

General  information  on  airborne  doppler  radar,  pulse 
doppler ,  pulse  doppler  clutter,  and  clutter  cancellation 
techniques  can  be  found  in  Merrill  Skolnik's  books 
Introduction  to  Radar  Systems  (Ref  56)  and  Radar  Handbook  (Ref 
57)  and  also  David  Barton's  Radar  System  Analysis (Ref  6). 
Specific  clutter  descriptions  and  clutter  models  can  also  be 
found  in  Barton's  Radar  Clutter,  Volume  5^.  This  book  contains 
a  compendium  of  articles  describing  clutter  and  clutter 
modeling.  Clutter  cancellation  techniques  and  clutter 
rejection  circuitry  can  be  found  in  Barton's  CW  and  Doppler 
Radar ,  Volume  and  D.  Curtis  Schleher's  MTI  Radar .  Both  of 
these  books  are  compendiums  of  articles  which  deal  with 


clutter  cancellation 


These  references  as  well  as  other 


support  material,  were  used  to  develop  the  foundation  of 
knowledge  that  is  presented  in  the  next  three  chapters. 
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Chapter  II  Theoretical  Background 


Introduction 

Before  developing  and  analyzing  return  clutter  waveforms 
and  the  receiver/processors  to  eliminate  clutter,  it  is 
necessary  to  establish  some  fundamental  concepts.  In  this 
chapter  I  develop  1 )  appropriate  time  and  frequency  signal 
representation  for  both  deterministic  and  random  signals;  2) 
signal  analysis  techniques;  3)  matched  filters  and  their 
usage;  and  4)  detection  and  estimation  criteria.  Since  signal 
analysis  in  radar  is  usually  accomplished  with  a  special 
purpose  digital  computer,  digital  domain  signal  representation 
and  digital  analysis  of  signals  and  systems  will  be  presented. 
It  is  not  my  intent  to  present  a  complete  and  thorough 
development  of  all  these  digital  concepts,  as  has  been 
elequently  accomplished  in  the  various  references  cited,  but 
rather  to  highlight  their  important  and  salient  features  and 
show  their  application  to  the  clutter  cancellation  problem. 
With  an  understanding  of  signal  representation  and  analysis, 
prudent  waveform  design  through  the  use  of  ambiguity  functions 
can  be  explained. 

Signal  Representation 

Signal  representation  is  probably  the  most  important  and 
fundamental  concept  needed  to  understand  both  the  analysis  and 
the  design  of  radar  engineering.  Information  about  the 
presence,  cross  section,  range,  velocity,  etc.  of  a  target  is 
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all  determined  from  the  received  reflection  of  the  transmitted 
signal.  More  efficient  signal  analysis  and  design  can  occur 
with  proper  waveform  representation.  For  the  clutter  cancel¬ 
lation  problem  a  way  of  representing  1 )  the  transmitted  and 
reflected  target  energy  2)  the  reflected  clutter  energy  and  3) 
noise  energy  is  needed. 

The  energy  transmitted  from  the  radar  antenna  and 

received  from  a  reflected  target  can  be  represented  by  a 

deterministic  signal  in  either  the  time  domain  or  in  the 

frequency  domain.  The  time  domain  representation  can  either 

take  a  real  or  complex  form  with  the  frequency  domain 

expressions  being  transforms  of  one  of  the  time  domain  forms. 

The  most  straightforward  representation  of  a  radar  signal  is 

as  a  real-valued  time  function  consisting  of  a  high  frequency 

carrier  modulated  in  amplitude  and/or  phase  by  functions  that 

are  generally  narrowband  (i.e.  B/f  £  .1  where  B  is  the  signal 

o 

bandwidth  and  f  is  the  carrier  frequency).  The  radar  signal, 
o 

denoted  s(t),  can  therefore  be  written  as 


s(t)  =  A(t)  cos[27Tf  t  +  Q(t)  +  ]  ( 1 1  —  1  ) 

o  o 


where  A(t)  is  the  amplitude  modulation  function,  9(t)  is  the 

phase  modulation  function,  f  is  the  carrier  frequency  in 

o 

hertz,  and  <f>  is  the  initial  phase  angle.  Throughout  the  rest 
o 

of  this  section  I  assume  4>  =0  with  no  loss  in  generality. 

o 

Making  use  of  the  trigonometric  identities  listed  in 
Appendix  C,  equation  II-1  can  be  presented  in  quadrature  form 


y ( t )  =  A(t)  sin  0 


In  the  development  so  far  we  have  simply  represented  the 
signal  in  a  different  form,  i.e.  as  quadrature  components, 
however  a  similar  quadrature  component  formation  is  used  in 
the  development  of  the  clutter  cancellation  techniques.  We 
will  see  that  the  formation  of  quadrature  components,  similar 
to  those  described  by  equation  II-3,  enhances  the  radar 
processor's  ability  to  recover  information  that  would 
otherwise  be  lost.  Quadrature  components  are  also  used  to 
enhance  the  detection  process. 

From  the  quadrature  component  representation,  the 
complex  exponential  representation,  u(t),  is  formed  from  the 
two  real  low  pass  signals  x(t)  and  y(t)  as 


u(t)  =  x ( t )  +  j  y { t )  =  A(t)  e 


jO(t) 


( II-4 ) 


where 


2  2  2 
A  (t)  =  x  (t)  +  y  (t) 


0(t)  =  tan  Iy(t)/x(t)J 


( II-5 ) 


In  equation  II-4,  u(t)  is  the  complex  envelope  of  the  signal 


K 


s(t).  Using  appropriate  trigonometric  identities  it  is  easy 
to  show  that  the  radar  signal  described  in  equation  II-1  can 
be  rewritten  as 

★ 

s(t)  =  .5  [u(t)  exp[j21TT  t]  +  u  (t)  exp[-j27Tf  t]]  (a) 

o  o 

( II-6 ) 

=  Re  [u(t)  exp[j277f  t ]  ]  (b) 

o 

where  Re  means  "real  part  of."  Equation  II-6b  provides  a 

form  where  only  the  complex  envelope  of  the  signal  need  be 

utilized  in  calculations,  thus  eliminating  the  need  for  any 

manipulations  of  the  carrier  frequency. 

A  different  complex  representation  of  the  radar  signal 

described  in  equation  II-1  can  be  obtained  by  application  of 

the  Hilbert  transformation.  Given  the  signal  s(t),  we  form 

the  analytic  signal  s  (t)  by 

a 

s  (t)  =  s(t)  +  j  §(t) .  .  ( 1 1  —  7 ) 

a 

The  function  §(t)  is  the  Hilbert  transform  given  by 

$(t)  =  ( 1  / 2rt)  j  s ( o'*  ) / ( t  -  o'*  )  d a*  ( II-8 ) 

where  f  denotes  the  Cauchy  principle  value  of  the  integral. 
(Note:  Integrals  without  labeled  limits  are  defined  to  have 

infinite  limits.)  From  equations  II-6b  and  II-7  we  see  that 


<  1 1  —  9  ) 


s(t)  =  Re  [s  (t)]  =  Re  [u(t)  exp[j277f  t ]  ] 
a  o 

therefore  the  analytic  signal  is  related  to  the  complex 
exponential  representation  by 

s  (t)  =  u(t)  exp[j27Tf  t]  (11-10) 

a  o 


Thus  for  analysis  of  radar  signal  we  can  use  the  complex 
exponential  signal  form  and  assume  it  is  the  same  as  the 
Hilbert  transform  representation  (Ref  20:20). 

When  using  the  complex  envelope  representation,  atten¬ 
tion  must  be  focused  on  the  relationship  between  the  energy, 
E,  of  the  signal  and  the  energy  of  the  complex  envelope.  From 
Parseval's  theorem  the  energy  of  a  signal  is 


(2  2  2 

E  =  | s ( t ) I  dt  =  A  (t)  cos  [277T  t  +  0(t)  +  4>  ]  dt 

s  )  o  O 

-  Ja  (t)  /  2  dt  (11-11  ) 


The  energy  in  the  complex  envelope  is 

(  2  (  jO(t)  *  -j®(  t) 

E  =  I  |u(t)|  dt  =  I  A  ( t )  e  A  (t)  e  dt 

=  |  A2(t)  dt  (11-12) 


Therefore  we  see  that  the  energy  in  the  complex  envelope 
contains  twice  the  energy  of  the  signal  it  represents. 


•  "•  v  *.  \  .  •  «, 
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As  mentioned  before,  the  radar  signal  can  also  be  des¬ 
cribed  by  its  frequency  content.  The  frequency  spectrum,  S(f), 


is  the  Fourier  transfrom  of  the  time  domain  function  s(t) 


S(f ) 


-j2TTft 

s  ( t )  e  dt 

/ 


(II-13) 


The  inverse  transform  is  thus  written  as 


s(t) 


j27TEt 

S(f)  e  df 

j 


(11-14) 


There  are  many  Fourier  transform  properties  which  aid  in 
the  radar  analysis  and  design.  (See  Appendix  C  page  C-2  for  a 
list.  Ref  45:  Chapter  3  or  Ref  46:  Chapter  2)  Two  particular¬ 
ly  useful  properties  are  the  symmetry  relation  and  the 
frequency  shifting  theorem.  If  s(t)  is  a  real-valued  function 
the  symmetry  relation  is 


★ 

S(f)  =  S  (-f) 


(H-15) 


* 

where  S  (f)  is  the  complex  conjugate  of  S(f).  The  shifting 
theorem  is  defined  as 


g(t)  exp[  j  2T7f  t]  o  G(f  -  f  ) 
o  o 


(11-16) 


Applying  equations  11-15  and  11-16  to  the  radar  signal, 
equation  II-1  can  be  written  as 


S(f ) 


(11-17) 


=  .5  [U(f  -  f  )  +  U(f  +  f  )] 
o  o 

where  U(f  -  f  )  is  the  Fourier  transform  of  u(t),  defined  in 
o 

equation  II-4. 

Taking  the  Fourier  transform  of  the  complex  signal  form, 
described  in  equation  II-6a,  we  get 

★ 

S(f)  =  . 5[U(f  -  f  )  +  U  (f  -  f  )]  =  U(f  -  f  )  (11-18) 

o  o  o 

From  equation  11-18  we  see  that  the  amplitude  frequency 
spectrum  of  the  complex  form  is  twice  the  amplitude  of  the 
first  component  of  the  signal  amplitude  spectrum.  Figure  II-1 

illustrates  the  spectrum  magnitude  for  the  real  ind  ,'cjjaplex 
signals  and  their  interrelationship. 


(b)  Complex  Signal  Magnitude  Spectrum 
Fig.  II-1  Real  and  Complex  Signal  Magnitude  Spectrum 


Up  to  this  point  when  describing  the  radar  signal 
frequency  spectrum  we  have  only  been  concerned  with  the 
amplitude  spectrum.  Two  other  useful  spectral  quantities  are 
the  energy  spectral  density  (ESD)  and  the  power  spectral 
density  (PSD),  defined  as 

2 

ESD  of  s ( t )  =  )  S  (  f ) }  (11-19) 

2 

PSD  of  s(t)  =  Lim  (T  -»  CO)  U(f)I  /T  (11-20) 


where  T  is  the  observation  time.  The  signal  energy  and 
average  power  can  thus  be  calculated  by  integrating  over  the 
signal  bandwidth.  The  signal  energy  is  thus 


E 


2 

!>  S  (  f  )  |  df  (Joules/Hz) 


and  the  average  power  is 


(11-21  ) 


(  2 

P  =  Lim  (T  ->C0)  |  S  ( f )  |  /T  df  (Watts/Hz)  (11-22) 

ave  ) 


The  signal  average  power  and  the  power  spectral  density 
provide  a  convenient,  and  most  useful,  way  of  describing 
clutter  and  noise  and  thereby  a  method  of  analyzing  the 
various  clutter  cancellation  schemes.  This  will  become 
evident  in  the  next  chapter. 

Because  of  the  randomness  in  amplitude  and  phase  of 


clutter  scatterers  and  the  random  nature  of  noise,  clutter  and 
noise  cannot  be  described  by  a  deterministic  signal,  but 
rather  by  a  stochastic  process.  A  stochastic  process  is 


defined  as  a  collection  of  time  functions,  x(t,c),  that  are 

the  outcomes  from  a  repeated  experiment.  The  outcomes  from 

each  experiment  repetition,  one  for  each  c  ,  form  a  family,  or 

1 

ensemble,  of  functions.  If  only  discrete  values  of  time  are 
considered,  then  x(t,c)  is  a  discrete  stochastic  process, 
otherwise  x(t,c)  is  a  continuous  stochastic  process.  (For 
notational  convenience  the  "c"  in  x(t,c)  will  be  suppressed  in 
future  discussions  when  it  is  clear  that  stochastic  processes 
are  the  issue. ) 


The 

probability 

theory  notation  to 

be 

used  will 

be 

described 

before  discussing  the  representation 

of  noise 

and 

clutter  as 

stochastic 

processes.  While 

the 

notation 

and 

statistical 

properties 

presented  assume 

continuous  random 

variables,  the  theory  applies  to  discrete  random  variables. 
To  obtain  the  discrete  form  simply  replace  the  integrals  with 
summations.  Given  two  events,  defined  by  the  random  variables 
X  and  Y  that  belong  to  a  sample  space,  the  following  probabi¬ 
lity  properties  are  defined: 

1)  The  probabilty  of  event  X  occurring,  or  the  relative 

frequency,  is  P(X)  =  Lim  (n  -»  «»  )  n  /n. 

x 

2)  The  joint  probability  of  both  events  X  and  Y  occurring 

is  P{ X , Y ) . 

3)  The  events  X  and  Y  are  independent  if  P(X,Y)  = 

P(X)P(Y) . 

4)  The  conditional  probability  of  event  X  occurring  given 

that  event  Y  has  occurred  is  P(X/Y). 

5)  Conditional  probability  can  be  written  as  P(X,Y)  = 


P(X/Y)P(Y)  =  P(Y/X)P(X). 


I 


The  distribution  of  a  function  of  the  random  variable  X 
is  defined 

P(x)  =  P [ X  £  x]  (11-23) 

defined  for  -  CO_i  x  <_  CO  (Ref  44:88).  The  probability  density 
function  (pdf)  of  X  is  the  derivative  of  the  distribution 
function . 

p ( x )  =  d  P(x)/dx  (11-24) 

where 

1)  p(x)  is  the  probability  that  the  value  of  x  lies 
between  x  and  x  +  dx. 

2)  p(x)  _>  0  for  all  values  of  x. 

3 )  J  p(x)  dx  =  1. 

4)  The  derivative  of  P(x)  is  assumed  to  exist. 

From  equations  11-23  and  11-24  we  see  that  . 

x 

P(x)=f'p(x)dx  (11-25) 

—  0O 

b 

P[a  ±  X  £  b]  =  P(b)  -  P(a)  =  J  P(x)  dx  (11-26) 

a 

where  P(b)  >  P(a).  On  page  C-3  of  Appendix  C  is  a  list  of 
the  pdfs  commonly  used  in  this  study. 

The  expected  value  (mean)  of  a  random  variable  with  a 


pd f  p ( x )  is 


E[x]  =  m 

x 


J  x  p( x )  dx 


(11-27) 


where  "E"  is  the  expectation  operator.  The  expectation  of  a 
function  f(x)  is  defined  as 


E[f(x) ] 


f 

f (x)  p(x)  dx 

j 


(11-28) 


n  n  th 

where  if  f(x)  =  x  ,  the  function  E[x  )  is  called  the  n 

moment.  Note  the  first  moment  of  p(x)  is  the  mean  of  p(x). 

The  second  moment  about  the  mean,  or  second  central  moment,  is 

the  variance  defined  as 


2  2  f  2 

a*  =  E  ( ( X  -  m  )  ]  =  l(X-m)  p(x)dx 

X  X  J  X 

2  2 

=  E [X  ]  -  m  (11-29) 

x 


The  square  root  of  the  variance  is  called  the  standard 

deviation.  The  third  central  moment  is  a  measure  of  the 

asymmetry  of  p(x)  about  m  and  is  called  the  skewness  of  the 

x 

distribution . 

The  statistical  properties  developed  for  a  single  random 
variable  can  be  extended  to  multiple  random  variables.  We  can 
think  of  multiple  random  variables  as  being  components  of  a 
random  point  in  a  multiple  dimension  space.  The  point  is 
called  a  vector  random  variable  or  random  vector.  This 
concept  will  be  exploited  in  detection  theory  when  detection 
is  based  upon  multiple  signal  returns.  Random  vector 
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statistics  is  also  utilized  when  a  continuous  signal  is 
discretized.  We  will  limit  the  discussion  of  multiple  random 
variables  to  two  variables  X  and  Y  with  no  loss  in  generality. 

If  X  and  Y  are  two  jointly  distributed  random  variables, 
with  joint  pdf  p(x,y),  the  following  properties  hold: 


cr  =  E[(X-m  )  ( Y-m  )]  =  (X  -  m  )(Y  -  m  )p(x,y)dxdy 
xy  x  y  ) )  x  Y 

11)  The  normalized  covariance  of  X  and  Y,  or  correlation 

coefficient,  is: 


E[(X  -  m  )(Y  -  m  ))  <r* 

x  y  xy 

p  =  -  =  - 

xy  2  2.5 

l  E  [  ( X  -  m  )  ]  E  [  (  Y  -  m  )  ]]  a*  0s 

x  y  x  y  . 

where  I  P  |  _£  1  .  ; 

xy 

12)  The  random  variables  are  uncorrelated  if  E[XY]  = 
E[X]E[Y].  If  the  random  variables  are  uncorrelated  then 
the  correlation  coefficient  is  zero.  If  the  two  random 
variables  are  independent  then  they  are  uncorrelated. 


13)  The  random  variables  are  orthogonal  if  E[XY]  =  0. 


A  final  statistical  operation  that  needs  to  be  explained 

is  a  method  of  obtaining  the  statistics  of  a  random  variables 

defined  indirectly  as  a  function  of  another  random  variable. 

This  occurs  for  example  when  the  quadrature  components  of  the 

received  voltage  is  known  and  the  pdf  of  the  magnitude  and 

phase  of  the  received  voltage  are  needed.  Given  a  N 

dimensional  space  X  with  joint  pdf  p(x  ...  x  ),  let  the  new  N 

1  N 

variables  be  defined  on  a  space  Y  be  defined  by  the  N  func¬ 
tions 


y  =  g  (x  ...  x  );  k  =  1,  ...  N  (11-30) 

k  k  1  N 


The  inverse  functions  are  defined  as 


X  =  h  (y  ...y);  k  =  1 ,  —  N 
k  k  1  N 


(H-31  ) 


The  new  joint  pdf  is  given  by 


P(y  . . .  y  )  = 
1  N 


p[x  =  h  (y  ...  y  );  k=1  ...  N]  U 
k  k  1  N 


(II-32) 


where  1JI  is  the  Jacobian  of  the  transformation  defined  as 


ah/dy  .  .  .  dh/^y 

11  N  1 


Ul  = 


(11-33) 


ah/dy  .  .  .  3  h  /  Jy 

IN  N  N 


(For  a  more  complete  discussion  of  the  probability  properties 
presented  consult  Ref  44:Chapters  1-8  and  Ref  45:Chapter  4.) 

The  statistical  properties  just  discussed  for  random 
variables  also  apply  to  stochastic  processes.  (For  a  complete 
explanation  of  stochastic  properties  consult  Ref  44:Chapters 
10  and  11  or  Ref  20:34-39.)  In  the  clutter  cancellation 
discussions  that  ensue,  we  will  assume  that  the  noise  and 
clutter  processes  are  wide-sense  stationary,  that  is  the  sta¬ 
tistical  properties  of  the  ensemble  do  not  change  under  a 
shift  in  the  time  origin,  but  are  a  function  of  time 
difference  T  (Ref  66:30)  and  that  the  ensemble  average  is  a 


constant 


For  a  random  process  we  define  an  ensemble  moments  and 
time  moments.  If  the  ensemble  moments  and  the  time  moments 
are  equal  then  we  say  the  process  is  ergodic  and  we  only  need 
analyze  one  set  of  moments,  usually  the  ensemble  moments.  The 
ensemble  first  moment,  or  average,  is  defined 


E  [  x  ( t )  ]  = 


x ( t  , q )  p(q , t  )  dq  = 
1  1 


x  p(x,t  )  dx  (11-34) 
1 


where  t  is  a  fixed  time.  The  time  average  is  defined  as 
1 


<  x  ( t )  >  =  Lim  (T  -»C0)  [1/2T]  ]  x(t)  dt 

-T 


(11-35) 


where  <  >  denotes  the  time  average. 

One  of  the  most  important  relationships  of  stochastic 
processes  is  that  of  correlation.  Correlation  is  a  measure  of 
interdependence  or  association  between  signal  values  (Ref 
60:233).  Again  we  have  both  an  ensemble  and  time  correlation 
function.  The  ensemble  autocorrelation  function,  the  correla¬ 
tion  between  a  function  and  a  delay  version  of  itself,  is 


★ 

R  (T)  =  E[x(t)x  ( t  +  T) J  = 
x 


JJ 


*  * 
x(t  ) X  (t  ) p [ X ( t  ) X  (t  )] 
12  12 


dx(t  )dx(t  ) 
1  2 


(11-36) 


where  t  =  t  and  t  =  t  +  T.  The  ensemble  cross-correlation 
1  2 

between  two  processes  x(t)  and  y(t)  is  defined  as 


...  .%  _  *  %  v'.  ^  .  %  . 


»  «.  ^  »  »_k  »  v  .*  »  '  »  ■  .■  »  r 


R  (T)  =  E [ x( t )y  (t  +  T)  ]  = 
xy 


x ( t  )y  (t  ) p [ x ( t  )y  (t  )]  dx ( t  )dy(t  )  (11-37) 

12  12  12 


R  (T)  =  E[y(t)x  (t  +  T)]  = 
yx 


r  f  * 

y ( t  )x 

JJ  1 


(t  )p[y(t  ) x  (t  )]  dy(t  )x(t  )  (11-38) 

2  12  12 


The  time  autocorrelation  function  is  defined  as 


IR  (T)  =  <x(t)x(t  +  T)  >  = 


Lim  (T  -»O0)  £  1  / 2T ]  x(t)x(t  +  T)  dt  (11-39) 


If  the  process  x(t)  is  ergodic  then  FR(T)  =  R(T) .  For  the 
clutter  problem  we  assume  that  the  noise  and  clutter  processes 
are  ergodic. 

The  correlation  function  has  several  properties  that  are 
useful  in  the  design  and  analysis  of  clutter  and  of 
cancellation  techniques. 


1)  For  complex  processes:  R  ( -T) '  =  R  (T) 


R  (T)  =  R  (T) 
xy  x 


2)  For  real  processes:  R  (T)  =  R  ( -T) 


R  (T)  =  R  (-T) 


3)  lR  (T) I  <  R  (0) 


1r  (T)|  1  [R  (0)  R  (0)]  1  . 5  [  R  (0)  +  R  (0)] 

x  x  y  x  y 

*  2 

4)  R  (0)  =  E  [  x  ( t )  x  (t)  ]  =  a* 


5)  R  (T)  =  x ( t )  *  x  (-t)  =  x (T)  x  (t-T)  dT 


6)  R  (T)  =  x  (-t)  *  y  ( t )  =  x  (t-T)y(T)  <JT 
xy  ) 

Analogous  to  narrowband  deterministic  signals,  a  random 
process  is  said  to  be  narrowband  if  its  spectral  density  is 
zero  except  for  a  narrow  frequency  region  around  a  high 
carrier  frequency  (Ref  66:75).  If  n(t)  is  a  real  narrowband 
WSS  random  process  with  zero  mean,  it  may  be  expressed  as 


n(t)  =  x(t)  cos27TT  t  -  y(t)  sin2Tff  t 


(11-40) 


where  the  quadrature  components  x(t)  and  y(t)  are  low  pass, 

zero  mean,  WSS  processes.  The  correlation  properties  listed 

above  also  apply  to  narrowband  processes  (Ref  66:75-78). 

The  last,  extensively  used  correlation  property  to  be 

discussed,  is  the  Wiener  -  Khinchine  theorem.  This  theorem 

states  that  the  power  spectral  density  of  a  WSS  process  x(t), 

with  autocorrelation  R  (T)  ,  is  defined  as  the  Fourier 

x 

transform  of  R  (T) .  Therefore 
x 

(  -  j  277fT 

S  (f )  =  R  (T)  e  dT  (11-41 ) 

x  J  x 
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and  it  follows  that 


S  (f)  =  S  (-f)  ;  Real  process 
xy  yx 

(11-46) 

* 

S  (f)  =  S  (f)  ;  Complex  process 
xy  xy 


As  will  become  evident  in  the  next  several  sections  and 
then  in  Chapter  III,  the  power  spectral  density  of  a 
deterministic  signal  or  a  random  process  provides  a  useful  way 
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of  describing  not  only  a  waveform,  but  also  a  way  to  describe 
and  develop  signal  analysis  techniques  that  exploit  the 
different  power  spectral  characteristics  of  signal  so  as  to 
distinguish  target  information  from  clutter  and  noise. 

Before  proceeding  on  to  the  development  of  signal 
analysis  techniques,  we  first  need  to  provide  a  descriptor  for 
white  noise.  We  will  delay  discussion  of  colored  noise,  i.e. 
clutter,  because  it  will  be  covered  in  Chapter  III.  In  the 
clutter  problem  we  are  dealing  with  bandpass  signals,  there¬ 
fore  white  noise  can  be  described  by  the  bandpass  signal 
defined  in  equation  11-45,  where  x(t)  and  y(t)  are  lowpass, 
zero  mean,  WSS  Gaussian  processes.  We  define  noise  to  have  an 
impulse  autocorrelation  function  and  an  unit  spectral 
response.  For  this  study  we  characterize  white  noise  as  show 
in  Figure  II-2. 


(a)  Autocorrelation  (b)  Spectral 

Fig.  II-2  Noise  Autocorrelation  and  Spectral  Characteristics 


Signal  Analysis 


As  would  be  expected,  since  deterministic  and  random 
signals  are  decribed  in  both  the  time  and  frequency  domain, 
signal  analysis  is  also  performed  in  both  domains.  The  signal 


analysis  techniques  used  to  cancel  clutter  and  obtain  target 
information  from  a  reflected  radar  signal,  are  based  in  one 
form  or  another  on  the  correlation  principle.  Correlation  is 
the  process  of  determining  the  mutual  relationship  which 
exists  between  two  signals  in  bandwidth,  time  duration,  shape 
and  energy  distribution,  etc.  (Ref  13:92).  As  shown  in  the 
last  section,  correlation  applies  to  both  deterministic  and 
random  functions,  therefore  the  correlation  concepts  will  be 
described  assuming  a  deterministic  signal.  In  the  cases  where 
s(t)  is  a  random  signal,  i.e.  clutter  or  noise,  then  moment 
calculations  must  be  applied.  Target  information  such  as 
range  and  range  rate  is  related  to  various  signal  parameters, 
i.e.  delay  and  frequency  shift.  Given  the  target  information 
desired,  the  appropriate  signal  parameter  from  the  reflected 
radar  signal,  is  compared,  or  correlated,  with  the  transmitted 
signal.  The  level  of  comparison  provides  an  indication  of  the 
value  of  the  signal  parameter  and  thus  the  desired  target 
information. 

From  the  time  domain  perspective,  correlation  is 
accomplished  via  the  concept  of  convolution  between  a  linear, 
time  invarient  system  and  a  signal.  From  the  previous  section 
we  saw  that  convolution  is  related  to  correlation  by 

* 

R  (t)  =  s ( t )  *  s  (-t)  (II_47) 

s 

If  s(t)  is  real  and  symmetric,  then  convolution  is 
correlation.  Given  a  signal  s(t),  a  linear  system,  or  filter. 
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can  be  derived  such  that  a  prescribed  response  occurs  for 
perfect  correlation.  Any  deviation  from  the  desired  response 
indicates  other  than  perfect  correlation.  The  response  of  the 
linear  filter  to  any  signal  is  given  by  the  convolution 
integral 

y(t)  =  s m  h(t  -  *n  dr  ai-48) 

The  convolution  between  the  signal  and  the  filter  is  denoted 
s(t)*h(t).  To  obtain  specific  target  information,  a  filter 
is  designed  to  provide  a  measure  of  correlation  of  the  appro¬ 
priate  signal  parameter  between  the  transmitted  and  reflected 
radar  signal.  The  filter,  operating  in  the  presence  of  noise 
and  clutter,  exploits  the  characteristic  differences  between 
the  signal  and  clutter  and  noise.  A  convenient  way  of 
designing  the  filter  to  perform  the  desired  correlation  on 
either  random  or  deterministic  signals  is  accomplished  via 
frequency  domain  analysis. 

Frequency  domain  analysis  can  be  approached  from  two 
points  of  view.  The  first  perspective  looks  at  frequency 
domain  analysis  as  transforming  a  signal  into  a  new  domain 
and  then  observing  the  amplitude,  energy, or  power  distribution 
over  the  frequency  spectrum  for  certain  distribution 
characteristics.  Unwanted  characteristics  are  removed  through 
selective  nulling,  or  filtering.  This  perspective  is  useful 
when  dealing  with  both  random  and  deterministic  signals.  The 
second  frequency  domain  perspective  views  frequency  analysis 


as  a  tool  to  conveniently  design  the  appropriate  correlation 


filters  that  are  suggested  in  the  time  domain.  Applying  the 
Fourier  transform,  equation  11-13,  to  the  convolution 
integral,  the  frequency  domain  representation  is 


Y(f ) 


f  -277ft  ( 

e  dt  s(T)h(t-T)  dT 


(11-49 ) 


Letting  p  =  t-T,  equation  11-48  can  be  rewritten  as 


Y(f ) 


'  -j27Tf(p+T) 

s  (T)e  dp 

4 


|  s(T)h(p)  dT 


( II-49a ) 


Rearranging  we  see 


-  j  27Tfi  (  -271fp 

Y ( f )  =  |  s(T)e  dT  |h(p)e  dp 


1 


=  S(f)  H ( f ) 


(11-50) 


that  the  time  domain  convolution  of  two  functions  is  simply 
the  multiplication  of  the  Fourier  transform  of  each 
deterministic  signals.  The  response  of  a  filter  to  a  random 
signal  is  described  in  terms  of  the  power  spectrum.  Driving  a 
deterministic  filter,  with  transfer  function  H(f),  with  a 
random  signal  x(t),  having  a  PSD  of  S  (f),  the  filter  output 


2 

S  (f)  =  S  (f)  1 H( f ) 1 


(11-51  ) 


(Ref  45:347).  Equation  11-50  also  illustrates  the  fact  that 
in  the  frequency  domain  the  necessary  filter  is  the  ratio  of 
the  desired  frequency  response  divided  by  the  signal 
frequency.  Both  perspectives  on  frequency  analysis  will  be 
utilized  throughout  this  study.  Spectral  characteristics  will 
be  used  in  describing  the  signal,  noise  and  clutter,  while 
frequency  domain  synthesis  techniques  (Refs  29,  44,  46)  will 

be  used  to  design  the  necessary  filters. 

Since  the  signal  can  conveniently  be  described  in  the 
frequency  domain,  and  analysis,  design,  and  synthesis  of 
filters  tend  to  be  easier  in  the  frequency  domain,  spectral 
analysis  will  be  heavily  utilized  in  describing  and  analyzing 
the  clutter  cancellation  problem.  Estimation  of  the  power 
spectral  density  of  a  signal  is  traditionally  calculated  by 
one  of  three  methods  (Ref  7:250):  1)  Bandpass  filtering  -  the 
signal  is  passed  through  a  bank  of  filters,  or  a  comb  filter, 
the  response  is  then  squared  and  scaled  by  the  observation 
time.  The  filter  bank  can  either  be  implemented  with  lumped 
constants  or  digitally.  Digital  filtering  will  be  discussed 
in  an  upcoming  section.  2)  Indirect  method  -  the  signal 
autocorrelation  function  is  estimated,  usually  by 

N-m-1  * 

R  (m)  =  ( 1 / (N-m) ]  2  s(nT  +  mT)s  (nT)  (11-52) 

s  n=0 

where  n  is  the  sample  number  and  m  is  the  frequency  index.  The 
PSD  is  calculated  by  taking  the  Fourier  transform  of  the 
correlation  function  usually  via  a  fast  Fourier  transform 


(FFT)  operation  such  as  the  one  suggested  by  Blackman  and 
Tukey 


M 

S  (f)  =  2  R  (m)  exp[ -j  27^mat  ]  (11-53) 

BT  m=-M  s 


(Ref  34:1383).  3)  Direct  method  -  the  Fourier  transform  of 
the  signal  is  calculated,  squared,  and  scaled,  usually  through 
a  FFT  algorithm.  FFT  algorithm  usage  will  be  discussed  in  the 
digital  filtering  section.  The  direct  method  can  also  be 
implemented  by  one  of  several  modern  estimation  methods,  the 
majority  of  which  are  based  on  the  rational  model  approach. 
In  the  rational  model  method  the  PSD  of  a  time  series  is 
assumed  to  be  modeled  as  a  rational  funtion 


jw 
S(e  ) 


-jw 

b  +  b  e  + 
o  1 


-jwn 

+  b  e 
n 


2 


- jw  - jwm 

a  +  ae  +  .  .  .  +  a  e 
o  1  m 


(11-54) 


where  m  >  n  (Ref  34:1387).  The  PSD  is  determined  through  the 
selection  of  the  coefficients.  The  use  of  nontrad itonal 
spectral  estimation  techniques  for  clutter  and  target 

estimation  will  not  be  considered  within  this  study. 

Having  established  the  fundamental  radar  signal  analysis 
techniques,  we  can  now  look  at  specific  filters  that  maximize 


the  signal  to  interference  and  allow  for  detection. 


Matched  Filter  Development 

In  the  last  section  we  showed  that  target  information  is 
obtained  by  correlating,  or  convolving,  the  received  signal 
with  a  linear  filter.  Proceeding  from  that  framework, 
effective  extraction  of  range  and  velocity  information  from 
the  echo  signal  can  be  accomplished  by  maximizing  the  signal 
to  interference  ratio,  S/I.  Interference  is  defined  as  the 
energy  from  white  noise  as  well  as  from  unwanted  scattering 
elements  (clutter  or  colored  noise).  It  can  be  shown  that  the 
signal-to-noise  ratio  is  indirectly  related  to  radar 

performance  measures  such  as  probability  of  detection, 

probability  of  false  alarm,  and  range  and  doppler  estimation 
accuracy  (Ref  22:143). 

A  radar  signal  processor  is  nothing  more  than  a 
composition  of  various  filters.  Before  developing  the  filter 
that  maximizes  the  S/I  ratio,  let’s  first  philosophically  look 
at  the  role  of  the  processor  to  aid  in  the  filter  development 
and  to  provide  insight  into  the  detection  and  estimation 
problem.  The  function  of  the  radar  signal  processor  is  to 
determine  whether  a  target  is  present  or  absent,  at  what  range 
and  at  what  velocity.  Detection  criteria  is  normally  based  on 
the  concept  of  direct  probability.  Direct  probability  is  the 
probability  which  describes  the  chance  of  an  event  happening 
before  the  event  occurs.  If  the  event  actually  happened,  then 
the  best  estimate,  or  hypothesis,  as  to  the  cause  of  the  event 
illustrates  the  concept  of  inverse  probability  (Ref  56:377). 
Detection  criteria  based  on  inverse  probability  is  particular- 


ly  useful  in  statistical  detection  and  statistical  parameter 
estimation.  To  illustrate  the  above  ideas,  take  for  example 
the  event  of  the  radar  receiver  output  voltage.  It  is  of 
interest  to  determine  whether  the  output  was  caused  by  noise 
or  by  a  target  in  the  presence  of  noise.  The  probability  that 
the  cause  of  the  event  was  noise  or  signal-plus-noise  before 
the  event  occurs  is  labeled  the  "a  priori  probability,"  while 
the  probability  that  the  receiver  output  voltage  was  caused  by 
noise  or  signal-plus-noise  after  the  event  occurred  is  called 
the  "a  posteriori  probability."  The  "a  posteriori"  probabili¬ 
ty  represents  the  state  of  information  obtained  as  a  result  of 
observing  the  output.  The  method  of  inverse  probability  uses 
the  "a  priori  probabilities"  associated  with  each  of  the 
possible  hypothesis  (i.e.  whether  the  receiver  output  was 
caused  by  noise  or  signal-plus-noise)  to  explain  the  cause  of 
the  event.  The  "a  priori  probabilities",  along  with  a 
knowledge  of  the  event,  are  used  to  compute  the  "a  posteriori 
probabilities."  A  separate  probability  is  computed  for  each 
hypothesis  and  the  largest  "a  posteriori  probability"  is 
selected  as  the  most  likely  cause  of  the  event  (Ref  56:377). 

The  best  any  receiver  can  do  is  to  form  the  "a 
posteriori  probabilty"  distribution  (Ref  67:31).  To 
illustrate  the  idea  of  inverse  probability,  the  "a  posteriori 
probability"  of  the  delay,  T,  of  the  transmitted  waveform  will 
be  formed.  As  stated  before,  these  signal  parameters  provide 
information  on  the  target's  range.  The  probabilty  distribu¬ 
tion  information,  most  conveniently  handled  in  logarithmic 
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form,  is  called  the  information  function.  The  information 
function  for  the  time  delay  distribution  is 

Q(T)  =  In  [P<Y/T)]  (11-55) 

where  Y  is  the  received  signal  voltage  with  a  time  delay 

associated  with  it,  and  P(Y/T)  is  the  probability  density  for 

the  received  signal  given  a  delay  T  (Ref  67:31). 

In  the  case  where  the  echo  return  of  the  waveform  has  a 

time  delay  of  T  the  received  signal  is 
o 

Y ( t )  =  s ( t  -  T)  +  n(t)  (11-56) 

where  s(t)  is  the  transmitted  signal  and  n(t)  is  the  added 
white  Gaussian  noise  with  probabilty  distribution  proportional 
to 

T  2 

k  exp[1/N  5  (n( t ) ]  dt]  (11-57) 

o  0 

where  N  is  the  mean  noise  per  unit  bandwidth  and  T  defines 
o 

the  observation  period  (Ref  67:31).  The  noise  function  can  be 
rewritten  as  n(t)  =  Y(t)  -  s(t-T).  Substituting  into  equation 
11-57,  we  obtain 

T  2 

k  exp[-1/N  $  t Y( t )  -  s ( t-T)  ]  dt  ]  (11-58) 

o  0 


Equation  11-58  is  proportional  to  the  probability  density  of 


the  noise  waveform  hypothesised  on  the  echo  signal  being 
delayed  T.  The  information  function  concerning  range  is  thus 

T  2 

Q(T)  =  -k/N  $  [Y(t)  -  s ( t-T)  ]  dt  (11-59) 

o  0 

The  actual  time  delay  is  obtained  by  maximizing  the 
information  function.  This  is  done  by  varying  T  to  obtain  the 
value  that  gives  the  least  mean  square  departure  from  the 
received  signal. 

Expanding  the  integral  in  equation  11-59,  we  get 

T  2  T  T  2 

S[Y(t)]  dt  -  2  $  Y  ( t )  s(  t-T)  dt  +  S  [  s  ( t--T)  ]  dt  (11-60) 

0  0  0 

Upon  reception,  Y(t)  is  known  and  can  be  absorbed  by  the  con¬ 
stant  k.  Integrating  over  one  period  makes  s(t-T)  independent 
of  time  delay  and  therefore  is  known  and  can  also  be  absorbed 
by  the  constant  k.  Thus  equation  11-60  becomes 

T 

Q(T)  =  2k/N  5  I Y ( t )  s(t-T)]  dt  (11-61) 

o  0 

Equation  11-61  is  the  information  function  which  in 
logarithmic  form  gives  the  ”a  posteriori  probability”  distri¬ 
bution  of  the  values  of  T.  From  this  function  we  obtain  the 
probability  that  the  received  signal  is  from  a  target  with  a 
time  delay  of  T.  A  similar  development  for  determining  the  "a 
posteriori”  distribution  for  the  received  waveform's  frequency 
shift  can  also  be  performed. 


The  information  function,  equation  11-61,  is  simply  the 
output  from  a  linear  filter,  which  apart  from  the  scaling 

factor  2k/N  ,  is  the  convolution  of  the  input  Y(t)  and  the 
o 

impulse  response  h(t)  =  s(-t).  The  filter  is  the  unique 
linear  filter  which  gives  the  maximum  peak  signal  to  noise 
ratio.  The  filter  has  the  same  frequency  response,  which  in 
magnitude  has  the  same  shape  as  the  amplitude  of  complex 
envelope  of  the  transmitted  signal  s(t),  but  in  phase  is  equal 
in  magnitude  and  opposite  in  sign  (Ref  67:32).  The  filter 
described  in  Woodward's  information  function  (equation  11-61) 
is  commonly  called  the  matched  filter. 

While  the  matched  filter  can  be  developed  in  the  time 
domain  (Ref  22:144-149),  it  is  more  easily  shown  in  the 

frequency  domain  by  applying  the  Schwartz  inequality  (Ref 

46:159-162) . 

In  general  the  received  signal  is  the  sum  of  the 

transmitted  waveform  s(t)  and  noise  n(t) 

y(t)  =  s(t)  +  n(t)  (11-62) 

where  n(t)  can  either  be  white  or  colored  noise.  The  signal 
y(t)  is  applied  to  the  receiver  filter  with  a  transfer 

function  H(f).  The  signal  component  at  the  filter  output  is 

r  j2‘77Jt 

s  (t)  =  S(f)  H ( f )  e  df  (11-63) 

o 

where  S(f)  is  the  Fourier  amplitude  spectrum  of  s(t).  If  t 


is  the  time  that  the  output  signal  is  maximum,  then  the  peak 
signal  power,  P,  is 


,  1 2  I 

(  , 

P  =  1  S(t)|  =  1 

I  S(f)  H(f)  exp[ j 277ft  ]  df  | 

(11-64) 


The  average  filter  output  noise  power  is 

2 


N  = 
o  j 


S  (f)  |  H( f ) |  df 
N 


(11-65) 


where  S  (f)  is  the  receiver  input  noise  density  spectrum.  The 
N 

signal  to  noise  power  to  be  maximized  is 


'1 


S(f)  H(f)  exp [277ft  ]  df 

o 


(11-66) 


N 


j  V£’  1 


H(f)|  df 


To  maximize  equation  11-70,  the  Schwarz  inequality  is  applied. 
If  A(f)  and  B(f)  are  general  complex  functions  of  the  real 
variable  f,  the  inequality  is 


|  jA(f)  B(f)  df|  ^  A(f)|  2df  |  |  B(f)| 


df  (11-67) 


where  the  equality  holds  only  when  B(f)  is  proportional  to 
A(f).  Making  the  substitutions 


.  .'"'O-Y'-  ^  I'V.V  .  -  .  .  . 


B(f)  =  S(f)  exp[277f  t  ]  /  (2  yjs  (f ) ) 

o  N 

into  equation  11-66  the  signal  to  noise  ratio  becomes 


S  / 
o 


N  <_  |  £(£  )| 

O 


/  s  (f)  df 

N 


(11-68) 


The  maximum  occurs  when  the  equal tiy  holds.  The  optimal 
transfer  function,  or  matched  filter,  is  therefore 


H  (f)  =  k  S  (f)  exp[-2TTft  ] 
opt  1  o 

S  (f) 

N 


(11-69) 


The  key  feature  of  equation  11-69  is  that  the  optimal  filter 
is  proportional  to  the  ratio  of  the  complex  conjugate  of  the 
signal  spectrum  to  the  noise  power  spectrum. 

If  the  noise  is  white,  the  matched  filter  becomes 


H  (f)  =  k  S  (f)  exp[ 277ft  ] 
opt  o 


(11-70) 


Taking  the  inverse  Fourier  transform,  the  impulse  response  is 


h  (t)  =  k  s  (t  -  t) 
opt  o 


(II-71  ) 


If  the  noise  spectrum  is  composed  of  white  and  colored  noise 
then 


s  (f)  =  N  /2  +  _ _  (11-72) 

N  o  2  2 

<*  +  (2m) 


where  *  and  Q  are  the  coeficients  of  the 

autocorrelation  function  R  (T)  =  of  exp[  -/SI  tIi- 
filter  is  thus 


colored  noise 
The  matched 


2  2  * 

2k  t(277f)  +/?  ]  S  (f)  exp[-j2TTEt  ] 

H  (f)  =  - £- -  o  (11-73) 

opt  2  2 

N  [  ( 2TTE )  +  u  ] 

o 

2 

where  u  -  +  4*^/N  (Ref  22:194). 

If  the  causality  delay,  t  ,  is  set  to  zero,  the 

o 

convolution  between  s(t)  and  its  matched  filter  is  simply  the 
s(t)  autocorrelation  function 


s(t)*h(t)  =  k  I  s(T)  s  (T-t)  d T=  R  (-T)  (11-74) 


and  the  crosscorrelation  between  s(t)  and  its  matched  filter 
is  simply  s(t)  convolved  with  itself 


* 

R(T)  =  s(t)*h  (-t) 


s  (T)  s(t-T)  d T 


(II-75) 


We  thus  see  that  the  matched  filter  response  provides  the 
correlation  measure  necessary  to  obtain  an  indication  of  the 
desired  target  information. 


The  majority  of  the  filters  in  the  radar  receiver/signal 


k  S  (f)  exp[-j277ft  ] 

H  (f)  =  _ o  (11-79) 

opt  «  I  2 

N  /2  +  p|S(f)| 


By  applying  the  shifting  theorem,  equation  11-16,  to  any 
of  the  matched  filter  transfer  functions  we  can  see  how  the 
matched  filter  can  be  used  to  constuct  the  filter  banks  that 
maximize  S/I  ratio  and  allow  for  doppler  detection.  The 
matched  filter  usage  in  the  detection  concepts  will  be  shown 
in  the  next  section. 


Detection  Processor  Development 

In  the  previous  section  the  idea  of  inverse  probability 
was  utilized  in  looking  at  the  radar  signal  processor's  func¬ 
tion.  The  role  of  the  processor  is  essentially  reduced  to 
computing  the  "a  posteriori  probabilities”  via  a  matched  fil¬ 
ter.  It  was  then  shown  that  this  filter  maximized  the  signal 
to  interference  ratio.  Since  radar  performance  measures  are 
directly  proportional  to  the  signal-to-interference  ratio, 
maximization  of  this  ratio  is  beneficial.  The  ability  of  a 
radar  to  detect  and  estimate  target  parameters  is  a  widely 
used  performance  measure.  Before  proceeding  with  the  detec¬ 
tion  processor  development,  the  interrelationship  between 
signal  analysis  and  the  use  of  matched  filters  to  maximize 


signal  to  interference  ratio  needs  to  be  clarified. 

The  term  receiver/processor  is  used  to  indicate  mathema- 


tical  operation(s)  performed  on  a  signal  corrupted  by 

interference  to  obtain,  or  format  to  obtain,  specific  informa¬ 
tion  (Ref  66:155).  The  "processor"  referred  to  in  the  last 
section,  dealt  with  the  mathematical  operations  that  are  used 
to  cancel  clutter  from  the  echo  waveform  and  then  prepare  the 
echo,  still  corrupted  by  noise,  to  be  analyzed  in  either  the 
time  or  frequency  domain.  It  will  be  seen  in  the  next  chapter 
that  this  type  of  processing  merely  segments  the  reflected 
radar  signal  into  range  and  doppler  bins,  discards  those  bins 
that  contain  clutter  and  then  passes  the  signal  to  be 
processed  so  as  to  extract  target  information.  Extraction  of 
target  information  can  be  divided  into  detection  and  esti¬ 
mation.  Detection  is  the  determination  of  a  target's 
presence,  while  estimation  is  the  process  of  obtaining  target 
range  and  velocity.  Once  a  target  is  detected,  the  echo 
signal  can  then  be  processed  more  efficiently  to  estimate 
range  and  velocity.  Since  this  study  is  limited  to  only  dop¬ 
pler  cancellation  techniques,  it  will  be  assumed  that  range 
and  velocity  estimation  is  obtained  as  a  result  of  target 
declaration  in  a  specific  doppler-range  bin.  For  this  reason 
target  estimation  techniques  will  not  be  discussed  and  the 
reader  is  referred  to  Ref  66:  Chapter  10  for  more  infomation 
on  estimation  theory.  The  remainder  of  this  section  will  deal 
with  the  detection  process. 

Since  it  is  assumed  that  detection  processing  occurs 
after  clutter  has  been  cancelled,  only  the  detection  processor 
in  the  precence  of  white  noise  need  be  discussed.  The  optimal 
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detection  processor  is  that  processor  that  best  satifies  given 
probability  of  detection  and  probability  of  false  alarm  cri¬ 
teria  under  a  given  set  of  detection  assumptions  (Ref  66:155). 
For  the  radar  detection  problem,  the  detector , utilizing  the 
concept  of  inverse  probability  previously  developed,  computes 
the  likelihood  ratio  for  the  observed  voltage 


p( v/H  ) 
l(v)  =  1 


p( v/H  ) 

o 


(11-80) 


where  p(v/H  )  is  the  "a  posteriori  probability"  that  the 
“  1 

voltage  output  from  the  previous  processing  was  caused  by  the 

reflection  from  a  target  and  p(v/H  )  is  the  "a  posteriori 

o 

probability"  that  the  output  voltage  was  caused  by  noise 
alone.  The  likelihood  ratio  is  then  compared  with  a  given 
threshold,  K, 


H 

p( v/H  )  1 

1(V)  =  1  >  K 

< 

p( v/H  )  H 
o  o 


(11-81 ) 


to  determine  if  a  target  is  present.  The  threshold  value  for 
the  radar  problem  is  determined  by  applying  the  Neyman  Pearson 
criteria.  This  criteria  maximizes  the  probability  of  detec¬ 
tion  for  a  given  false  alarm  probability.  From  equation  11-81 
we  see  that  if  the  likelihood  ratio  exceeds  the  threshold 

value  K  then  hypothesis  H  is  assumed  true  and  a  target's 

1 


presence  is  declared,  likewise  if  the  ratio  value  is  less  than 

the  threshold  then  hypothesis  H  is  assumed  true  implying  that 

o 

the  voltage  was  the  result  of  noise.  Knowledge  of  the  likeli¬ 
hood  ratio  thus  permits  specification  of  the  test  statistic 
which  describes  the  structure  of  the  optimal  detection 
processor.  The  performance  of  the  processor  can  also  be 
evaluated  by  analyzing  the  statistical  properties  of  the  test 
statistic  (Ref  22:291). 

The  likelihood  expression  for  a  single  pulse,  with  known 

signal  amplitude  and  phase,  in  the  presence  of  white  Gaussian 

noise  will  be  derived.  The  likelihood  expressions  for  the 

cases  where  the  phase  and  amplitude  are  unknown  and  also  the 

pulse  train  cases  easily  follow  from  the  simple  case  that  will 

be  presented.  These  various  signal  cases  are  extensively 

covered  in  Ref  22: Chapters  9,10,  and  11. 

To  derive  the  likelihood  ratio  it  is  convenient  to 

employ  the  sampling  theorem  waveform  representation  instead  of 

the  more  general  Karhunen-Loeve  expansion  (Ref  22:291).  Given 

a  low  pass  bandlimited  (-f  ,f  )  signal,  v(t),  through  applica- 

c  c 

tion  of  the  sampling  theorem,  sampling  every  1/2f  seconds, 

c 

v(t)  can  be  completely  specified  by  the  samples  using  the 
expression 


v(t)  = '  2  v(k/2f  )  sinc[ 277f  (t  -  k/2f  )] 
k=-«*  c  c  c 


(11-82) 


(Ref  22:53).  From  the  sampling  theorem  it  can  be  shown  that 
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(11-83) 


where  T  is  the  signal  duration  (Ref  22:56). 

Assuming  the  signal  is  known  and  bandlimited,  (-f  ,f  )» 

c  c 

with  bandlimited  noise,  the  detector  processor  input  waveform, 
v(t),  is  written  in  vector  form  as 


V* 


m  * 


i 


rv 


v  =  s  +  n 


(11-84) 


where  v  =  v(k/2f  ),  s  =  s(k/2f  ),  and  n  =  (k/2f  )  for 

k  c  k  c  k  c 

k=1,2,...,2f  T.  The  conditional  joint  density  function  that 
c 

describes  the  hypothesis  H  (signal  plus  noise)  is  represented 

1 

by 


1 


p( v/H  ) 
1 


N/2 

(277)  (N  f  ) 
o  c 


N 

~exp[-1/f  N  2 
N  co  k=1 


(v  -s  )  ] 
k  k 


(11-85) 


where  N=2f  T  and  the  white  noise  PSD  is  N  /2.  Similarly  the 
c  o 

conditional  joint  density  that  describes  the  hypothesis  H 

o 

(noise  alone)  is  represented  by 


1 


p( v/H  )  = 


N/2 

(2m  (N  f  ) 
o  c 


N  2 

—  exp[-1/2f  N  2  (v  )  ) 
N  co  k=1  k 


(11-86) 
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From  equation  11-81  we  see  that  the  likelihood-ratio 
expression  is 


A(v) 


p(v/H  ) 
1 


exp[-1/2f  N  2  (v  -s  )  ] 
c  o  k=1  k  k 


p( v/H  ) 
o 


expt -1 /2f  N 
c  o 


N 

2 

k=1 


2 

(v  )  ] 
k 


N  2  N 

=  exp[-1/2f  N  (2  (s)  -2  2  vs)]  (11-87) 

c  o  k=1  k  k=1  k  k 


Using  equation  11-83 ,  equation  11-87  may  be  rewritten  as 


T  2  T 

A(v)  =  exp[-1/N  S  s  (t)dt  +  2/N  S  v(t)s(t)dt  _>  K  (11-88) 

o  0  o  0 


Taking  the  natural  logarithm  of  equation  11-88  will  not  alter 
the  inequality  of  equation  because  the  function  In  x  is  a 
monotonically  increasing  function  of  x,  thus  equation  11-88  is 


H 

T  1 

2/N  S  v(t)s(t)  dt  > 
o  0  < 

H 


o 


K  +  E/N 

o 


(11-89) 


The  test  statistic,  the  cross  correlation  between  v(t)  and 
s(t),  is  realized  via  a  matched  filter  as  suggested  in  the 
last  section.  If  the  output  of  the  matched  filter  exceeds  the 
threshold  K  +  E/N  then  a  target  is  detected  and  tagged  for 


further  processing. 


The  detection  process  just  developed  was  for  the  case  of 
white  Gaussian  noise.  If  the  noise  is  assumed  to  be  nonwhite, 
e.g.  clutter,  the  inverse  probability  concept  can  again  be 
applied  to  form  the  structure  of  the  optimal  detection 
processor.  The  processor  based  upon  the  likelihood-ratio  to 
effectively  cancel  clutter,  is  an  example  of  nondoppler 
processing.  In  an  attempt  to  limit  the  size  of  the  study, 
nondoppler  processing  will  not  be  examined,  however  a  deployed 
system  must  use  a  combination  of  nondoppler  and  doppler 
processing  to  cancel  clutter. 

It  must  be  remembered  that  the  detection  processor 
developed  is  applied  to  each  range-doppler  bin  to  determine  a 
target's  presence.  As  mentioned  before  the  filters  to  perform 
all  the  processing  are  generally  matched  filters,  matched  to 
each  test  cell.  Because  of  the  large  number  of  bins  a 
convenient  way  of  implementing  the  filters  other  than  with 
lumped  constants  must  be  obtained.  With  the  advance  of 
microelectronics,  special  purpose  digital  computers  are  used 
to  implement  the  filters  to  perform  the  time  and/or  frequency 
domain  analysis.  The  next  section  examines  digital  domain 
signal  representation  and  analysis. 


Digital  Signal  Processing 


Following  the  pattern  previously  established  when 
describing  continuous  signals  and  systems,  this  section  first 


covers  the  representation  of  digital  signals  and  then 

describes  digital  signal  analysis  techniques.  Because  we  are 

dealing  with  the  radar  problem,  we  assume  that  the  digital 

signals  and  filters  described,  represent  sampled  versions  of 

continuous  analog  signals  and  systems. 

In  the  strict  sense,  a  digital  signal  is  a  signal  which 

has  discrete  values  in  time  and  in  amplitude.  However,  in  the 

ensuing  discussion,  digital  signals  will  denote  dicrete  time 

values  of  a  continuous  amplitude  valued  signal.  Since 

amplitude  quantinization  levels  are  the  result  of  system 

choices,  this  limitation  will  not  restrict  the  digital 

processing  dicussions.  This  limitation  also  eliminates  the 

need  for  discussion  of  quantinization  noise,  which  can  be 

found  in  Ref  29: Chapter  4  and  Ref  43: Chapter  9. 

A  continuous  analog  signal  is  described  digitally  by 

applying  the  baseband  sampling  theorem.  This  theorem  states 

that  given  a  lowpass,  bandlimited  ( -f  ,f  )  signal  v(t),  the 

c  c 

signal  can  be  uniquely  characterized  over  all  time  by  the 

samples  of  v(t),  if  the  samples  are  taken  at  intervals  of 

1/2f  (Ref  22:53).  The  frequency  2f  ,  is  the  Nyquist  frequen- 
c  c 

cy  and  sampling  at  less  than  the  Nyquist  rate  results  in 

aliasing.  The  sampling  can  be  accomplished  by  1)  natural 

sampling  :  v(nT)  =  v  (t)*sinc[2T7E  t),  2)  flat  top  sampling  : 

a  c 

v(nT)  =  v  ( t )*rect[ t/T] ,  or  3)  impluse  sampling  :  v(nT)  = 
a 

v  (t)*  (t);  where  v(nT)  is  the  sampled  version  of  the  analog 
a 

signal  v  (t).  There  are  N  samples  taken  every  T  seconds, 
a 

While  the  digital  processing  discussions  will  assume  the 


impulse  sampling  method,  in  practical  situations  natural  or 
flat  top  sampling  would  be  used.  (Note:  In  the  limit  as  T  -*  0 
both  flat  top  and  natural  sampling  reduce  to  impulse  sam¬ 
pling.)  Figure  II-3  illustrates  the  impulse  sampling  of  a 
periodic  continuous  signal  and  the  corresponding  frequency 
domain  representation. 
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(Ref  59:243) 


Fig.  II-3  Impulse  Sampling  and  its  Fourier  Transform 


The  relationship  between  the  sampling  period,  T,  and  the 


sampling  frequency,  f  , 


T  =  1/f 


(11-90) 


and  the  relationship  between  the  observation  time,  t  , 

P 

(generally  t  corresponds  to  the  signal  period  for  a  periodic 
P 

signal)  and  the  frequency  domain  frequency  interval,  F, 


t 


1  /F 


(11-91  ) 


P 

both  illustrated  in  Figure  II-3,  will  be  used  when  describing 
the  discrete  Fourier  transform  ( DFT ) . 

From  Figure  II-3  it  is  also  seen  that  the  sampled  signal 
can  be  represented  by 


v(nT)  =  2  v(nT)  5<fc  -  nT)  (11-92) 

n=-*> 


The  Fourier  transform  of  equation  11-92  is 


-  (  -jzmt 

v(nT)  =  2  v(nT)  5(^~nT)  e  dt 

n=-«  J 

•o  -  j  2*771nT 

=  2  v(nT)  e  (11-93) 

n=-« 

jZfPET 

Defining  the  complex  variable  z  =  e  ,  the  z-transform  of 

the  function  v(t)  is  defined  as 


«•  -n 

V(Z)  =  2  v(nT)  z  (11-94) 

n=-«» 

The  z-transform  is  the  frequency  domain  representation  for 
discrete-time  domain  signals  and  systems,  analogous  to  the 
Fourier  transform  for  continuous  time  domain  signals  and 
systems.  Difference  equations,  such  as  equation  11-92,  are 
changed  into  algebraic  equations  in  the  complex  z  domain 
through  the  use  of  the  z-transform  (Ref  28:511).  The  z- 
transform  is  defined  for  all  real  and  complex  functions  for 


which  the  series  defined  by  equation  11-94  converges,  there¬ 
fore  the  unique  z-transform  of  a  signal  function,  is  obtained 
by  applying  equation  11-94  and  specifying  the  region  of 
convergence.  As  with  the  Fourier  transform,  the  z-transform 
has  many  useful  properties  that  aid  in  analysis  and  design. 
(See  Appendix  C  page  C-7  for  a  list  or  Ref  29:Chapter  2  or  Ref 
43:Chapter  2.)  Although  the  discussion  thus  far  has  dealt 
with  digital  representation  of  signals,  the  development  also 
applies  to  systems. 

The  signals  discussed  in  this  section  were  assumed  to  be 
deterministic,  however  the  theory  also  applies  to  random  sig¬ 
nals.  As  mentioned  previously  the  statistical  properties 
developed  earlier  for  continuous  signals  applies  equally  well 
to  discrete  time  signals  by  changing  the  integrals  to  summa¬ 
tions.  Although  the  probability  theory  for  discrete  time 
signals  does  not  need  to  be  redeveloped,  the  effect  of  sam¬ 
pling  a  stochastic  signal  must  be  examined. 

Given  a  WSS  random  signal,  n(t),  with  zero  mean  and 

autocorrelation  R  (T) ,  the  statistics  of  the  sampled  version, 

n 

n(mT),  need  to  be  determined.  Assuming  independence  from 
sample  time  interval  to  sampling  interval,  E[n(mT)]=0.  To 
compute  the  autocorrelation  of  the  sampled  process,  we  assume 
flat  top  sampling.  The  autocorrelation  of  the  sampled  n(t)  is 

( i+1 )T  ( i  +  1 )T 

R  (T)  =  E[[1/T  <  n(t  )dt  ][1/T  $  n(t  )dt  ]] 


2  ( i+1  )T 

1/T  ^  R  (t  ,t  )  dt  dt 

IT  n  1  2  12 


(H-95) 


For  the  WSS  process,  equation  11-95  can  be  converted  from  a 
double  integral  to  a  single  integral,  thus 


2  (i+1  )T 

R  (f)  =  1/T  $  R  (T)  d*T 

n(mT)  iT  n 


=  R  (T)  /  T 
n 


(11-96) 


In  the  limit,  as  T  -*  0,  we  have  impulse  sampling  and  thus 


R  (T)  =  R  (T) 
n(mT)  n 


(11-97) 


For  impulse  sampling,  the  statistics  of  the  sampled  signal  are 
the  same  as  the  continuous  signal. 

The  z-transform  just  developed  was  for  infinite  length 


sequences , 


For  the  case  in  which  the  sequence  to  be 


represented  is  of  finite  duration,  the  discrete  Fourier 
transform  (DFT)  is  used.  The  DFT  is  a  sequence  rather  than  a 
continuous  function,  and  corresponds  to  samples,  equally 
spaced  in  frequency,  of  the  Fourier  transform.  The 
development  and  interpretation  of  the  DFT  representation  of  a 
finite  duration  sequence  is  based  upon  the  relationship 
between  the  finite-length  and  periodic  sequences  (Ref  43:87). 

The  complex  Fourier  series  representation  of  periodic 
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signal  x(t),  is 


x(t)  =  1/t  2  X(k)  e 

p  n=-<» 


j2'T7*t 


(11-98) 


where 


t  /  2 

pr  -jk2TIft 

X(k)  =  )  x(t)  e  dt 

-t  /2 
P 


(11-99) 


and  t  is  the  observation  period  as  defined  in  Figure  II-3. 
P 

If  x(t)  is  sampled  every  T  seconds,  such  that  the  Nyquist 
criterion  is  satisfied,  equation  11-99  becomes 


X(k)  =  2  x(nT)  exp[-j2TTknT/t  ] 

-N/2  p 


(11-100) 


where  N  is  the  total  number  of  points  in  the  sequence 

equalling  t  /T.  Substituting  1/N=T/t  into  equation  11-99, 
P  P 

and  also  noting  that  X(k)  and  x(n)  are  symmetric  for  positive 
and  negative  values  of  N,  the  DFT  of  x(nT)  is  defined 


X(k)  = 


x(nT)  e 


-  j  2TTkn/N 


N-1  kn 

=  2  x(nT)  W 

n=0 


(11-101 ) 


where  W  =  e 


.  The  inverse  DFT  is  thus 


-jZTT/N 


N-1  -kn 

x(nT)  =  1/N  2  X(k)  W 

k=0 


(11-102) 


1  ■*.' 


■  v* 

•  *  m 
»  ■  * 
vv, 

v'. 


The  variable  k  in  equations  11-101  and  11-102  represents 
the  digital  frequency  index.  It  is  related  to  the  analog 
frequency  variable  f,  by 


f  =  k  /  NT 


(11-103) 


The  periodic  sequence  X(k)  can  be  interpreted  as  samples, 
equally  spaced  in  angle,  on  the  unit  circle  of  the  z- 
transform,  for  one  period  of  the  periodic  sequence  x(nj.  Thus 
the  relationship  between  the  z-transform  and  the  DFT  is 


X(k)  =  X(Z ) 


z=  exp[j277k/N]  =  W 


-k 


(11-104) 


(Ref  43:89-90).  Because  the  DFT  is  a  modified  z-transform, 
the  DFT  properties  are  essentially  the  same  as  those  for  the 
z-transform. 

Because  the  radar  signals  and  systems  can  be  represented 
by  finite  sequences,  the  DFT  finds  extensive  usage  for  radar 
signal  design  and  analysis.  Although  these  radar  sequences 
are  generally  aperiodic,  by  zero  padding  the  sequences  can  be 
made  to  appear  (for  computational  purposes)  periodic  with  very 
long  periods. 
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Since  the  digital  signals  described  represent  continuous 
analog  signals,  it  is  quite  natural  that  the  digital  analysis 
of  these  signals  be  based  upon  the  correlation  principles 
previously  established  for  analog  signals.  The  correlation 
principle  is  implemented  by  a  digital  filter.  The  digital 
filter,  assumed  to  be  linear,  causal,  and  time  invarient,  is 
usually  an  approximation  of  a  desired  analog  filter. 

Digital  filtering,  like  analog  filtering,  can  be 
developed  in  either  the  time  domain  or  frequency  domain.  In 
the  discrete  time  domain,  digital  filtering  is  accomplished  by 
direct  convolution 

y (nT)  =  x(nT)  *  h(nT) 

oO 

=  2  x(kT)  h(nT  -  kT)  (11-105) 

k=-°° 

oo 

=  2  x(nT  -  kT)  h (nT) 

k=-  «d 

In  the  frequency  domain,  filtering  is  accomplished  either  by 
linear  recursive  or  nonrecursive  equations  or  by  FFT  (fast 
Fourier  transform)  operations.  Because  recursive  equations 
cannot  be  implemented  via  FFT  operations  and  also  due  to  the 
transient  period  needed  for  the  recursive  filter  to  stabilize, 
this  class  of  filters  will  not  be  discussed  for  use  in  clutter 
cancellation  schemes.  This  is  not  to  say  that  recursive 
filters  do  not  find  use  in  radar  processing,  for  recursive 


filters  digitally  implement  the  common  Butterworth,  Chebyshev 
and  elliptical  analog  filters,  however  they  will  not  be  used 
in  this  study. 

Nonrecursive,  or  finite  impulse  response  (FIR),  filters 
are  of  the  form 


a 


N-1  -n 

H( Z )  =  2  h(nT)  z  (11-106) 

n=0 

indicating  that  nonrecursive  filters  are  formed  using  input 
sequence  values.  The  convolution  of  a  signal  x(t)  and  filter 
response  h(t)  can  be  accomplished  in  the  frequency  domain  by 
multiplication  of  the  z-transform  of  the  FIR  filter  described 
in  equaton  11-106  with  the  signal  z-transform.  Thus 

Y ( Z )  =  X(Z)  H(Z)  (11-107) 

If  the  signal  and  filter  sequences  are  represented  by  DFT's, 
the  convolution  is  also  accomplished  by  multiplication  of  the 
two  DFT's.  Due  to  the  periodic  nature  of  the  DFT  calculation, 
if  the  sequences  are  aperiodic,  the  two  sequences  need  to  be 
zero  padded  so  as  to  appear  periodic  as  described  before. 
Both  sequences  should  be  of  period  L,  where  L  is  the  smallest 
power  of  two  greater  than  M  +  N,  with  M  and  N  being  the  number 
of  discrete  points  of  the  two  sequences  (Ref  61:379). 

The  most  straightforward  approach  to  obtain  a 
nonrecursive  filter  is  to  truncate  an  infinite  duration 
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impulse  response  sequence.  Given  the  desired  filter  response 
H{ jw)  then 

-jwnT 

H(jw)  =  2  h(nT)  e  (11-108) 

n=- 


where 


h(nT)  =  t  /27T 
P 


TT/t 


P 


TT/t 


P 


H(  jw) 


jwnT 

e  dw 


(11-109) 


Using  the  inverse  DFT,  equation  11-102,  a  zero-order 
approximation  to  h(nT)  is  obtained 


M  jTTTkn/M 

h(nT)  =  1 /2k  2  H(k)  e  -N  <  n  <  N  (11-110) 

k=-M  ~  “ 


where  M  is  generally  larger  than  N. 

Truncation  of  the  ideal  impulse  response  causes  the  FIR 
filter  to  experience  sidelobes.  These  filter  sidelobes  appear 
to  a  digital  clutter  canceler  system  as  clutter,  thus 
degrading  performance.  By  tapering  the  truncation  window,  the 
filter  sidelobes  can  be  reduced  from  -21  db  to  -74  db 
depending  on  the  tapering  scheme  (Ref  43:250).  Table  II-1 
illustrates  various  common  tapering  methods  and  their  asso¬ 


ciated  characteristics. 


Table  II-1 


Tapering  Methods 


Typ< 

Lag  window  (k{r)) 

Spectral  window  (k\f)) 

Equivalent  bandwidth  (Bt) 
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where :  Tmmhm*  lolal  time  lag  for  window,  h  »  sampling  interval,  and  m  -  number  of  samples  in  lag  window. 


The  modified  filter  thus  becomes 

h  (nT)  =  h(nT)  w(nT)  0  <_  n  <_  N-1  (11-111) 

m 

At  the  beginning  of  the  digital  filtering  section  it  was 
stated  that  nonrecursive  filters  can  be  implemented  via  a  FFT 
operation.  The  fast  Fourier  transform  is  simply  a  method,  or 
algorithm,  used  to  efficiently  compute  the  discrete  Fourier 


transform  (DFT)  of  a  discrete  time  series  of  data  samples  (Ref 
17:313).  Equation  11-110  illustrates  the  fact  that  nonrecur¬ 
sive  filters  can  be  implemented  using  the  DFT  of  the  desired 
frequency  response,  thus  the  nonrecursive  filter  can  indeed  be 
implemented  via  a  FFT. 

The  unique  feature  of  FFT  algorithms  is  that  they  elimi¬ 
nate  the  redundancy  in  DFT  calculations  by  reducing  the  DFT  to 
shorter  and  smaller  DFT  sequences.  Because  of  the  binary 
operations  of  computers,  FFT  algorithms  are  most  efficient  if 
N  is  a  power  of  two.  Reduction  in  calculations  are  achieved 
by  decimating  the  DFT  in  either  time  or  frequency.  The 
elegance  of  the  DFT  reduction  techniques  and  the  variety  of 
methods  available  to  accomplish  the  FFT  calculations  can  be 
found  in  Ref  15,  Ref  29  Chapter  6  and  Ref  43  Chapter  6. 
Equations  11-112  and  11-113  summarize  the  FFT  calculations 
after  the  DFT  has  been  time  and  frequency  decimated.  Given 
the  sequence  f(nT),  the  time  decimation  FFT  is 

N-1  2  nk  k  N-1  2  nk 

F(k)  =  2  g ( nT ) ( W  )  +  W  2  h ( nT ) ( W  )  (11-112) 

n=0  n=0 

where  g(nT)  =  f(2nT)  and  h(nT)  =  f((2n+1)T).  The  frequency 
decimation  FFT  of  sequence  f(nT)  is 

N/2-1  2  nk 

F(2k)  =  2  [ g(nT)  +  h(nT)](W  ) 

n=0 

(11-113) 

N/2-1  n  2  nk 

F( 2k+1 )  =  2  [[g(nT)  -  h(nT) ]W  ] (W  ) 

n=0 
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By  repeated  application  of  the  FFT  algorithms  described  by 
equations  11-112  and  11-113,  the  filtering  operations 
discussed  for  continuous  signals  can  be  accomplished.  The 
convolution  between  a  signal  and  a  filter  can  be  performed  by 
fast  Fourier  transforming  both  the  signal  and  filter, 
multiplying  the  two  FFT's  and  then  inverse  Fourier 
transforming  the  result.  The  inverse  FFT  is  accomplished  by 
appropriate  modification  of  the  arguments  of  either  11-112  or 
11-113.  The  FFT  operation  thus  provides  a  means  of 
accomplishing  the  signal  analysis  techniques  described  for 
continuous  radar  signals. 


Ambiguity  Function  Development 

Although  the  detectability  of  a  target  in  the  presence 
of  interference,  whether  noise  or  clutter,  is  dependent  on  the 
signal  strength,  the  ambiguity  removal  and  resolution 
accuracy  between  multiple  targets  is  a  function  of  signal 
parameters,  most  notably  the  waveform  duration  and  bandwidth. 
A  high  clutter  environment,  such  as  the  one  present  in  the 
spaceborne  problem,  can  be  thought  of  as  a  collection  of 
multiple  point  scatterers.  The  detection  of  a  target  in  the 
presence  of  clutter  is  thus  heavily  influenced  by  waveform 
selection.  Because  the  reflections  of  point  targets  vary  in 


range  and  doppler  from  the  target,  a  response  is  induced  in 
the  matched  filter  radar  receiver,  matched  to  a  specific  delay 
and  doppler,  that  is  no  longer  the  autocorrelation  function  of 


the  transmitted  signal.  The  response  of  secondary  targets, 
clutter,  at  a  slightly  different  range  and  velocity  may  appear 
at  the  matched  filter  output  when  the  desired  target  response 
is  at  its  peak.  As  a  result  of  signal  overlap,  due  to 
multiple  targets,  a  special  set  of  mathematical  functions 
called  time-frequency  autocorrelation  functions,  or  ambiguity 
functions  evolved  (Ref  41:283).  The  ambiguity  function 
describes  the  matched  filter  response  to  a  complex  signal 
envelope  as  a  function  of  radar  target  radial  velocity  and 
range.  This  function  is  determined  solely  on  the  transmitted 
waveform  and  receiver  filter  characterisitcs  (Ref  57:3-10). 
It  provides  an  indication  of  the  ability  of  a  waveform-filter 


combination 

to  resolve  multiple 

radar  reflectors 

at 

arbitrarily 

different  range  and 

velocities. 

The 

name 

"ambiguity 

function"  is  somewhat 

misleading 

in  that 

the 

function  describes  more  than  just  waveform  ambiguities. 

Woodward  used  the  word  to  show  that  the  total  volume  under 

2 

this  function  is  a  constant,  equal  to  (2E)  ,  independent  of 
the  shape  of  the  waveform  (Ref  56:42). 

To  obtain  the  time-frequency  autocorrelation  function, 
first  the  range  ambiguity  function  and  delay  resolution 
constant  will  be  determined,  then  the  doppler  ambiguity 
function  and  velocity  resolution  constant  will  be  determined, 
and  finally  the  combined  range  and  velocity  ambiguity  function 
will  be  developed. 

Using  the  complex  signal  notation  developed  earlier,  the 
time  domain  transmitted  function  is 


s(t)  =  Re  [u(t)  exp[j27Tf  t]]  =  Re  [s  (t)]  (11-114) 

o  a 


Given  two  stationary  point  targets  which  differ  only  in  range, 

the  received  signal  from  target  one  is  s  (t-t  )  and  from 

a  o 

target  two  s  (t-t  -T)  where  t  is  the  round  trip  delay  time 
a  o  o 

from  target  one  and  is  the  difference  in  time  delay  from 
target  one  and  two.  To  resolve  the  range  difference,  or  T,  we 
integrate  the  square  of  the  difference  between  the  two  signals 
over  all  time 


2  2 
5  =  | s  (t)  -  s  (t-T)  I  dt 

/ 


(11-115) 


Since  we  are  only  interested  in  the  time  difference  we  set  t 

2  o 

=  0  with  no  loss  in  generality.  Expanding  e  we  obtain 


2  2 
=  | s  ( t ) |  dt  +  |  s  (t-T)l  dt 

)  a  )  a 

’  *  * 
s  (t)s  (t-T)  dt  -  s  ( t ) s  (t-T)  dt 
)  a  a  J  a 


(11-116) 


a  a 


Because  of  the  periodicity  of  s  (t)  we  observe  that 


J|s  (t)|2dt  =  |s  (t-T)|2dt  =  |  |  u(  t )  | 


dt  (11-117) 


Multiplying  s  (t)s  (t-T)  and  s  (t)s  (t-T)  out  and  summing  the 
a  a  a  a 

products  we  obtain 
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s  (t)s  (t-T)  +  s  ( t ) s  (t-T)  = 
a  a  a  a 


2  Re  [s  (t)s  (t-T)] 
a  a 


(11-118) 


Combining  equation  11-117  and  11-118  into  equation  11-116  we 
get 


=  2 


u( t)  dt  - 


2  Re[exp[-j2’fTf  t] 
o 


j  u  (t)u(t- 


T)  dt]  (11-119) 


The  measure  e  shows  the  separation  between  s  (t)  and  s  (t-T) 

2  a  a 

therefore  if  Tjl  fl,  e  should  be  made  as  large  as  possible. 

The  first  term  in  equation  11-119  is  constant  and  equals  4E  . 

s 

The  second  term  flucuates  because  exp[-j27Tf  t],  which  implies 

o  2 

a  high  degree  of  range  ambiguity.  To  maximize  e  the  magni¬ 
tude  of  the  second  term  needs  to  be  minimized.  The  magnitude 
of  the  second  term  is  defined  as  the  range  ambiguity  function 


|c(T)|  =  u  (t)  u(t-T)  dt 


(11-120) 


We  see  that  |  c(T)  J  is  the  envelope  of  the  output  of  a  filter 
matched  to  the  signal  s(t).  The  matched  filter  output  forms 
the  crosscorrelation  between  the  received  signal  corrupted  by 
Interference  and  a  replica  of  the  transmitted  signal.  The 
replica  of  the  transmitted  signal  is  built  into  the  filter 
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through  the  frequency  response  of  the  filter  (equation  11-71). 
If  the  input  signal  is  the  same  as  that  for  which  the  filter 
was  designed  (the  filters  can  be  designed  for  specific  time 
delays  and  doppler)  the  output  is  the  autocorrelation  function 
(Ref  47 jXVII-5) . 

A  delay  constant,  T  ,  is  defined  as  the  total  area 

i  i2  r 

under  the  curve  c(T)  normalized  to  one  at  T=  0.  The 

constant  is 


(11-121 ) 


where  U(f)  is  the  Fourier  transform  of  u(t).  The  inverse  of 
the  delay  constant  is  defined  as  the  effective  bandwidth. 


B  =  |  f|u(f )|  2  df|2 

e  - L_  (11-122) 

'  4 

|u(f)J  df 

/ 


A  range  resolution  constant  is  also  defined  as 


ar  =  cT  /  2 
r 


(11-123) 


Thus  to  improve  range  resolution  (i.e 
effective  bandwidth  must  be  increased. 


reduce  ar)  the 


Similar  to  the  development  of  the  range  ambiguity 
function ,  the  mean  square  difference  is  used  to  measure  the 
velocity  difference  between  two  targets.  Working  in  the 
frequency  domain 

2  2  * 

e  =2  (|s  <f)|  -  s  (f)s  (f-f  )]  df  (11-124) 

f  J  a  a  a  f 

2 

Again  to  minimize  e  the  second  term  in  equation  11-124  must 

f 

be  minimized.  The  doppler  ambiguity  function  is  thus  defined 


2  &  v  /  <  F 

r  r 


(11-127) 


To  obtain  the  combined  range  and  velocity  resolution 


function  the  same  transmitted  waveform  as  in  the  two  previous 
cases  is  used.  We  assume  the  form 

s  (t)  =  u(t)  exp[ j2T7f  t]  (11-128) 

a  o 

The  time  delayed  and  doppler  shifted  signal  is 


s  (t-T)  =  u(t-T)  exp[j2T7U  -f  )(t-T)] 
a  o  d 


(11-129) 


Using  the  mean  square  difference  as  before 


2  f  ,  2 

e  =  ||  s  (t)  -  s  (t-T)!  dt 
/  a  a 


(11-130) 


Expansion  of  the  integrand  yields 


e  =  2  I  I  u ( t ) I  dt  -  2  Re  [  exp  [  j  2’7T(  f  -f  )T] 

o  d 

t  * 

u(t)u  (t-T)  explj2'T7t  t]  dt 
J  d 


(11-131  ) 


The  combined  range  and  velocity  ambiguity  function  is  then 


|x(T,f  ) 
d 


|  =  Ju(t)u 


(t-T)  exp[ j 277f  t]  dt  (11-132) 

d 


Thus  we  see  that  the  ambiguity  function  describes  the  complex 
signal  envelope  at  the  matched  filter  receiver  output  as  a 
function  of  target  range  and  radial  velocity. 


A  combined  resolution  constant  can  be  obtained  as  in  the 


doppler  and  range  developments 


f  2 

|x<T,f  )|  d-Tdf 

LI _ ±- - & 

. 2 

|x(0,0) 


(11-133) 


Assuming  | X ( 0 , 0 ) |  is  normalized  to  one,  the  delay  and  doppler 
resolution  constants  can  be  written  respectively  as 


T  =  1/B  =  |  X(T, 0)|  dT 

re/ 

F  =  1/T  =||  X( 0 , f  )|  df 

r  e  /  d  1  d 

The  plot  of  the  square  magnitude  of  the  ambiguity 

function  is  called  the  ambiguity  diagram.  This  ambiguity 

surface,  in  the  T  and  f  plane,  is  used  as  a  design  tool  to 

d 

indicate  the  resolution  and  ambiguity  of  a  waveform-filter 
combination.  Ideally  the  ambiguity  diagram  should  consist  of 
a  single  spike  of  infinitesimal  thichness  at  the  origin  and 
zero  everywhere  else.  The  single  peak  eliminates  ambiguities, 
and  its  infinitesimal  thichness  at  the  origin  permits 
simultaneous  determination  of  frequency  and  echo  time  delay 
with  a  high  degree  of  accuracy.  The  narrow  spike  permits  the 
resolution  of  any  two  targets  no  matter  how  close  they  are 
together.  In  reality  the  ideal  ambiguity  diagram  does  not 
exist  and  depending  on  the  specific  application  a  suboptimal 
waveform  can  be  chosen.  The  width  of  the  ambiguity  diagram 
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(11-134) 

(11-135) 


central  peak  indicates  the  resolution  ability  of  multiple 
targets  in  delay  and  doppler.  The  sidelobes,  or  secondary 
peaks  provide  an  indication  of  the  amount  of  self-induced 
clutter  and  thus  the  masking  of  targets  by  mutual  interference 
(Ref  49:25).  For  a  given  waveform-filter  pair,  the  clutter 
power  out  of  the  filter,  designed  for  a  specific  delay  and 
doppler,  is 


=  ||p(T/fd)  |xC^fd>|2  dTdf 


(11-136) 


where  p(T, f  )  is  the  distribution  of  clutter  scattering  points 
d 

in  delay  and  doppler.  Equation  11-136  illustrates  the  fact 

that  the  clutter  distribution  as  well  as  the  waveform 

influences  the  amount  of  clutter  power  return.  Waveforms  hav- 
2 

ing  |  XC^,  f  )|  — >0  in  the  *T,  f  region  where  clutter  exists,  can 
d  d 

be  designed  to  reduce  the  amount  of  induced  clutter. 
Ambiguity  functions  used  in  various  applications  are  listed  in 
Ref  57:3-15  -  3-30. 

Having  gained  an  understanding  of  how  to  represent  and 
analyze  signals  in  both  time  and  frequency  domain,  and  also 
having  gained  insight  into  methods  of  implementing  the  signal 
analysis  techniques,  clutter  and  clutter  cancellation 
techniques  can  now  easily  be  developed  and  described. 


Ill  Clutter  Models  and  Cancellation  Techniques 

Introduction 

The  fundamental  problem  of  radar  design  is  the  extrac¬ 
tion  of  target  range  and  velocity  information  from  the  return 
radar  signal  that  is  corrupted  by  clutter  and  noise  energy. 
Target  information  is  contained  in  the  radar  echo  signal's 
amplitude,  frequency  and  phase.  To  extract  the  target  infor¬ 
mation  buried  in  the  echo  signal  the  signal  must  be  processed, 
which  is  essentially  accomplished  via  comparison  techniques, 
i.e.  comparing  the  echo  signal  with  the  transmitted  signal 
and/or  comparing  the  echo  signals  with  other  echo  signals 
which  are  varied  in  time  or  space.  Depending  on  the  type  of 
information  desired  (mere  indication  of  a  target's  presence, 
target  range,  relative  elevation  or  azimuth  angles,  target 
velocity  or  target's  physical  characteristics)  the  echo  signal 
is  processed  differently. 

Effective  extraction  of  range  and  velocity  information 
from  the  echo  signal  is  accomplished  by  first  cancelling 
clutter.  To  effectively  develop  clutter  cancellation  tech¬ 
niques,  it  is  necessary  to  first  gain  an  understanding  of  the 
nature  of  clutter.  Using  the  fundamentals  that  were  esta¬ 
blished  in  Chapter  II  we  will  in  this  chapter  first  describe 
clutter,  then  briefly  explain  a  typical  MTI  radar,  and  con¬ 
clude  with  a  description  of  various  clutter  cancellation  tech¬ 
niques  available  to  be  employed  on  spaceborne  radar  systems. 


Elements  of  Clutter 


Radar  clutter  is  the  unwanted  return  or  group  of  returns 

of  reflected  radar  energy  which  results  from  distibuted 

scattering  elements  such  as  buildings,  chaff,  birds,  insects, 

precipitation,  land  and/or  sea.  The  distributed  scattering 

centers  produce  individual  re-radiated  electrical  fields  which 

are  of  vector-phasor  nature  having  amplitude  A  and  phase  4>  . 

i  i 

The  individual  vectors  (having  both  horizontal  and  vertical 
polarization)  of  all  the  illuminated  scatterers  may  be  summed 
to  produce  the  unwanted  electrical  clutter  field  that  is 
received  by  the  radar.  There  are  several  models  which  attempt 
to  explain  the  scattering  of  the  transmitted  electromagnetic 
waves  produced  by  non-uniform  surfaces,  (Refs  19,  33,  35-37, 

40,  55,  68)  however  they  are  not  sufficiently  well  developed 

to  be  associated  with  every  type  of  terrain  of  interest  (Refs 
19:9,57:25-4).  Attempts  have  also  been  made  to  characterize 
terrains  probabalistically  and  then  determine  the  return  echo 
field.  These  models  have  met  with  some  success  (Ref  36:124, 
57:25-4).  While  the  exact  interaction  of  the  electromagnetic 
wave  with  various  terrains  is  not  known,  experimental  results 
show  that  ground,  sea  and  weather  clutter  are  dependent  on:  1 ) 
radar  parameters  such  as  wavelength,  power,  polarization, 
pulse  width,  antenna  scan  rate,  frequency  agility,  and  resolu¬ 
tion  cell;  2)  geometry  parameters  such  as  grazing  angle, 
direction  of  illumination  and  range;  3)  environmental  condi¬ 
tions  such  as  wind  velocity,  moisture  content  and  season  of 
the  year  and;  4)  type  of  clutter,  described  by  surface  rough- 


ness,  number  and  orientation  of  scatterers  in  resolution  cell, 
roughness  of  surface  depth,  surface  conductivity  and 
permittivity  (Refs  47:XV-9,  55:63,  57:25-4).  The  fundamental 

relationship  between  all  of  the  above  listed  parameters  and 
the  clutter  returns  are  graphed  in  Skolnik's  Radar  Handbook 
(Ref  56)  and  Long's  book  Radar  Reflectivity  of  Land  and  Sea 
(Ref  35).  Since  there  are  a  large  number  of  parameters  which 
affect  clutter  returns,  and  the  interactive  mechanisms  between 
the  transmitted  radar  signal  and  the  scatterers  has  not  been 
firmly  established,  some  manageable  way  of  describing  clutter 
must  be  obtained. 

The  clutter  signal  is  of  the  form 

N  j^(t) 

c(t)  *  2  A(t)  s(t-f)  e  (III-1) 

i 

where  A(t)  is  the  amplitude,  is  the  time  delay,  tf(t)  is  the 
phase  and  N  is  the  number  of  scattering  centers.  While  these 
quantaties  may  at  times  be  considered  deterministic,  they  are 
generally  described  as  random  variables.  Clutter  can 
therefore,  under  some  conditions,  be  thought  of  as  quasirandom 
noise.  Experimental  data,  however,  indicates  that  fluctua¬ 
tions  in  clutter  returns,  due  to  environmental  or  geometrical 
conditions,  cause  clutter  to  be  more  frequently  considered  as 
a  nonstationary  process.  Clutter  returns  exhibit  time 
dependence  within  resolution  cells  and  spatial  dependence  from 
resolution  cell  to  resolution  cell  (Ref  47:XV-8).  Thus  clut¬ 
ter  can  be  described  by  not  only  average  returns,  but  also  by 


statistical  descriptions  of  temporal  and  spatial  variations. 

For  a  given  set  of  radar,  geometry  and  environmental 
conditions,  clutter  is  usually  described  by  any  one  or  all 
three  of  the  following  methods:  1)  an  average  radar  cross 
section  per  unit  area,  2)  an  amplitude  probability 
distribution  and  3)  a  power  spectral  density. 

Since  clutter  comes  from  extended  areas  containing  many 
scattering  elements,  clutter  can  be  described  statistically  in 
terms  of  the  average  radar  backscatter  cross  section  per  unit 
area,  denoted  a4*  ,  or  clutter  reflectivity  factor.  Clutter 
reflectivity  is  a  measure  of  the  reflective  ability  of  an  area 
of  the  earth.  Theoretically,  clutter  reflectivity  is  given  by 

o  2  2 

a*  =  [lim  R  fl)  4T7fc  IE  /E  1  ]  /A  (III-2) 

r  i  c 

where  E  and  E  are  the  reflected  and  incident  field  strength 
r  i 

respectively  and  A  is  the  area  illuminated  by  the  radar  beam. 

c 

(Appendix  B  derives  the  area  calculations.)  Because  of  the 
spactial  and  temporal  variations  of  the  reflecting  earth,  the 
backscatter  coeffcient  is  generally  calculated  experimentally 

O 

by  applying  the  radar  equation  and  solving  for  a4  given  the 
observed  power  returns.  The  clutter  reflectivity  is  then 
calculated  by  solving  equation  (III-3) 

a*  =  o*  *  (A  )  ( III-3 ) 

c  c 

The  clutter  reflectivity  factor  is  used  in  signal  to  clutter 
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calculations  similar  to  those  done  in  noise  calculations.  The 
signal  to  clutter  ratio  is 

S/C  =  o-1  /  cr*0  A  ( III-4 ) 
t  c 

In  many  cases  the  clutter  power  is  far  greater  than  the  signal 
power  (i.e.  when  average  backscatter  reflectivity  is  large  or 
the  area  of  illumunation  is  large)  thus  masking  the  presence 
of  targets.  Clutter  cancellation  techniques  discussed  later 
will  show  how  the  target  signal  can  be  extracted  from  the 
clutter  when  the  S/C  ratio  is  less  than  one. 

Amplitude  probability  distributions  are  used  to  describe 
clutter  so  as  to  be  able  to  predict  the  false  alarm 
probabilities  due  to  clutter  (Ref  47:XV-11).  Knowledge  of  the 
amplitude  statistics  is  important  in  designing  constant  false 
alarm  receivers,  CFAR.  Because  the  backscatter  from  clutter 
shows  both  temporal  and  spatial  variation,  a  time  varying 
probability  distribution  or  random  process  is  really  needed  to 
describe  the  clutter  return  statistics.  The  probability 
models  which  have  found  widest  use  in  predicting  clutter 
fluctuation  statistics  assume  stationarity  and  independence  of 
scatterers  within  a  resolution  cell  (Ref  47:XV-11). 

The  simplest  probability  model  is  the  Rayleigh  or 
Gaussian  model,  which  assumes  a  uniform  spatial  distribution 
of  independent  scatterers  with  equal  amplitude  and  relatively 
uniform  phase  distribution  from  0-2*7T  (Ref  47:XV-14).  The 
Gaussian  envelope  clutter  model  is  used  for  land  clutter  with 


2  ft 


vegetation  that  consists  if  non-dominant  tree  trunks  or  for 


sea  clutter  when  large  individual  waves  are  not  numerous  at 
grazing  angles  greater  than  five  degrees  and  pulse  widths 
greater  than  .5  microseconds  (Ref  52:38).  The  predetection 
clutter  process  in  terms  of  its  voltage  envelope  v(t)  and 
phase  Q(t)  is 


c  (t)  =  v(t)  cos[w  t  -  Q(t)]. 
G  c 


(III-5) 


Expanding  the  cosine  term,  c(t)  can  be  represented  by  its 
quadrature  components  as 


c  (t)  =  x(t)  cos[w  t]  -  y(t)  sin[w  t] 
G  c  c 


( III-6 ) 


where  x(t)  and  y(t)  are  zero  mean,  identical  distibuted  low 

2 

pass,  independent  normal  processes  with  variance  cr*  and  w  is 

c 

the  carrier  frequency.  The  voltage  envelope  of  c(t)  is 


/ 2  2 
v(t)  = V  x  (t)  +  y  (t) 


( III-7 ) 


and  the  phase  of  c(t)  is 


0(t)  =  tan  [y(t)/x(t)J 


( III-8 ) 


The  joint  probability  density  funtion  of  x(t)  and  y(t)  is 


q(x,y)  =  q(x)q(y)  = 

2  2 
y  )  /2cr*  ] . 


2  2 
( 1  /277bJ  )  exp[  -(x  + 


(HI-9) 


V, 

>1 


To  find  the  density  in  terms  of  v(t)  and  0(t)  we  need  to  make 
a  coordinate  transformation  and  express  the  joint  probability 
density  function  of  x(t)  and  y(t)  in  the  v,  ©  coordinate 
system.  The  joint  density  function,  q(x,y)  is  thus 

q{x,y)dxdy  =  q(v  cos©,  v  sin©)  vdvd©  (III-10) 

since  the  element  of  the  area  dxdy  in  the  x,y  plane 
corresponds  to  the  element  of  area  vdvd©.  We  state  that 

q(x,y)dxdy  =  p(v,©)dvd©  (III-11) 

and  then  see  that 

p  ( v ,  0 ) dvd©  =  q( v  cos©,  v  sin©)vdvd©.  (III-12) 

Making  the  appropriate  subsitutions  into  equation  (III-8)  we 
obtain 


2  2  2 

p(v,0)  =  (v/2T7b')  exp[-(v  cos  ©  + 

2  2  2 
v  sin  ©) /2 cr*  ] 


2  2  2 

p(v,©)  =  (v/2{^o,)  expt-v  /2oJ]  v  >.  0 

=0  v  <  0. 


(HI-13) 


Equation  III-13  implies  that  ©  is  uniformly  distributed  and 


that  the  marginal  density  function  of  the  phase  is 


p<0)  =  1  /  2*7T  0  1  6  2.  27T  (III-14) 

The  marginal  density  for  the  envelope,  v(t),  can  thus  be  found 
by  integrating  over  p{v,9)  with  respect  to  0.  Therefore  we 
have 


2*T T  2  2  2 

p(v)  =  j  (v/2TJcrJ  )  exp[-v  f2a*]  d© 

0 

2  2  2 

p(v)  =  (v  /  cr* )  exp[-v  / 2cr* ]  v  _>  0 

=  0  v  <  0 


(III-15) 


(III-16) 


which  is  a  Rayleigh  probability  density  function  that  is 

independent  of  time  due  to  stationarity  and  is  only  a  function 

2 

of  clutter  power,  P  =  2a*  (Ref  52:38).  P  is  the  mean 

c  c 

square  value  of  the  envelope. 

The  power  spectral  density  for  the  Gaussian  model  is 

thus 


2  2 

S  (f)  =  (P  exp[-f  I2a*\  (HI-17) 

G  c  f  f 


where  P  is  the  clutter  power.  The  term  cr*  is  the  standard 
c  f 

deviation  of  the  spectrum  and  is  related  to  the  RMS  velocity 

clutter  spread  a*  by  a*  =  2CT*  /  \ . 

v  f  v 

When  the  clutter  process  c(t)  is  sampled  by  a  pulsed 


radar,  the  probability  density  of  the  sample  is  given  by  the 


multidimensional  Gaussian  distribution 


n  -.5  T  -1 

p(x  ...  x  )  =  [(2  ) IMI ]  exp[-.5x  M  x]  (III-18) 

In  x 


where  T  is  the  transpose  of  the  vector  x  and  the  covariance 
matrix  is 


2 

M  =  c r* 

x  x 


1  p  p 

12  13 

P 

21 


P 

In 

P 

2n 


(III-19) 


p  ...  1 

nl 


2 

and  p  =  exp[-(i-j)  2'TEr*T/2],  For  analytical  purposes  the 

ij  f 

clutter  process  is  assumed  to  be  a  Gaussian  Markov  process. 

The  Rician  model  is  similar  to  the  Gaussian  model  except 
that  a  predominant  steady  scatterer  (S)  is  added  to  the 
distributed  clutter.  The  random  clutter  process  is  given  by 


C  (t)  =  (S  -  x(t))  cos(w  t)  -  y(t)  sin(w  t)  (III-20) 
R  c  c 


with  voltage  amplitude  probability  density  function 


2  ^  2  2  „  2 
p(v)  =  (1  +  m  )/  ?  exp[-m  ]  exp[(-v  /P)(1  +  m 


I  (  2msj ( 1  + 


2  2 
m  )  v  /P) 


)] 

(III-21  ) 
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2  2 
where  m  is  the  ratio  of  the  steady  (S  )  to  distributed  (P  ) 

ws  2  o 

power  and  total  power  P  =  S  +  P  (Ref  52:41).  The  power 

o 

spectral  density  is 


2  2  2 
S  (f)  =  m  P  S(f)  +  P  exp[  (-f/2oJ  ) /\f2fTcri  ] 
R  o  o  f  f 


( III-22 ) 


The  log-normal  model  has  been  used  to  model  clutter 
obtained  from  high  resolution  (pulse  width  less  than  .5  micro¬ 
seconds)  sea  data,  where  the  sea  is  observed  at  grazing  angles 
less  than  five  degrees.  The  log-normal  model  has  also  been 
fitted  to  ground  clutter  data  observed  at  low  grazing  angles 
(Ref  52:42).  The  log-normal  random  clutter  process  is 


c  (t)  =  v(t)  cos(w  t  +  ©) 
LN  C 


( III-23 ) 


where  v(t)  is  the  voltage  envelope  as  in  equation  III-6  and  9 

is  a  uniformly  distributed  random  phase  variable.  If  the 

quadrature  components  are  zero  mean,  uncorrelated  .  and  have 

identically  band  limited  spectral  densities,  then  c  (t)  is 

LN 

wide  sense  stationary  and  the  voltage  amplitude  probability 
distribution  is  given  by 


_  2  2 

p(v)  =  1  / (7s Jrn&’w)  exp[(ln  v/v  )  /  2cr*  ] 


( III-24 ) 


a 

n 


where  v  _>.  0/  o'*  is  the  standard  deviation  of  the  normal 

v 

distribution,  and  v  is  the  median  value  of  the  distribution 

m 
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(Ref  52:45).  The  log-normal  model  characterizes  well  the 
clutter  from  specular  reflection  associated  with  wave  crests 
and  from  large  directive  scatterers  that,  depending  on  viewing 
aspect,  produce  bright  ridges  and  shadowing  (Ref  52:42).  In 
the  presence  of  sea  clutter,  the  Gaussian  pfd  tends  to  under¬ 
estimate  the  range  of  voltage  amplitudes  of  the  real  clutter 
while  the  log-normal  pdf  tends  to  overestimate  the  range  of 
values  (Ref  56:479). 

The  Weibull  clutter  model  has  properties  that  lie 
between  the  Gaussian  and  log-normal  clutter  models.  This 
model  has  the  potential  to  represent  real  clutter  over  a  wide 
range  of  conditions  (Ref  52:48).  It  is  used  to  model  both 
land,  sea,  and  atmospheric  clutter.  The  Weibull  clutter  model 
voltage  amplitude  out  of  a  clutter  envelope  detector  is 


<x  -  1 

p(v)  =  In  2  (v/v  )  exp[-ln  2 (v/v  )]  (III-25) 

m  m 


where  °<  is  a  parameter  that  relates  to  the  skewness  of  the 

distribution  and  v  is  the  median  value  of  the  distribution 

m 

(Ref  56:480). 

Several  other  distributions  have  been  applied  to 
empirical  clutter  data.  These  include  the  chi-square 
distribution,  the  K-distributions  and  the  IIT  model  (complex 
overland  airborne  radar  derived  model).  These  models  have 
not,  however,  had  widespread  use. 

The  power  spectral  density  which  describes  clutter 
returns  is  probably  the  most  useful  for  analyzing  clutter 


S'  /*/ V  - 
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cancellation  techniques.  The  most  widely  used  discription  of 


the  clutter  spectrum  is  that  proposed  by  Barlow.  Barlow, 
analyzing  experimental  data  made  from  ground  targets,  sea 
echoes,  and  rain  clouds,  proposed  a  Gaussian  model  for  the 
frequency  distribution  of  clutter  (Ref  5:351).  The  power 
spectrum  from  fixed  targets  is  thus  approximated  by 


2 

W(f)  =  W  exp[-a(f/f  )  ]  ( III-26 ) 

o  t 


where  W  is  the  mean  value  of  the  power  density,  f  is  the 
o  t 

radar  transmitting  frequency  and  "a"  depends  on  the  nature  of 

clutter  and  weather  conditions.  Values  of  "a"  range  from 
19  15 

3.9x10  for  sparsely  wooded  hills  on  a  calm  day  to  2.8x10 

for  rain  clouds  (Ref  47:XV-25).  Figure  III-1  show  the 

frequency  spectra  of  various  types  of  clutter. 


Curve  Target  a 

17 

1  Heavily  wooded  hills,  2.3x10 
20  mph  wind  blowing 

19 

2  Sparsely  wooded  hills,  3.9x10 
calm  day 

16 

3  Sea  echo,  windy  day  1.41x10 

15 

4  Rain  clouds  2.8x10 

16 

5  Window  "jamming"  1.0x10 


FREOUENCY  IN  GPS. 

(Ref  5:351  ) 

Fig.  III-1  Frequency  spectra  of  various  types  of  fixed  targets 


Equation  III-26  can  be  written  into  the  more  familiar  Gaussian 
form 

2  2 

W(f)  =  W  exp[-f  /2crJ]  ( III-27 ) 

o  c 

2  2 

where  a  =  f  /2cr*  and  cr1  is  the  standard  deviation  of  the 
t  c  c 

clutter  spectrum. 

While  the  frequency  spectrum  for  clutter  is  modeled  as 
Gaussian,  it  is  in  actuality  a  filtered  version  of  the  trans¬ 
mitted  signal  spectrum  envelope.  The  return  radar  signal 
resulting  from  clutter  can  be  thought  of  as 

2 

C(f)  =  1H( f ) 1  S ( F )  ( III-28 ) 

where  C(f)  is  the  clutter  power  spectrum,  S(f)  is  the  signal 
power  spectrum  and  H(f)  is  the  transfer  function  that  des¬ 
cribes  the  clutter  response.  Due  to  the  effects  of  antenna 
tapering,  smearing  by  the  convolution  of  the  IF  spectrum  with 
itself  (Ref  11:484),  system  instabilities,  movement  by  the 
radar  and  the  random  nature  of  the  clutter,  the  fine  grain 
structure  of  the  transmitted  signal  spectrum  is  generally  lost 
and  the  clutter  power  spectrum  can  be  approximated  by  the 
Gaussian  function  as  proposed  by  Barlow. 

The  spectrum  of  the  clutter  return  results  from  the 
convolution  of  thu  Gaussian  spectrum  with  the  transmitted 
frequency  and  all  harmonics.  An  example  of  a  CW  radar 
spectrum  is  shown  in  Figure  III-2a  while  the  spectrum  of  a 


pulse  doppler  radar  is  shown  in  Figure  III-2b. 

r1  ls(f)| 


a  b 
Fig.  III-2  Clutter  Spectrum 

Figure  III-2b  illustrates  why  it  is  important  for  a  doppler 
radar  to  use  a  high  pulse  repition  frequency  so  as  to  prevent 
aliasing  between  the  harmonics,  thereby  eliminating  velocity 
ambiguities. 

The  clutter  spectrum  shown  in  Figures  III-2a,2b  is 
spread  primarily  due  to  four  effects:  1)  antenna  scanning 
modulation;  2)  radar  platform  motion;  3)  transmitter  frequency 
and  4)  internal  clutter  motion. 

A  scanning  antenna  induces  a  modulation  that  is 
characterized  by  a  long  pulse.  Given  a  Gaussian  antenna 
pattern,  the  contribution  to  clutter  spread  due  to  antenna 
scanning  is 

CT*  =  «/  oJl'TTO  )  Hz  ( III-29 ) 

S  B 


where  oc  is  the  scan  rate,  and  ©  is  the  two  way  beamwidth 

B 


The  platform  motion  has  two  effects  on  the  clutter 
spectrum.  First  there  is  a  translation  of  center  frequencies 
due  to  the  doppler  shift  induced  by  the  platform.  This 
doppler  shift  is 

f  =  2  (v  /X)  cos^  cose  ( III-30 ) 

d  p 

where  O  is  the  azimuth  angle  between  the  platform  velocity  and 
the  direction  of  the  radar  beam  and  $5  is  the  depression  angle 
to  clutter  patch  (Ref  53:210).  Besides  the  shifting  in  the 
center  frequency  of  the  clutter,  the  spectrum  is  widened  due 
to  a  relative  rotation  induced  by  the  moving  observer  with 
respect  to  a  fixed  point.  The  platform  motion  contribution  to 
clutter  spread  is 

a*  =  ( . 6/D)  v  cos  <!  sin  0  Hz  (III-31) 

P  P 

where  v  is  the  radar  velocity,  the  angles  are  defined  as 
P 

above  and  D  is  the  diameter  of  the  antenna  (Ref  53:212). 

The  transmitter  frequency  instabilities  cause  a  spread¬ 
ing  of  the  clutter  spectrum  because  of  the  decorrelation  from 
pulse  to  pulse.  The  contribution  to  clutter  spreading  due  to 
transmitter  frequency  instability  is 

o'*  =  (.42  /  B)  f’(t)  Hz  ( III-32 ) 

f 

where  B  is  the  transmitter  pulse  bandwidth  and  f(t)  is  the 
transmitter  drift  rate  (Ref  53:216). 
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The  scatterers  often  move  relative  to  one  another 
because  of  wind,  thus  causing  a  spread  in  the  clutter 
spectrum.  The  standard  deviation  of  the  internal  motion  is 


cr*  =  2  V  /  X  Hz 
v  RMS 


(III-33) 


Typical  values  of  cr*  are  shown  in  Table  III-1 

v 


TABLE  III-1 

Characteristics  of  Clutter  Spectra 


Source  of 
clutter 


Wind 
speed 
(knots ) 


Sparse  woods  calm 
Rocky  terrain  1 0 
Wooded  hills  1 0 

ii  ii  nn 


Sea  echo 


ii 

II 

r- 

ii 

II 

8-20 

p  %m» 

ii 

II 

(windy) 

Chaff 

0-10 

% 

II 

25 

Rain  clouds 

ii  ii 


Ratio 

2 

m 


30 

5.2 


Barlow's 

a 


3.9x10 


7.2x10 


2.3x10 


9x10 


1  .1x10 

1 

2.4x10 

1 

( 1 -2 )x1 0 


( 0 . 6-2 . 6 )x1 0 

1 

1.4X10 

1 

0  ( 1 . 4-8 ) xl 0 


7x10 


0  ( 0 . 7-3 )x1 0 

1 

2.8x10 


c  v 

(cm/sec)  ( ft/ sec) 

3.5  0.057 


8  0.13 

45  0.78 

23  0.38 

65  1.06 

140  2.3 

165-205  2. 5-3. 3 

100-220  1.5-3. 5 

183  3.0 

75-180  1. 2-3.0 


250 

370-800 

410 


4.1 

6-13 

6.7 


(Ref  6:100) 


The  four  standard  deviations  combine  together  to  give  an 
overall  clutter  spread  of 


2  2  2  2  2 

o'1  =  oJ+cr,+oJ+oJ  (III-34) 

c  s  p  v  f 


2  2  2 
cr*  =  ( <*  /  K/T'Tro  ))  +  (  ( .  6v  /D)  cosji  sinO)  + 

c  B  p 

2  \  2 

( (  .  42/B )  f'(t))  +  (2  V  /X)  .  (III-35) 

RMS 


For  a  spaceborne  radar  system  the  effects  of  the  side- 
lobes  must  be  added  to  the  spreading  of  the  clutter  spectrum 
since  the  sidelobes  can  illuminate  a  large  clutter  area.  The 
sidelobes  in  a  spaceborne  radar  system  will  produce  clutter 
returns  which  extend  from  2v/\on  either  side  of  the  transmit¬ 
ter  frequency  (f  )  because  the  sidelobes  extend  in  different 

o 

directions  around  the  main  beam.  The  clutter  returns  result¬ 
ing  from  sidelobes  that  point  directly  below  the  spacecraft 

will  produce  a  strong  return  at  f  because  the  radial  velocity 

o 

at  that  point  is  zero.  A  typical  clutter  spectrum  of  a  space¬ 
borne  radar  system  is  shown  in  Figure  III-3. 


Frequency 


e 


(Ref  56:146) 


Fig.  III-3  Spaceborne  Radar  Clutter  Spectrum 


Two  other  spectral  models  have  found  use  to  describe 
clutter  returns  in  certain  situations.  These  models  are  the 
Markoff ian  model  and  the  cubic  model.  The  Markoff ian  model 
describes  the  frequency  spectrum  of  clutter  that  is 
bandlimited  white  noise  like.  The  model  is 


2  2  2 

W<f)  =  W  [f  /  (f  +  f  )] 
o  c  c 


f  >  0  ( III-36  ) 


where  W  is  the  mean  value  of  the  power  density  and  f  is  the 
o  c 

half  power  frequency  (Ref  47:XV-26).  The  cubic  model  has  been 

proposed  by  Fishbeir  for  clutter  returns  at  X-band.  The  model 

is 


3 

W(f)  =  1  /  [1  +  (f/f  )  ]  ( III-37 ) 

c 


/?v  -1 

where  f  =  k  e  ,  with  k  =  1.334,  .356  (knots)  and  v  is 

c 

the  windspeed  (knots)  (Ref  47:XV-27). 


Typical  Radar  Receiver 

The  clutter  cancellation  techniques  analyzed  in  the  next 
section  assume  a  pulsed  doppler  radar  system.  Shown  in  Figure 
III-4  is  a  generalized  block  diagram  of  a  coherent  pulsed 
amplifier. 
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(Ref  56:105) 

Fig.  III-4  Block  Diagram  of  a  MTI  Radar  Receiver 


Either  a  pulsed  oscillator  system  or  a  pulse  amplifier  system 
can  be  used  in  a  pulse  doppler  radar.  Since  the  operation  of 
the  pulsed  amplifier  is  very  similar  to  that  of  the  pulsed 
oscillator  and  for  a  spacecraft  application  a  pulsed  amplifier 
system  would  be  used,  only  the  pulsed  amplifier  operation  will 
be  described.  An  analysis  of  the  pulsed  oscillator  operation 


can  be  found  in  Ref  6:192-195. 

The  system  shown  in  Figure  III-4  is  one  form  of  a  co- 


herent  MTI  radar,  in  which  the  phase  reference  is  maintained 


in  the  radar  itself  (Ref  6:192).  Coherent  reference  is  main¬ 
tained  by  a  stable  oscillator  called  the  COHO,  or  coherent 
oscillator,  operating  at  IF.  Besides  providing  a  reference 

signal  from  pulse  to  pulse,  the  output  of  the  COHO,  f  ,  is 

c 

mixed  with  the  stable  local  oscillator,  STALO,  at  frequency, 
f  .  The  function  of  the  STALO  is  to  provide  the  necessary 
frequency  translation  from  the  IF  to  transmitted  RF  frequency. 
The  transmitted  signal  is  thus  a  combination  of  the  COHO  and 
STALO  frequencies. 

The  RF  echo  signal  is  heterodyned  with  the  STALO  signal 
to  produce  the  IF  signal.  The  reference  signal  from  the  COHO 
is  then  mixed  and  low  pass  filtered  with  the  IF  echo  signal  in 
the  phase  detector.  The  phase  detector  output  is  proportional 
to  any  phase  difference  between  the  COHO  and  IF  echo  signal 
(Ref  56:105).  The  zero  frequency  output  of  the  phase 
detector,  sometimes  called  the  video,  is  then  processed  to 
remove  clutter. 


Cancellation  Techniques 


Clutter  masses  can  be  considered  as  stationary  targets, 
thus,  relative  to  the  radar  antenna,  targets  and  clutter 
masses  exhibit  different  relative  speeds.  The  differences  of 
target  and  clutter  velocities,  conveyed  by  doppler  shifts  in 
the  transmitted  signal,  are  used  to  seperate  moving  targets 
from  undesired  clutter.  A  radar  that  uses  the  doppler 
frequency  shifts  as  a  way  of  distinguishing  moving  targets 


from  clutter  is  called  a  MTI  (moving  target  indication) 
system  (Ref  56:101).  A  MTI  system  essentially  employs  fil¬ 
tering  networks  (or  their  equivalent)  which  null  the  returns 
from  fixed  targets,  which  have  no  doppler  shift,  while  passing 
the  frequency  shifted  echoes  of  moving  targets.  Clutter  can¬ 
cellation  by  this  means  is  limited  by  the  extent  that  the 
energy  spectrum  from  the  clutter  mass  is  distinguishable  from 
that  of  the  target  (Ref  24:82).  As  illustrated  in  the  pre¬ 
vious  section,  the  clutter  power  spectrum  of  a  spaceborne 
radar  occupies  a  range  of  frequencies,  which  at  times  is  in 
the  vicinity  of  components  of  the  target  radial  velocity. 
This  makes  it  difficult  to  ditinguish  between  clutter  and 
target  energy.  The  ability  of  the  MTI  system  to  reject  clut¬ 
ter  and  detect  low  velocity  targets  is  in  part  based  upon  the 
clutter  cancellation  scheme  employed. 

A  MTI  radar  system  allows  for  operation  in  a  high  clutter 
environment  with  the  minimum  amount  of  dwell  time  for  a  given 
probability  of  detection.  This  is  a  necessary  feature  for  a 
spaceborne  surveillance  and  tracking  system.  A  pulse  doppler 
radar  is  a  type  of  MTI  system  that  uses  a  high  pulse  repeti¬ 
tion  frequency  to  allow  unambiguous  doppler  measurements  while 
also  making  it  possible  to  make  range  measurements.  Figure 
III-5  Illustrates  an  example  of  a  pulse  doppler  radar  MTI 

IThe  term  MTI  can  either  connote  any  radar  system  which  uses 
velocity  as  a  method  of  discriminating  between  targets  and 
clutter  or  can  connote  a  clutter  cancellation  scheme  which 
uses  delay  line  cancelers.  While  most  clutter  cancellation 
techniques  employ  some  form  of  a  delay  line  canceler,  a 
distiction  in  the  type  of  cancellation  arrangement  will  be 
made  and  the  first  MTI  definition  will  be  used. 
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system.  This  general  pulse  doppler  MTI  system  will  be  used  in 
the  description  for  this  study.  The  MTI  system  basically 
solves  the  binary  decision  problem  of  whether  a  target  is 
present  or  not.  Given  that  a  target  is  present,  the  system 
then  provides  an  estimate,  or  mark,  of  target  velocity  and 
range.  Depending  on  the  task  of  the  processor,  cost  and 
complexity  desired  and  method  of  processing,  certain  element 
blocks  may  be  added,  eliminated,  or  just  rearranged. 

The  receiver  block,  discussed  in  the  previous  section, 
performs  the  front  end  processing  on  the  signal  by  converting 
the  RF  signal  to  either  an  IF  or  video  signal.  A  possible 
modification  to  the  receiver  diagram  shown  in  Figure  III-4  is 
to  add  an  adaptive  velocity  compensation  mechanism  which 
compensates  for  the  platform  motion.  The  compensation  is 
accomplished  by  adjusting  the  frequency  of  the  COHO  as  a 
function  of  the  platform  doppler  frequency  and  sighting  vec¬ 
tor. 

If  a  modulated  waveform  is  transmitted,  a  pulse  compres¬ 
sion  filter,  in  the  form  of  a  filter  matched  to  the 
transmitted  signal,  is  used  to  compress  the  receiver  output 
signal.  In  delay  line  canceler  systems,  pulse  compression  is 
usually  done  after  clutter  processing.  Following  pulse  com¬ 
pression,  the  signal  is  range  gated  to  eliminate  excess 
receiver  noise  and  permit  target  range  measurement  (Ref  56:19- 
13).  The  number  of  range  gates  depends  on  the  number  of  range 

bins,  M,  where  M  *  2R  /cT  with  R  being  the  slant  range,  c 

s  s 

being  the  speed  of  light  and  T  being  the  compressed  pulse 


width.  If  pulse  compression  is  accomplished  digitally  then 
the  A/D  conversion  occurs  before  the  pulse  compression,  other¬ 
wise  the  signal  from  each  range  gate  is  A/D  converted. 

The  range  gated  signal  is  divided  into  quadrature  and 
in-phase  channels  to  eliminate  the  effects  of  blind  phases. 
Blind  phases  occur  when  pulse  sampling  appears  at  the  same 
point  in  the  doppler  cycle  at  each  sampling  instant  (Ref 
56:120).  The  quadrature  channel  is  especially  utilized  when 
delay  line  cancelers  are  used  as  a  clutter  canceler.  The 
digital  signal  in  the  I  and  Q  channels  are  then  digitally 
filtered  to  suppress  clutter  energy  and  enhance  the  target 
energy.  Various  clutter  cancellation  and  filtering  techniques 
will  be  discussed  later.  The  filtered  I  and  Q  channels  are 
then  combined  to  produce  a  signal  with  no  loss  (Ref  56:121). 

Following  the  combination  of  the  I  and  Q  channels,  the 
signal  from  each  filter  is  integrated,  either  coherently  or 
incoherently,  and  then  compared  in  the  CFAR  with  an 
established  threshold  level  to  determine  the  presence  or 
absence  of  a  target.  The  CFAR  circuitry,  designed  using  the 
principles  described  in  the  Detection  Processor  Development 
section  in  Chapter  II,  maintains  a  constant,  satisfactory 
number  of  false  alarms  for  the  radar  system.  The  threshold 
value  is  generally  calculated  by  applying  the  Neyman-Pearson 
criterion  using  either  fixed  or  real  time  clutter  and  noise 
statistics.  If  real  time  (updated  clutter  and  noise)  statis¬ 
tics  are  used  in  modifying  the  threshold  value,  then  the  CFAR 
system  is  considered  adaptive.  Signals  which  exceed  the 


threshold  value  are  then  passed  on  to  the  data  processor  where 
target  velocity  and  range  estimation  is  made,  and  monopulse 
calculations,  ambiguity  removal  and  track  initiation  are  all 
accomplished.  The  difference,  A  #  and  auxilary,  AUX,  radar 
channels  shown  in  Figure  III-5  may  or  may  not  be  employed. 
When  employed  they  are  used  for  monopulse  calculations,  jammer 
suppresion  and  the  adaptive  loop  calculations  such  as  the  one 
for  the  filters  or  the  CFAR  circuits.  These  channels  contain 
block  elements  whose  functions  are  similar  to  those  just 
described  for  the  sum,  2,  channel. 

In  designing  the  clutter  cancellation  techniques  care 
must  be  taken  to  maintain  perspective  on  the  relationship  of 
the  variety  of  methods  and  techniques  used  for  suppresing 
clutter.  Doppler  clutter  cancellation  is  basically 
accomplished  by  spectrally  analyzing  each  range  gated 
signal.  It  can  be  seen  that  filter  design  and  arrangement 
become  the  heart  of  any  cancellation  scheme.  There  are  two 
problem  areas  in  designing  the  filters.  The  first  is 
designing  the  optimal  filter  so  as  to  maximize  the  signal  to 
interference  (clutter  +  noise)  ratio  as  discussed  in 
References  21,  51,  58,  62.  The  second  problem  is  to 
approximate  the  optimal  filter  as  close  as  possible  digitally 
so  as  to  minimize  filter  sidelobes  which  introduce  ambiguous 
interference.  This  second  problem,  discussed  extensively  in 
Reference  29,  43,  and  60,  essentially  deals  with  choosing  the 
proper  windowing  of  the  signal,  proper  weighting  to  implement 
the  filter,  whet"  to  implement  the  filter  with  a  finite 
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impulse  response  (FIR)  circuit  or  an  infinite  impulse  response 
(HR)  circuit  or  use  a  FFT  software  package.  Windowing  and 
weighting  choices  are  influenced  by  the  clutter  spectrum  which 
must  be  cancelled,  while  choice  of  implementation  method  is 
governed  by  system  constraints. 

The  least  complex  and  most  widely  used  cancellation 
technique  is  the  delay  line  canceler.  Figure  III-6  shows  a 
single  delay  line. 


Fig.  III-6  Single  Delay  line  Canceler 

While  Figure  III-6  shows  the  delay  line  canceler  following 
pulse  compression  and  range  gating,  the  delay  canceler  may, 
and  often  does,  follow  the  phase  detector.  The  video  output  of 
the  phase  detector,  for  a  specific  range,  for  one  complete 
interpulse  period  is  fed  into  the  delayed  channel  which  has  a 
delay  time  memory  equal  to  the  transmission  interpulse  period. 
During  the  next  transmission  period,  the  output  of  the  phase 
detector  is  applied  to  the  two  channels  of  the  canceler.  At 
the  summer  junction  the  video  output  from  the  first  interpulse 
period,  which  has  been  delayed,  is  subtracted  from  the  most 
recent  undelayed  video  signal.  Theoretically,  if  there  is  no 
built-in  system  instability  or  radial  motion  caused  by  the 
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platform,  the  video  output  from  fixed  targets,  clutter,  will 

be  unchanged  between  successive  interpulse  periods  thus 

allowing  for  complete  cancellation  of  the  output  of  the 
summer.  In  contrast,  moving  targets  detected  during 
successive  interpulse  periods  will  have  a  finite  phase 

difference  causing  a  variance  in  signal  level  at  the  input  of 
the  subtracting  unit  thus  producing  a  finite  residue  output. 
The  signal  present  in  the  canceler  due  to  a  target  return  at 
time  t  is 

E  (t)  =  A  sin(  2TTE  t  +  0  +  tf)  (III-38) 

1  d 

where  0  is  the  initial  phase  and  is  the  phase  shift  due  to 
the  target  at  range  4TTR/\  .  The  signal  at  an  interpulse 
period  later  is 

E  (t)  =  A  sin[277f  (t,  +  T)  +  ©  +  «S]  (III-39) 

2  d 

.  Thus  the  output  of  the  canceler  is 

E  (t)  *  E  (t)  -  E  (t)  =  2A  sin fTTf  T) 

0  2  1  d 

cos  1 27Tf  (t  +  .  5T)  +  ©  +  tf]  ( III-40 ) 

d 

From  equation  III-40  we  can  see  that  the  magnitude  of  the 
phase  imbalance  is  thus  proportional  to  the  doppler  frequency 
shift  of  the  moving  target  By  envelope  detecting  and  doppler 
filtering  the  signal  described  by  equation  III-40,  the  target 
doppler  can  be  determined.  From  equation  III-40  we  can  also 


V*. 


see  that  system  instabilities  such  as  pulse  time  jitter,  pulse 
width  jitter,  pulse  amplitude  jitter,  pulse  distortion,  and 
frequency  modulation  within  the  pulse  all  effect  the 
cancellation  ability  of  a  delay  line  canceler. 

The  delay  line  canceler  acts  as  a  periodic  filter 
notching  out  fixed  clutter  and  signal  energy  around  the  dc  and 
the  PRF  harmonic  frequencies.  The  transfer  function  is 


I H (  f )  1 .  =  2  sin  f7Tf  T) 


( III-41  ) 


This  filter,  periodic  in  the  frequency  domain,  is  sometimes 
described  as  a  "comb  filter."  Figure  III-7  shows  the  response 
of  a  delay  line  canceler  or  comb  filter. 
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Fig.  III-7  Delay  Line  Frequency  Response 


The  motion  of  the  spaceborne  radar  platform  induces  a 
frequency  translation  of  the  clutter  spectrum.  To  compensate 
for  the  doppler  shift  in  clutter  the  frequency  in  the  COHO  can 
be  modified  to  shift  back  the  clutter  spectrum.  A  phase 
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shift,  inserted  in  one  branch  of  the  delay  line  canceler  that 

shifts  the  null  of  the  response  to  match  the  translated 
clutter  spectrum,  can  also  compensate  for  the  translational 
effects  of  the  platform  motion.  The  appropriate  phase  shift  is 
a  function  of  the  mainbeam  sighting  vector  and  platform 
velocity.  The  mean  doppler  shift  due  to  platform  motion  can 
be  removed  by  phase  locking  the  receiver  to  the  clutter  as  is 
done  with  time  average  clutter  coherent  airborne  radar 

(TACCAR)  (Ref  3:566).  As  in  the  airborne  application,  the 

! 

clutter  doppler  shift  for  the  space-based  radar  will  be  range 
dependent  since  the  doppler  frequency  is  a  function  of  the 
sighting  angle  from  the  radar  to  the  clutter  cell.  The  aver¬ 
age  doppler  frequency  of  the  clutter,  computed  by  averaging 
the  clutter  signal  phase  change  from  several  range  intervals, 
is  used  to  cause  the  mean  doppler  frequency  to  coincide  with 
the  null  of  the  delay  line  filter  response  over  the  rest  of 
the  range  of  observation  (Ref  56:142).  This  system  is  an 
example  of  a  clutter  lock  scheme.  The  clutter  lock  addition 
simply  compensates  for  the  doppler  shift  and  then  uses  any 
conventional  clutter  canceler  to  remove  the  remaining  clutter 
spectrum. 

The  clutter  locking  technique  becomes  ineffective  when 
more  than  one  type  of  clutter  is  present  in  the  range  of 
interest  because  of  the  wide  disparity  in  the  mean  velocities. 
Furthermore,  in  the  absense  of  clutter,  the  system  compensates 
for  only  the  target  velocity,  thus  rejecting  the  target.  Hence 
it  is  necessary  to  bypass  the  clutter  locking  circuits  when 


clutter  power  Is  small (Ref  40:  328). 

The  single  canceler  in  many  instances  does  not  have  a 
sufficiently  broad  enough  null  to  reject  the  extended  clutter 
spectrum  which  results  from  a  moving  platform.  To  reduce  the 
clutter  residue,  the  filter  rejection  notches  can  be  widened 
by  cascading  several  single  cancelers.  The  filter  transfer 
function  for  a  cascaded  delay  line  system  is 


n  n 

4H(f )  |  =  2  sin  (Iff T) 


. (III-42) 


where  n  is  the  number  of  cascaded  elements.  Figure  III-7  also 
shows  the  response  of  multiple  cancelers.  Care  must  be  taken 
when  widening  the  notch  that  target  energy  in  the  vicinity  of 
the  notch  is  not  also  rejected. 

The  multiple  delay  line  pulse  canceler  is  an  example  of 
a  transversal  or  comb  filter.  The  transversal  filter  finds 
much  use  in  delay  line  cancelers  because  of  its  ease  of 


digital  implementation. 


Figure  III-8  illustrates  a 


nonrecursive  or  transversal  filter, 


(Ref  56:110) 


Fig.  III-8  Transversal  Filter 


The  output  of  the  filter  can  be  written  as 


N 

E  (NT)  =  2  w  E  (kT)  ( III-43 ) 

0  k=1  k  I 


where  E  (NT)  is  the  output  of  the  filter  and  E  (kT)  is  the 
0  I 

input  to  the  filter.  Since  the  filter  response  is  controlled 

by  adjusting  the  weights  w  ...  w  ,  much  effort  has  been  spent 

1  n 

to  optimally  choose  the  weights  and  design  methods  to 
adaptively  modify  the  weights  as  a  function  of  the  clutter 
return  (Ref  2,  14,  24:91,  54).  The  transversal  filter  with 

alternating  binomial  weights 


k-1 

w  =  (-1) 
k 


k  *  1 , 2 ,  ...  n+ 1 


( III-44 ) 


closely  approximates  the  filter  which  maximizes  the  average 
improvement  factor  (Ref  56:111).  By  using  alternate  weighting 
schemes,  i.e.  Hanning,  Butterworth,  a  different  filter 
response  results  which  maximizes  performance  criteria  other 
than  the  improvement  factor. 

Since  the  transversal  filter  simply  implements  a 
cascaded  delay  line  canceler  network,  as  the  number  of  delays 
increases,  the  notches  at  dc  and  all  the  PFR  harmonics  are 
widened  and  cancel  an  increasing  amount  of  target  energy. 
Ideally  the  filter  response  should  have  notches  at  dc  and  all 
PRF  harmonics  which  match  the  clutter  spectrum.  The  filter 
response  can  be  shaped  to  approximately  match  the  clutter  by 


employing  recursive,  or  feedback,  filters.  The  z-plane  fre¬ 


quency  analysis  techniques  discussed  in  Chapter  II  are  used  to 
sysnthesize  the  desired  frequency  response.  While  recursive 
filters  offer  a  better  steady  state  response  than  nonrecursive 
filters,  they  exhibit  severe  ringing  in  the  presence  of  large 
clutter  (Ref  56:114)  and  require  a  number  of  pulses  to  be 
processed  before  the  steady-state  condition  is  met.  For  wide 
area  surveillance  application  time  is  limited.  Since  the  time 
required  for  the  transcient  response  of  a  recursive  filter  to 
die  is  large,  recursive  filters  will  not  be  considered  for  use 
in  the  delay  line  cancelers  for  this  study. 

In  a  delay  line  canceler  system,  target  doppler  is 
obtained  indirectly  from  the  received  signal.  That  is,  the 
doppler  information  describing  either  a  target  or  clutter  is 
obtained  from  the  comparison  of  several  interpulse  video  sig¬ 
nals.  The  relative  change  from  one  interpulse  period  to  the 
next  indicates  the  object's  doppler.  Low  doppler  objects  are 
then  eliminated  through  filtering.  The  next  clutter  can¬ 
cellation  technique  obtains  object  doppler  information  direct¬ 
ly  by  determining  the  frequency  shift  of  the  received  signal 
from  that  of  the  transmitted  signal  frequency. 

In  this  second  cancellation  scheme,  fixed  and  moving 
targets  are  discriminated  through  the  use  of  a  single  or  a 
bank  of  contiguous  doppler  filters.  In  this  scheme  a  typical 
clutter  cancellation  system  would  consist  of  a  bandpass  filter 
which  passes  the  expected  frequency  response  corresponding  to 
targets  and  suppressing  the  response  from  fixed  targets.  A 
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bank  of  contiguous  narrowband  f liters , which  would  extend  over 
the  expected  target  doppler  frequencies,  may  follow  the 


bandpass  filter.  The  contiguous  doppler  filters  are  more 
complex  than  simple  narrowband  filters.  These  filters  are 
designed  to  better  match  the  characteristics  of  the  clutter 
spectrum.  The  complex  narrowband  filters  permit  detection, 
measurement  of  target  velocity,  and  resolution  between  targets 
having  different  velocities  within  the  same  beamwidth  (Ref 
32:16).  Because  the  clutter  spectrum  is  predominantly  at  low 
doppler  frequencies,  the  low  doppler  filters  of  the  contiguous 
bank  of  filters  can  be  designed  differently  than  the  high 
doppler  filters.  The  low  doppler  filters  can  be  of  the  form 
suggested  by  Stutt  and  Spafford  (Refs  58,  62)  and  Delong  and 
Hofstetter  (Ref  21  ) ,  where  the  "mismatched"  filters  are  a 
function  of  the  clutter  statistics.  The  major  problem  with 
these  filters  is  their  implementation,  thus  filters  similar  to 
those  suggested  by  Urkowitz  can  be  used  (Ref  22:170). 

Figure  III-9  shows  a  system  block  diagram  and  a 
typical  response  for  a  doppler  cancellation  scheme. 
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Fig.  III-9  (a)  Block  Diagram  (b)  Typical  Response 


The  placement  of  the  bandpass  filter  within  the  receiver  is 
optional.  The  filter  can  be  placed  just  beyond  the  maximum 
anticipated  clutter  spectrum  with  the  knowledge  that  low 
radial  velocity  targets  will  be  rejected.  It  is  desirable 
when  using  a  narrowband  filtering  scheme  as  just  descibed  to 
use  a  high  PRF  waveform  to  avoid  doppler  ambiguities. 

The  cancellation  techniques  discussed  so  far  work  well 
when  there  is  a  large  differential  between  clutter  and  target 
doppler  frequencies.  A  spaceborne  radar  system's  ability  to 
distinguish  low  velocity  targets  from  clutter  targets  is 
however,  severly  impaired  because  the  high  platform,  velocity 
widely  spreads  the  clutter.  An  additional  compensation  or 
cancellation  technique  is  thus  required  to  reduce  the  spread 
of  clutter  in  a  space-based  radar  if  low  velocity  doppler 
targets  are  to  be  detected.  A  clutter  spread  compensation 
method,  mostly  employed  on  phased  array  radars,  is  the 
Displaced  Phase  Center  Antenna  (DPCA)  MTI  technique.  The  DPCA 
MTI  technique  attempts  to  remove  the  effects  of  platform 
motion,  thereby  improving  the  subclutter  visibility  of  the  low 
velocity  targets.  The  elements  which  compose  the  DPCA  MTI 


I  CSC 


system  are  shown  in  Figure  III-10. 
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Fig.  III-10  DPCA  MTI  Canceler 


The  linear  array  consists  of  two  distinct  sub-aperatures  or 
phase  centers.  In  the  DPCA  MTI  technique,  energy  from  a 

single  aperature  is  transmitted  at  a  time  t  .  All  the  return 

1 

energy  is  processed  in  the  coherent  receiver  as  in  a 
conventional  canceler.  Returns  from  the  initial  interpulse 
period  are  stored  in  the  array  delay  which  is  equal  to  the 
interpulse  period.  As  the  platform  moves  forward  during  ^t, 
where  £t  corresponds  to  the  interpulse  period,  the  phase 
center  of  the  array  is  shifted  rearward  and  a  second  transmis¬ 
sion  is  made  from  theoretically  the  same  point  in  space. 
Fixed  targets  from  the  two  successive  transmissions  will  have 
identical  phase  and  amplitudes  thus  being  cancelled  in  the 
array  delay  processor.  Moving  targets  would  have  changing 


phase  and  amplitudes.  The  difference  between  the  two  periods 
will  thus  produce  a  residue  voltage  out  of  the  processor  that 
corresponds  to  their  velocity.  Theoretically,  for  perfectly 
matching  aperature  patterns,  compatability  between  switching 
time  of  the  beam  and  the  system  PRF,  and  no  vehicle  motion 
effects,  the  DPCA  system  should  provide  high  rejection  of  the 
mainbeam  and  sidelobe  clutter,  resulting  in  good  subclutter 
visibility  and  detection  of  low  velocity  targets  (Ref  36:27). 
Matching  of  system  parameters  and  vehicle  rotation  have  how¬ 
ever  limited  the  performance  of  the  DPCA  MTI  technique.  The 
output  from  the  array  cancelling  is  then  further  processed  by 
one  of  the  cancellation  techniques  previously  discussed. 

An  improvement  to  the  DPCA  MTI  concept  is  obtained  by 
using  a  multiple  phase  center  array  combined  with  pulse 
compression  and  long  integration  time.  This  technique  is 
employed  on  the  Multiple  Antenna  Moving  Target  Surveillance 
Radar  (MASR) .  The  operation  is  identical  to  the  DPCA  MTI 
system  except  1 )  a  number  of  discrete  phase  centers  from  which 
beams  are  transmitted  allow  for  multiple  look  integration  2) 
the  output  of  the  coherent  receiver  undergoes  pulse 
compression  before  entering  the  MTI  delay  line  canceler  and 
3)  video  integration  is  performed  on  the  outputs  of  the 
canceler  (Ref  32:33).  The  system  limitations  discussed  with 
the  DPCA  MTI  system  are  now  compounded  because  of  the  multiple 
phase  centers. 

Both  the  DPCA  and  MASR  compensation  techniques  need 
multiple  "looks"  to  limit  clutter  spreading.  Multiple  looks 


reduce  the  area  search  rate  of  the  radar.  Employing  a 
technique  similar  to  that  used  in  monopulse  tracking,  clutter 
spreading  can  be  compensated  for  without  reducing  the  search 
rate.  The  monopulse  anomolous  nulling  clutter  compensation 
technique  uses  the  normal  sum  pattern  from  the  antenna  and  the 
difference  pattern  from  the  antenna  to  cancel  the  clutter  in 
the  mainbeam  without  canceling  the  target.  A  target  in  the 
mainbeam  is  near  the  peak  of  the  sum  pattern,  thus  having  a 
large  amplitude  while  having  a  low  amplitude  in  the  difference 
pattern.  The  clutter  in  the  mainbeam  is  doppler  spread  by  the 
platform  motion  as  described  before.  The  target  doppler  shift 
competes  with  the  clutter  doppler  shift  from  a  different  angle 
within  the  mainbeam.  Essentially  the  clutter  is  off  the  peak 
of  the  beam.  The  clutter  on  the  other  hand  has  a  large  gain 
in  the  difference  pattern.  Modifying  the  gain  in  the 
difference  pattern  to  have  the  same  amplitude  as  the  sum 
pattern,  the  difference  clutter  pattern  is  subtracted  from  the 
sum  clutter  pattern,  thus  canceling  out  the  clutter.  Since 
the  target  amplitude  in  the  difference  pattern  is  nearly  zero, 
while  the  target  amplitude  in  the  sum  pattern  is  large,  the 
resulting  target  amplitude  is  not  sufficiently  reduced  when 
differencing  occurs.  Following  the  differencing,  the  signal 
is  then  processed  further  by  one  of  the  previously  described 
cancellation  techniques  to  eliminate  the  remaining  clutter 
energy  and  provide  for  range  and  velocity  measurements.  While 
this  method  enhances  the  S/C  ratio  of  the  target  in  the 
mainbeam,  sidelobe  clutter  is  limited.  This  method,  like  the 


DPCA  and  MASR  techniques,  is  limited  by  the  total  system 
stability. 

Various  clutter  cancellation  techniques  have  been 
presented  in  this  chapter  which  take  advantage  of  the 

different  spectral  characteristics  between  targets  and 
clutter.  While  this  is  not  an  exhaustive  list  of  all  the 

techniques  available,  the  cancellation  schemes  presented  are 
representative  of  the  manner  in  which  clutter  for  a  spaceborne 
radar  system  would  be  cancelled.  Essentially  each  of  these 
techniques  distinguish  between  targets  and  clutter  doppler 
characteristics  by  using  some  form  of  a  delay  line  canceler  or 
a  bank  of  contiguous  doppler  filters.  To  enhance  the  radar's 
ability  to  detect  low  doppler  targets  DPCA,  MASR  or  monopulse 
anomolous  nulling  can  be  added  to  compensate  for  the  clutter 
spreading  due  to  the  high  platform  velocity.  Because  the 

clutter  spreading  compensation  techniques  involve  total  radar 
system  aspects,  they  will  not  be  added  in  the  modeling  and 
analysis  of  the  clutter  cancellation  techniques  found  in  the 
next  chapter. 


IV  Computer  Models 


Approach 

The  modeling  of  the  clutter  cancellation  problem  can  be 
divided  into  three  separate  areas  or  subprograms.  The  calcu¬ 
lation  of  the  clutter  power  spectral  density,  for  given  radar 
parameters  and  given  surveillance  scenerio  geometry  is  the 
first  area.  The  routines  used  to  calculate  clutter  PSD  can  be 
exercised  to  obtain  insights  and  understanding  into  the  magni¬ 
tude  and  extent  of  the  clutter  PSD  for  various  surveillance 
situations  and  clutter  statistics.  The  second  computer 
modeling  area  is  the  processing  of  the  radar  echo  signal  so  as 
to  cancel  clutter  and  thus  obtain  target  information.  Since 
the  analysis  of  Chapters  II  and  III  showed  that  the  modeling 
of  a  clutter  cancellation  processor  essentially  becomes  a 
matter  of  implementing  a  particular  digital  filtering  scheme, 
the  second  subprogram  simulates  a  digital  filter.  The  final 
computer  modeling  area  deals  with  the  calculations  to  obtain 
the  attenuation  factor  for  the  various  clutter  models  -  fil¬ 
tering  scheme  combinations.  The  output  from  this  subprogram 
provides  a  quantitative  measure  of  the  clutter  cancellation 
ability  of  various  filtering  schemes  against  different  clutter 
backgrounds. 

The  computer  simulation  written  for  this  study  is 
structured  such  that  variations  in  the  clutter  background, 
or  cancellation  scheme,  can  be  easily  implemented  by  simple 
substitutions  of  routines.  The  main  program,  CLUTCAN,  calls 


the  subprograms  that  model  each  of  the  areas  described  above. 
These  subprograms  in  turn  call  the  subroutines  that  comprise 
each  subprogram.  The  CLUTCAN  program  simulates  the  entire 
clutter  cancellation  problem,  from  the  calculation  of  clutter 
PSD,  the  filtering  schemes  that  cancel  selected  clutter  spec¬ 
tral  components,  to  the  calculation  of  the  attenuation  factor 
of  the  cancellation  scheme  used.  Figure  IV-1 ,  shown  on  the 
next  page,  shows  the  flowchart  for  the  CLUTCAN  program.  Al¬ 
though  the  subprograms  are  exercised  by  the  CLUTCAN  program, 
they  each  can  easily  be  used  in  separate  programs  to  analyze 
specific  problems.  For  example,  the  effects  of  various  clut¬ 
ter  statistics  or  variations  in  surveillance  geometries  on 
clutter  PSD  can  be  obtained  by  writing  a  seperate  program  that 
only  calls  the  clutter  PSD  calculation  subprogram. 

Since  the  main  program,  CLUTCAN,  is  only  an  initializa¬ 
tion  and  calling  program,  and  will  not  be  dicussed  in  detail. 
The  CLUTCAN  program  listing  can  be  found  in  Appendix  D.  The 
next  three  sections  will  discuss  the  three  subprograms  called 
by  CLUTCAN.  The  subprograms'  flowcharts  and  program  listings 
can  also  be  found  in  Appendix  D. 

Clutter  PSD  Calculations 

The  subprogram  CLUT  calculates  the  power  spectral 
density  that  is  received  by  a  space-based  doppler  radar  system 
given  radar  system  parameters  and  the  viewing  geometry 
parameters  for  a  specific  surveillance  scenerio.  This 
subprogram  computes  the  clutter  PSD  by  one  of  two  methods. 


INITIALIZE  VARIABLES 
DECLARE  CONSTANTS 


The  first  method  simply  assumes  a  Gaussian  power  spectral 
distribution  function,  as  described  by  equation  III-27.  The 
PSD  mean  power  value  is  then  computed  by  equation  IV-1 


W  =  P  G2  X2  o' 
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where  the  spectrum  variance,  cr4  ,  is  computed  via  equation 

c 

III-35.  The  clutter  spectrum  is  defined  to  be  centered  at 
zero  frequency  because  it  is  assumed  that  platform  induced 
frequency  translation  is  compensated  for  by  adjustments  in  the 


COHO. 


The  second  clutter  PSD  calculation  method  is  based  upon 


the  summation  of  power  reflected  from  range  and  range  rate 
distributed  clutter.  The  region  of  distributed  clutter  that 
is  of  interest  is  defined  by  the  area  of  the  earth  illuminated 
by  the  radar  for  a  pair  of  depression  and  azimuth  angles. 
These  angles  are  measured  relative  to  the  velocity  vector  of 
the  spaceborne  platform  as  illustrated  in  Figure  IV-2.  The 
clutter  level  for  an  arbritrarily  selected  range-doppler  cell 
within  the  illuminated  area,  where  the  clutter  cell  range 
width  is  defined  by  the  range  resolution  constant,  equation 
11-123,  and  the  clutter  cell  doppler  width  is  defined  by  the 
doppler  resolution  constant,  equation  11-126,  is  calculated  by 
multiplying  the  average  backscatter  coefficient  with  the  cell 
area. 


,  -  . 


„  v  • 

\  iW, 


§S  •  % 

s!!. 


Fig  IV- 2  Situation  Geometry 


For  the  high  PRF  case,  the  clutter  power  from  the  ambiguous 
range  cells  are  summed  for  an  unambiguous  doppler  cell  to 
obtain  the  clutter  power  for  a  specific  doppler. 

Projected  beneath  the  radar  is  a  locus  of  points  which 
form  a  sphere  centered  at  the  radar  nadir.  This  locus  of 
points  are  instantaneously  at  a  fixed  range,  R,  from  the 
radar.  The  intersection  of  the  spheres  with  the  earth  forms  a 
circle  called  an  isorange,  all  points  of  which  contribute  to 
the  range  gate  containing  the  range  R.  The  earth  isorange 
circles  are  also  illustrated  in  Figure  IV-2.  A  range  gated 
return  will  contain  clutter  energy  from  the  region  between  two 
isorange  circles  (Ref  50:375). 

In  the  frequency  domain,  the  locus  of  fixed  points  in 


I  • 


space,  from  which  a  return  would  have  a  given  doppler  shift, 

f  ,  is  a  half  cone  whose  vertex  is  at  the  radar  and  whose  axis 
o 

is  the  velocity  vector  of  the  radar.  The  half  cone  arises 
because  the  instantaneous  relative  speed  between  the  radar  and  a 
fixed  point,  is  the  radar  speed  multiplied  by  the  cosine  of 
the  angle  formed  by  the  velocity  vector  and  the  position 
vector  of  the  point  relative  to  the  radar.  The  half  cone 
formed  by  the  set  of  points  with  fixed  doppler  is  illustrated 
in  Figure  IV-3. 


-2 


(Ref  47:XVI-7 ) 

Fig.  IV-3  Generation  of  Constant  Doppler  Contours 

The  intersection  of  the  half  cone  with  the  earth  forms  a  conic 

section  called  an  isodop  all  points  of  which  contribute  to  the 

doppler  cell  containing  the  doppler  f  .  Figure  IV-4  shows  the 

o 

projection  of  the  isoranges  and  isodops  on  a  flat  surface. 
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(Ref  47:XVI-8) 

Fig.  IV-4  Isodop  and  Isorange  Points 

The  doppler  cell  thus  contains  the  clutter  returns  from  the 

region  between  two  isodop  conic  sections  (Ref  50:375). 

To  compute  the  PSD,  the  illuminated  range  cell  region, 

shown  in  Figure  IV-4  between  the  two  isodops Af  a  part,  is 

d 

first  determined.  This  area  is  the  clutter  area  for  an 

arbitrarily  selected  range-doppler  cell.  This  clutter,  area, 

.1 

multiplied  by  oJ  is  used  in  the  radar  range  equation  (along 
with  the  appropriate  two  way  antenna  gain  defined  along  the 
observation  line  of  sight)  to  calculate  the  region's 
contribution  to  the  platform  motion  power  spectrum.  Figure 
IV-5  illustrates  the  region  which  generates  the  clutter  for  a 
range-doppler  cell. 


0  0 

Since  o r*  varies  with  problem  geometry,  the  c r*  value  which 

corresponds  to  the  problem  grazing  angle  must  be  determined 
and  used. 


Fig.  IV-5  Range-Doppler  Clutter  Cell 


The  complete  power  spectra  can  be  computed  by  integrating 
over  all  isodop  regions  of  interest.  The  high  PRF  used  for  a 
space-based  radar  causes  ambiguities  in  range  which  necessi¬ 
tates  integration  over  the  entire  ambiguous  range  region  for 
each  unambiguous  doppler.  The  clutter  PSD  calculation  for  a 

doppler  cell,  f  ,  is  summarized  by  equation  IV-2 
o 
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where  R  is  the  first  unambiguous  range  between  0  and  c/2PRF 


i 

on  which  the  range  gate  is  centered  and  R  forms  the  set  of 


discrete  ranges 


ik 


R  =  R  +  c/2PRF 


( IV- 3 ) 


The  range  zone  width,  AR»  is 


AR  =  c/2PRF  ( IV-4 ) 

The  approach  just  described  to  calculate  the  clutter  PSD  is 
based  upon  the  clutter  modeling  developed  in  Refs  4,  25-27, 
and  50.  Chapter  2  of  Ref  27  provides  a  complete  mathematical 
description  of  this  clutter  PSD  calculation  method,  therefore 
it  will  not  be  repeated  in  this  text.  The  mathematical  equa¬ 
tions  are,  however,  applied  in  the  computer  program. 

Because  of  the  iterations  needed  to  calculate  the 
clutter  PSD  and  the  ease  with  which  variations  in  the  problem 
can  be  accomplished,  subroutines  were  used  extensively  in 
CLUT,  especially  in  the  second  method  of  clutter  computations. 
Method  two  of  CLUT  calls  several  geometry  calculation  routines 
to  compute  the  angle  and  range  relationships  between  the  radar 
coordinate  system  and  the  coordinate  system  of  the  clutter 
scattering  region.  The  subroutine  SIGMAO  computes  the  mean 
backscatter  coefficient  for  different  grazing  angles.  The 
subroutine  ANTG  computes  the  antenna  gain  as  seen  by  the 
doppler-range  cell.  A  short  description  of  each  subroutine's 
usage  is  provided  with  the  code  listing  found  in  Appendix  D. 
The  development  of  the  equations  used  in  the  subroutines  are 
straightforward  geometric  relationships,  therefore  will  not  be 
covered  here.  The  calling  sequence  of  these  routines  is  shown 
in  CLUT's  flowchart,  also  found  in  Appendix  D. 
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Clutter  Cancel er  Simulator 

The  cancellation,  or  in  some  cases  the  reduction  of 
clutter  PSD,  is  accomplished  by  digital  filtering  the  PSD 
calculated  by  subprogram  CLUT.  The  subprogram  CANCEL  simu¬ 
lates  the  clutter  cancellation  via  application  of  the  digital 
filtering  techniques  discussed  in  the  "Digital  Signal 
Processing"  section  of  Chapter  II,  primarily  the  application 
of  equation  11-51 .  Equation  11-51  is 

j  277fT  2 

S  (f)  =  IH(e  )i  S  (f)  ( IV-5 ) 

y  c 

j2fl7£T 

where  H(e  )  is  the  transform  of  the  filter  h(nT).  The 

subroutine,  FFT2C,  from  the  Cyber  IMSL  library  is  used  to 

compute  the  filter  transform  of  the  real  valued  sequence 

h(nT).  The  subroutine  requires  as  inputs  the  h(nT)  vector  of 

sequence  values  and  N,  the  number  of  data  points  to  be 

transformed,  where  as  mentioned  before,  N  should  be  a  power  of 

two.  The  subroutine  outputs  the  filter  transform  as  a  complex 

vector  of  length  N.  The  magnitude  of  the  components  of  H(k), 

taken  at  each  digital  frequency,  k,  is  obtained  by  taking  the 

squared  magnitude  of  the  subroutine  returned  vector.  The 

filtered  clutter  PSD  is  obtained  by  application  of  equation 

IV-5.  The  clutter  power  spectra,  S  (f),  is  converted  from  the 

c 

analog  frequency  axis,  f,  to  the  digital  frequency  axis,  k,  by 
k=fNT.  The  filtered  clutter  PSD,  S  (k),  is  converted  back  to 

y 

the  analog  frequency  axis  by  equation  11-103. 


The  subprogram  CANCEL  is  designed  such  that  either  a 
nonrecursive  filter  with  one  of  the  following  weighting 
schemes  -  binomial,  Hanning,  Hamming,  Blackman,  uniform, 
Butterworth  or  a  filter  impulse  function,  h(t),  is  utilized  to 
perform  the  filtering  that  accomplishes  the  clutter 
cancellation.  The  vector  of  filter  sequence  values  is  ob¬ 
tained  by  sampling  either  the  weighting  functions  or  h(t) 
every  T  seconds  such  that  N  samples  are  taken  with 
ns=0, 1  ,2, . .  .N-1 . 

The  subprogram  CANCEL  returns  to  the  main  program, 
CLUTCAN,  the  filtered  clutter  PSD  and  the  coefficients  of  the 
transformed  filter.  The  flowchart  and  program  listing  for 
CANCEL  can  be  found  in  Appendix  D. 

Attenuation  Factor  Calculation 

The  clutter  attenuation  factor,  defined  in  Chapter  I  as 

CA  =  C  /C  (IV- 6) 

in  out 

provides  the  quantitative  measure  for  evaluating  the  merits  of 
various  cancellation  schemes.  The  subprogram,  CLTATT , 
computes  the  improvement  factor  by  using  the  results  of  the 
subprograms  CLUT  and  CANCEL. 

The  clutter  power  at  the  filter  input  is  obtained  by 
simply  summing  the  clutter  PSD  values  from  the  subprogram  CLUT 
and  likewise  the  clutter  power  at  the  output  of  the  filter  is 
obtained  by  summing  the  filtered  clutter  PSD  values  from 
subprogram  CANCEL.  The  clutter  attenuation  factor  is  thus 


determined  by  substitution  of  the  clutter  power,  both  at  the 
input  and  output  of  the  filter,  into  equation  IV-7.  The 
flowchart  and  program  listing  for  CLTATT  can  also  be  found  in 
Appendix  D. 

By  utilizing  the  clutter  cancellation  problem  simulator, 
CLUTCAN,  the  clutter  rejection  capability  of  various 
cancellation  schemes  against  different  clutter  environments, 
different  target  parameters,  and/or  with  different  radar 
parameters  can  be  tested.  In  the  next  chapter,  various  test 
cases  and  the  associated  clutter  cancellation  results  will  be 
presented  and  analyzed. 


V  Results  and  Analysis 


Test  Case 


Using  the  computer  simulation  of  the  clutter  cancella¬ 
tion  problem  outlined  in  the  previous  chapter,  there  are  two 
areas  in  which  analysis  was  performed.  The  first  analysis 
area  deals  with  the  extent  and  magnitude  of  clutter  PSD  as 


seen  by  a  space-based  radar. 


Clutter  PSD  curves  were 


generated  for  different  depression/azimuth  angle  combinations 

using  both  clutter  PSD  calculation  methods.  Due  to  system 

constraints,  (Ref  4:4-29)  the  grazing  angle  was  limited 
o  o 

between  2  and  70  ,  thereby  limiting  corresponding  depression 

o  o 

angles  to  between  35  and  75  .  Since  the  calculation  of 

clutter  returns  is  symmetric  in  each  of  the  four  quadrants  of 

the  isodop-isorange  plane  (except  for  doppler  sign  changes) 

o  o 

the  azimuth  angle  is  restricted  between  0  and  90  .  Clutter 

PSD  calculations  are  computed  with  the  depression  angles 

o  o 

having  an  initial  value  of  35  and  incremented  1 0  and  the 

o 

azimuth  angle  having  an  initial  value  of  0  and  incremented 
o  o 

30  with  the  exception  of  a  data  point  taken  at  45  .  Thus, 

there  are  twenty-five  depression  angle/azimuth  angle  combina¬ 
tions  used  in  calculating  clutter  PSD  curves. 

t 

The  clutter  PSD  curves  generated  assumed  a  land 

background.  For  comparison  purposes,  the  clutter  returns  for 

sea  and  snow/ ice  were  also  generated.  Since  the  backscatter 

coefficient,  or**,  is  not  azimuth  dependent,  but  rather 

o 

depression  angle  dependent,  an  azimuth  angle  of  45  was  used 
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O  O  O  O 

with  depression  angles  of  35  ,  55  ,  and  75  .  In  all  cases  c r* 

was  assumed  constant. 

The  second  area  of  analysis  deals  with  the  cancellation 

ability  of  the  different  filter  systems.  This  study  used  a 

nonrecursive  filter  with  a  varying  number  of  taps.  Seven 

different  weighting  schemes  were  used:  Bartlett,  Hanning, 

Hamming,  Blackman,  Binomial,  Uniform  and  Butterworth. 

Although  these  weighting  schemes  may  not  represent  optimal 

weighting  methods,  they  do  illustrate  the  cancellation  ability 

of  easily  generated  weights.  Through  a  subroutine  change, 

more  complicated,  better  "matched"  filters,  similar  to  those 

suggested  by  Delong  and  Hofstetter  (Ref  21 )  and  Spafford  (Refs 

58  and  62)  could  be  implemented  and  tested.  The  cancellation 

ability  of  the  seven  filter  schemes  were  tested  against  ground 

clutter,  generated  using  the  Gaussian  function,  with  azimuth 
o  o  o  o  o 

angle  of  0  ,  45  ,  and  90  and  depression  angles  of  35  ,  55  , 

o 

and  75  . 

Since  this  study  was  not  intented  to  be  a  radar  system 

design  study,  the  radar  parameters  described  in  Appendix  A 

were  used  for  all  test  cases.  A  nonfluctuating  target,  with  a 

13  dbsm  average  RCS,  at  a  range  of  2350  kilometers  was  used 

for  analysis.  This  range  corresponds  to  a  target  at  an 

o 

altitude  of  4.5  kilometers  with  a  45  depression  angle 
relative  to  the  radar  system.  The  target  was  assumed  to  have 
a  velocity  of  617  m/s. 


Clutter  Comparisons 


In  this  section  the  clutter  spectrum,  seen  by  a  space 
based  radar  system  is  illustrated.  Figures  V-1  through  V-5 
show  the  clutter  PSD  for  various  depression  angle/azimuth 
angle  combinations  generated  using  the  Gaussian  function 
approximation  method.  The  various  curves  in  each  figure 
represent  different  azimuth  boresight  angles.  The  clutter 
returns  from  sea  and  snow/ice  backgrounds  are  illustrated  in 
Figures  V-6  and  V-7.  The  different  curves  in  each  figure 
represent  different  boresight  depression  angles.  The  sea  and 
snow/ice  clutter  returns  assumed  a  Gaussian  function. 


B6Ti0  ik  m  iOS  3 '0.00  ,40.00  sfc.  m 

FREQUENCY  CHZ)  *10l 


Fig.  V-1  Ground  Clutter  PSD  for  35°Boresight  Depression  Angle 
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10.  00  20.  00  30.  00  ,  40.  00  50.  0 

FREQUENCY(HZ)  M01 _ 

Fig.  V-5  Ground  Clutter  PSD  for  75°  Eoresight  Depression  Angle 
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From  Figures  V-1  through  V-5,  we  observe  that  by 

increasing  the  antenna  boresight  azimuth  angle  the  clutter 

spectrum  spread  increases.  Using  Figure  V-3  for  example,  the 

o 

3db  power  point  occurs  at  6.25  Hz  for  0  azimuth  angle,  at 
o  o 

73.75  Hz  for  30  azimuth  angle,  at  100.98  Hz  for  45  azimuth 

o 

angle,  at  124.96  Hz  for  60  azimuth  angle  and  at  142.5  Hz  for 
o 

90  azimuth  angle.  This  is  as  expected  from  equation  III-35, 

which  explained  the  causes  of  clutter  spread. 

Increasing  the  depression  angle,  which  likewise 

increases  the  grazing  angle,  effects  the  clutter  spectrum  in 

two  ways.  First,  the  greater  the  depression  angle  the  less 

the  clutter  spread,  again  per  equation  III-35.  This  is 

o 

illustrated  in  Figures  V-1  through  V-5.  Using  the  45  azimuth 
angle  as  a  comparison  point,  Table  V-1  illustrates  the 
decrease  in  frequency  where  the  3  db  power  point  occurs  for 
increasing  depression  angles. 

Table  V-1 

Depression  Angle  Influence  on  Clutter  Spread 


Depression 

Frequency  of 

Angle 

3  db  Point 

o 

35 

148  Hz 

o 

45 

127.5  Hz 

o 

55 

105  Hz 

o 

65 

66.25  HZ 

o 

75 

47.50  Hz 

Prom  Table  V-1  we  see  a  4.92  db  reduction  in  frequency 
°  o 

spread  for  a  75  depression  angle  relative  to  a  35  depression 

angle.  Because  MTI  target  detection  is  based  upon  doppler 

differential,  it  is  important  to  understand  the  amount  of 

clutter  spread  associated  with  different  viewing  geometries. 

Surveillance  strategies  should  try  to  operate  with  look  angles 

so  as  to  minimize  clutter  spread. 

The  second  effect  depression  angle  has  on  clutter 

returns  is  an  increase  in  clutter  power  with  increasing 

depression  angle.  Table  V-2  shows  the  increase  in  received 

clutter  power  for  different  depression  angles  over  the  clutter 

o 

power  received  for  a  35  depression  angle. 


Table  V-2 

Depression  Angle  Influence  on  Received  Clutter  Power 


Depression 

Anqle 

Clutter  Power 
Increase 

o  i 

45 

1 .87  db 

o 

55 

4.61  db 

o 

65 

9.25  db 

o 

75 

_ _ _ ! 

15.76  db 

The  increase  in  power  is  because  the  range  to  the  clutter 

patch  is  reduced  with  increasing  depression  angle. 

Figure  V-6  and  V-7  show  the  clutter  spectrum  for  sea  and 

snow/ice  backgrounds  respectively.  Comparing  the  curves  of  V- 

o 

6  and  V-7  with  the  45  azimuth  angle  curves  of  Figures  V-1 


m 
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% 
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through  V-5,  we  see  very  little  change  in  clutter  spread  for 
different  backgrounds.  Although  the  clutter  spread  from  sea 
and  snow/ice  backgrounds  differs  little  from  ground  back¬ 
grounds,  there  is  however,  a  reduction  in  amount  of  clutter 
power.  Tables  V-3  and  V-4  indicate  the  reduction  in  clutter 
return  from  sea  and  snow/ice  as  compared  to  ground  clutter 
returns . 


Table  V-3 

Comparison  of  Ground  /  Sea  Clutter  Power 


Depression 

Anqle 

Ground /Sea 
Power  Ratio 

o 

35 

-23.26  db 

o 

55 

-12.83  db 

o 

75 

-7.64  db 

Table  V-4 

Comparison  of  Ground  /  Snow-Ice  Clutter  Power 


Depression 

Angle 

Ground/ Snow-Ice 
Power  Ratio 

o 

35 

-8.5  db 

o 

55 

-3.86  db 

o 

75 

-4.0  db 

The  reduction  in  sea  and  snow/ice  clutter  power  return  is  due 
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to  the  increase  in  diffuse  reflection  that  occurs  in  sea  and 
snow/ice  backgrounds. 

Due  to  a  limitation  in  time,  the  second  method  for 
clutter  PSD  calculation  was  not  executed  and  compared  with  the 
Gaussian  clutter  PSD  calculation  method.  Had  this  second 

method  been  run,  we  would  have  seen  a  greater  spread  in 
clutter  as  compared  to  the  first  method  and  an  increase  in  the 
amount  of  clutter.  This  is  because  the  second  method 
calculates  the  clutter  contribution  of  all  illuminated  earth 
area  and  not  just  the  area  illuminated  by  the  mainbeam. 


Cancellation  Scheme  Comparison 

The  ability  of  a  particular  cancellation  system  to 
eliminate  clutter  is,  as  seen  from  the  development  in  Chapter 
III,  based  upon  the  ability  of  a  filter,  or  set  of  filters,  to 
cancel  the  spectral  components  which  correspond  to  clutter, 
while  preserving  target  spectral  energy.  Regardless  of 
whether  the  spectral  information  is  obtained  indirectly,  as  in 
a  delay  line  canceler  system,  or  directly,  an  indication  of 
the  cancellation  ability  of  any  doppler  clutter  cancellation 
method  can  be  obtained  by  analyzing  the  filters. 

Because  the  transversal  filter  can  be  used  to  implement 
a  delay  line  canceler  system  or  a  filter  bank  system,  it  is 
used  in  this  study  as  a  basis  for  comparison  of  different 
clutter  cancellation  systems.  That  is,  it  is  assumed  that 


different  configurations  of  the  transversal  filter  represent 


different  cancellation  systems.  The  filters  analyzed  are 
intended  to  represent  front  end  filters,  cancelling  the  clut¬ 
ter  spectrum,  while  allowing  for  signal  detection  in  seperate 
filter  circuitry.  To  maintain  perspective  as  to  where  these 
filters  fit  into  the  radar  signal  processing  chain,  the 
transversal  filters  analyzed  would  be  used  at  the  doppler 
filter  bank  block  in  Figure  III-5. 

Before  analyzing  the  clutter  cancellation  ability  of  a 
multitude  of  transversal  filters  against  the  various  clutter 
spectrum  described  in  Figures  V-1  through  V-5,  a  pragmatic 
approach  to  limit  the  number  of  filters  studied  was  necessary. 
The  ideal  filter  response  for  cancelling  clutter,  as  described 
in  Chapter  III,  was  one  that  possesed  a  notch  at  dc  and  all 
PRF  harmonics  wide  enough  to  cancel  clutter  yet  not  too  wide 
so  as  to  eliminate  targets.  There  are  three  parameters  in  a 
transversal  filter  which  can  be  varied  to  produce  different 
spectral  responses  -  1)  the  weighting  on  each  tap,  2)  the  time 
delay  between  taps,  and  3)  the  number  of  taps.  Figures  V-8 
through  V-1 4  show  the  filter  magnitude  response  due  to  Butter- 
worth,  Bartlett,  Hanning,  Hamming,  Blackman,  Binomial  and 
uniform  weighting  schemes.  Figures  V-1 5  through  V-1 8  show  the 
filter  response  of  a  different  dela^  'ne  time.  Finally 
Figures  V-1 9  through  V-22  show  the  effect  of  varying  the 
number  of  filter  taps.  The  filters  illustrated  in  all  the 
figures,  Figures  V-8  through  V-22  all  have  a  periodic  spectral 
response,  with  period  1/T.  T  is  the  time  delay  between  taps. 

The  number  next  to  the  curves  in  each  figure  indicates 
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Fig.  V-9  Transversal  Filter  with  Hanning  Weights 
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Fig.  V-JO  Transversal  Filter  with  Hanning  Weights 
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Fig.  V-l  1  Transversal  Filter  with  Hamming  Weights 


Fig.  V-12  Transversal  Filter  with  Blackman  Weights 
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V-13  Transversal  Filter  with  Binomial  Weights 
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19  Butter»orthr.Weighted  Filter  (Even  number  of  Taps) 
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•21  Butterworth  Weighted  Filter  (Odd  number  of  Taps) 


HANNING  WEIGHTS 

1  /fYs.  (N  =  128) 


3>5‘ 


* 

A  VV 
(9  V  f%\ 
§  -.A 


00  60,  00  120.  00  180.  00,  240. 00  3 

FREQUENCY  (HZ)  *10* 


•AT  7 


,  % 

*.  % 


BINOMIAL  WEIGHTS 
(N  =  128) 


2* 


.*  f+-  i 


4^  ,  , 

'uSS  6k  00  1^3. 00  lbs.  00,  240. 00  3&0. 00 

FREQUENCY  <HZ>  *10* 


Fig.  V-23  Binomial  Weighted  Filter  with  Different  Number  of  Taps 


Comparison  of  Figures  V-8  through  V-14  shows  a  great 
similarity  between  Bartlett,  Hanning,  Hamming,  Blackman  and 
uniform  weighting  schemes.  This  is  to  be  expected  since  the 
weight  values  for  these  five  schemes  are  similar,  varying  from 
between  +_  1  .  The  Butterworth  and  Binomial  responses  are 
characterized  by  a  broad  mainlobe  response,  broader  than  the 
mainlobe  response  of  the  other  five  weighting  schemes,  and 
essentially  nonexistent  sidelobes.  (The  Butterworth  filter 
has  a  3  db  bandwidth  of  7125  Hz.  the  Binomial  filter  has  a  3 


db  bandwidth  of  8325  Hz  and  the  Hanning  filter  has  a  3  db 
bandwidth  of  only  4500  Hz.)  The  low  sidelobes  and  broader 
mainlobes  are  due  to  the  dynamic  range  of  the  Butterworth  and 
Binomial  weight  values.  The  slight  sidelobes  which  the 
Butterworth  filter  exhibits  is  due  to  the  asymmetric  nature  of 
the  Butterworth  weights. 

Figures  V-15  through  V-18  show  the  effect  of  changing 
the  time  delay  between  taps.  Increasing  the  delay  time  has 
three  effects  on  the  filter  response  :  1 )  reduces  the  band¬ 
width  of  the  filter,  2)  reduces  the  depth  of  the  sidelobe 
nulls,  and  3)  reduces  the  sidelobe  levels.  Because  the  filter 
response  with  delay  times  greater  than  33.333E-6  (reciprocal 
of  the  possible  target  doppler  bandwidth)  do  not  span  the 
entire  spectrum  in  which  a  target  could  be  present,  and  delay 
times  less  than  33.333E-6  span  a  spectral  region  greater  than 
necessary,  thereby  increasing  the  amount  of  noise,  only  delay 
time  equal  to  33.333E-6  seconds  needs  be  used.  If  a  bank  of 
narrowband  filters  is  to  be  implemented,  then  delay  times 
greater  than  33.333E-6  would  be  necessary,  with  the  delay 
dependent  on  the  bandwidth  of  the  filter  desired.  The  trends 
that  are  seen  with  a  filter  having  a  delay  time  equal  to 
33.333E-6  will,  however,  apply  to  a  filter  with  different 
delay  times. 

Figures  V-19  through  V-22  illustrate  Binomial, 
Butterworth  and  Hanning  weighting  schemes  with  different 
number  of  tap  weights.  The  Bartlett,  Hamming,  Blackman,  and 
uniform  weighting  schemes  were  eliminated  because  of  the 


similarity  with  the  Hanning  weighting  scheme.  The  mainlobe  of 
the  filter  response  using  either  the  Butterworth,  Hanning,  or 
Binomial  weighting  schemes,  broadens  as  the  number  of  taps  is 
reduced.  This  is  again  due  to  the  reduction  in  dynamic  range 
of  the  weights.  The  Butterworth  filter  and  the  Hanning 
filter  with  an  odd  number  of  taps,  both  exhibit  nonzero  filter 
responses  at  frequencies  less  than  400  Hz  where  clutter  is  the 
greatest.  This  undesired  high  response  worsens  as  the  number 
of  taps  decreases.  Because  of  the  high  filter  response  in  the 
low  doppler  regions,  the  Butterworth  filters  and  the  Hanning 
filter  with  an  odd  number  of  taps,  can  be  eliminated  from 
further  analysis. 

After  systematic  analysis  of  various  filters,  only 

filters  with  Binomial  weights  and  Hanning  weights  with  an  even 

number  of  taps,  need  to  be  analyzed  for  clutter  cancellation 

ability.  It  is  interesting  to  note  that  a  Hanning  filter  with 

two  taps  has  the  same  response  as  a  Binomial  with  two  taps, 

therefore  the  two  tap  Hanning  filter  can  also  be  eliminated 

from  further  analysis.  The  curves  in  Figures  V-1  through  V-5 

show  that  the  most  stressing  clutter  case  is  when  the  antenna 

o  o 

boresight  has  a  35  depression  angle  and  a  90  azimuth  angle. 
Table  V-5  illustrate  the  clutter  attenuation  factor  of  a 
transversal  filter  with  Binomial  weights  with  different  number 
of  taps  and  a  filter  with  Hanning  weights  with  four  and  six 


Table  V-5 


Clutter  Attenuation  Factor  for  Different  Cancellation  Schemes 


Weighting 

Scheme 

Number  of 

Taps 

Clutter  Attenuation 
Factor 

Binomial 

2 

30.697  db 

Binomial 

4 

76.623  db 

Binomial 

5 

96.904  db 

Binomial 

7 

134.681  db 

Binomial 

8 

152.363  db 

Hanning 

4 

24.694  db 

Hanning 

6 

21.201  db 

From  Table  V-5  we  see  there  is  a  122  db  improvement  in 
clutter  attentuation  by  using  an  eight  tap  Binomial  filter 
instead  of  a  two  tap  Binomial  filter,  however,  this  is  at  the 
expense  of  detecting  targets.  Comparing  the  Hanning  and 
Binomial  weighting  schemes  we  observe  a  32  db  reduction  in 
attenuation  using  a  four  tap  Hanning  filter  instead  of  a  four 
tap  Binomial  filter  and  a  94  db  reduction  in  attenuation  using 
a  six  tap  Hanning  instead  of  a  six  tap  Binomial. 

Table  V-5  provides  attenuation  figures  for  only  one 
particular  viewing  geometry.  With  different  viewing 
geometries  the  clutter  attenuation  will  vary  because  the 
clutter  spread  varies,  thereby  changing  the  ratio  of  the 
clutter  into  and  out  of  the  filter.  Table  V-6  shows  the 
clutter  attenuation  factor  for  different  viewing  geometries 
using  a  four  tap  Binomial  filter. 


Table  V-6 


Clutter  Attenuation  Factor  for  Different  Viewing  Geometries 


Azimuth 

Angle 

. 

Depression 

Angle 

Clutter  Attenuation 
Factor 

o 

45 

o 

35 

86.144  db 

o 

o 

90 

35 

76.623  db 

o 

o 

45 

55 

94.857  db 

o 

o 

90 

55 

86.374  db 

o 

o 

45 

75 

157.264  db 

o 

o 

90 

75 

118.190  db 

The  average  clutter  attenuation  factor  provides  a  better  indi¬ 
cation  of  the  cancellation  ability  for  a  surveillance  radar 
since  the  viewing  geometries  continually  change.  The  mean 
clutter  attenuation  for  a  four  tap  Binomial  transversal 
filter  is  103.24  db. 

From  observation  of  Figures  V-8  through  V-22  and  Table 
V-5,  it  is  apparent  that  the  Binomial  weighted  transversal 
filter  provides  the  best  cancellation  ability  of  any  of  the 
filters  tested.  Variations  on  the  Binomial  weighting  scheme 
were  tested,  but  with  no  increase  in  attenuation  factor. 


Chapter  VI  Conclusions  and  Recommendations 

Conclusions 

To  prudently  handle  the  clutter  cancellation  problem  a 
method  of  mathematically  describing  the  target  signal  and 
clutter  return  was  needed.  It  was  found  that  the  most 
efficient  method  for  description  was  in  the  frequency  domain, 
in  terms  of  the  power  spectral  density.  In  Chapter  II  it  was 
proven  that  the  power  spectral  density  could  be  determined  for 
a  deterministic  (i.e.  target)  and  random  (i.e.  clutter) 
signal.  Since  both  clutter  and  target  returns  were  easier  to 
describe  in  the  frequency  domain,  signal  analysis  was  also 
described  in  the  frequency  domain.  It  turned  out  that  signal 
analysis  essentially  reduced  to  designing  the  matched  filter 
which  would  accomplish  the  correlation  principle,  normally 
implemented  via  a  convolution  integral.  The  design  and 
analysis  of  the  matched  filter  were  easier  in  the  frequency 
domain  because  the  transformation  between  the  time  and 
frequency  domains,  accomplished  by  the  Fourier  transformation, 
eliminates  the  need  to  perform  convolution  integral  and 
instead  performs  simpler  algebraic  manipulations.  If  the 
clutter  environment  statistics  are  known,  then  with  the  aid  of 
the  ambiguity  function  an  "optimal"  waveform-filter  pair  could 
be  defined. 

The  survey  of  doppler  clutter  cancellation  methods 
showed  that  clutter  cancellation  is  accomplished  by 
distinguishing  clutter  returns  from  target  returns  by  the 


velocity  differential  between  the  two.  This  velocity 
differential  manifests  itself  as  power  spectra  located  at 
different  frequencies.  The  cancellation  occurs  by  eliminating 
those  spectral  components  which  correspond  to  clutter.  If  the 
environment  and  the  geometry  of  the  particular  surveillance 
situation  are  given,  the  PSD  of  the  clutter  can  be 
predetermined,  thus  allowing  for  proper  spectral  nulling.  This 
relatively  straightforward  approach  to  eliminating  clutter  is 
complicated  when  the  target  and  clutter  velocity  differential 
is  such  that  the  target  power  spectral  components  overlap  the 
spectral  components  of  clutter.  For  a  space-based  radar 
system,  the  situation  of  a  low  velocity  differential  between 
clutter  and  target  spectrum  arises  either  because  of  the 
geometry  of  the  particular  surveillance  situation  and/or 
because  of  the  radar  system's  platform  motion.  In  the  pre¬ 
vious  chapter  we  saw  that  the  greatest  clutter  spread  occurs 
at  small  depression  angles  and  large  azimuth  angles.  We  also 
observed  that  clutter  energy,  after  clutter  locking  compensa¬ 
tion,  is  centered  about  zero  Hz  and  spreads  out  to  greater 
than  400  Hz.  The  cancellation  of  clutter  in  the  low  doppler 
(velocity)  regions  is  the  essence  of  the  clutter  cancellation 
problem. 

If  the  system  compensation  methods  are  ignored,  i.e. 
DPCA  or  monopulse  anomolous  nulling,  the  clutter  cancellation 
problem  reduces  to  a  filter  design  problem.  Because  of  the 
large  area  which  the  space-based  radar  must  survey,  there  is  a 
vast  number  of  detection  cells  which  must  be  analyzed.  This 


144 


necessitates  that  the  filters  be  implemented  digitally.  The 
choice  must  be  made  whether  the  filters  should  be  implemented 
via  recursive  equation,  thereby  taking  advantage  of  the  vast 
experience  in  the  design  of  analog  filters  or  implementing  the 
filters  via  nonrecursive  equations,  which  can  be  computed 
using  one  of  the  various  FFT  algorithms.  Because  recursive 
equations  cannot  be  implemented  via  FFT  operations,  and  also 
because  the  transient  period  for  stabilization  tends  to  be 
longer  for  recursive  filters,  nonrecursive  filters  were 
utilized  for  various  clutter  cancellation  schemes. 

Since  the  transversal  filter  can  be  used  to  represent 
either  a  delay  line  canceler  or  a  filter  in  a  filter  bank, 
analyzing  the  spectral  nulling  of  different  transcient  filters 
provides  an  indication  of  the  clutter  cancellation  ability  of 
different  cancellation  schemes.  In  Chapter  V  the  spectral 
response  of  transversal  filters  with  seven  different  weighting 
schemes,  with  different  tap  delay  times,  and  with  different 
numbers  of  taps  were  all  examined  to  determine  which 
transversal  filter  configuration  provided  sufficient  clutter 
cancellation.  It  was  shown  that  a  transversal  filter  using 
the  Binomial  weighting  scheme  and  a  tap  delay  time  equal  to 
the  reciprocal  of  the  target  doppler  bandwidth  provided  the 
best  clutter  cancellation.  Increasing  the  number  of  taps 
broadened  the  filter  spectral  null,  thereby  improving  the 
clutter  cancellation  at  the  expense  of  eliminating  more 
targets.  A  four  tap  Binomial  weighted  transversal  filter 
provided  103.24  db  of  average  clutter  attenuation,  a  value 


sufficiently  adequate  to  cancel  clutter  and  obtain  target 
range  and  velocity  information. 

Recommendations 

The  primary  objective  of  this  study  was  to  organize, 
define  and  develop  the  nature  and  description  of  the 
spaceborne  radar  clutter  cancellation  problem.  Besides  the 
development  of  the  clutter  problem,  a  secondary  objective  was 
to  develop  a  computer  tool  that  would  1 )  implement  the  equa¬ 
tions  and  concepts  which  describe  the  clutter  problem  so  a 
fuller  understanding  of  the  clutter  problem  could  be  gained 
and  2)  provide  a  vehicle  to  test  various  clutter  cancellation 
methods.  This  study  concentrated  on  only  doppler  cancellation 
methods,  which  digitally  filtered  the  radar  echo  signal  by 
application  of  one  of  several  FFT  algorithms.  Although  this 
study  has  concisely  tied  together  the  clutter  cancellation 
problem,  there  are  still  many  areas  which  require  further 
research. 

In  this  study,  the  spectrum  of  the  echo  signal  was 
canceled  by  a  nonrecursive  filter,  implemented  by  FFT 
algorithms.  There  are  several  other  methods  of  estimating  the 
spectra  based  on  modeling  the  clutter  as  a  rational  function 
and  then  determining  the  coefficients  for  the  specific 
situation.  The  relative  merits  of  these  other  spectral 
estimation  techniques  should  be  explored  and  compared  in  terms 
of  the  clutter  attenuation  factor  and  operational  speed.  Ref 
34  provides  an  excellent  reference  to  initiate  the  analysis. 


The  continuing  change  in  the  clutter  statistics, 
especially  due  to  the  platform's  motion  results  in  filter 
mismatches.  The  use  of  an  adaptive  algorithm  such  as  the  one 
suggested  by  Applebaum,  Ref  1,  or  Reed  and  Brennen,  Refs  8  and 
10  to  modify  the  matched  filters  should  be  explored.  This 
effort  should  compare  the  clutter  attenuation  factor  of  each 
of  the  algorithms  with  that  of  a  fixed  matched  filter. 

Several  system  platform  motion  compensation  methods  were 
discussed,  but  not  analyzed.  The  effects  of  these  methods  on 
the  total  radar  system  and  the  gain  in  attenuation  factor 
should  be  determined. 

To  limit  the  size  of  this  study,  nondoppler  clutter 
cancellation  techniques  were  not  considered.  The  nondoppler 
cancellation  techniques  should  be  analyzed,  and  their  clutter 
reduction  ability  quantitatively  determined.  Normally,  the 
benefits  of  nondoppler  techniques  have  been  described  in  terms 
of  probability  of  detection.  A  performance  measure  similar  to 
the  clutter  attenuation  factor  should  be  developed  and  then 
the  doppler  and  nondoppler  cancellation  techniques  should  be 
compared . 

Finally,  the  simulation  of  the  clutter  cancellation 
problem  was  accomplished  through  the  application  of  the  radar 
range  equation.  This  somewhat  simplistic  approach  to  the 
problem  is  based  upon  approximations  to  the  power  levels. 
A  theoretically  more  accurate  and  thorough  approach  to  the 
modeling  of  the  clutter  problem  would  be  to  perform  a  video 
signal  simulation,  simulating  the  interaction  of  the 


transmitted  video  signal  with  the  target  and  the  clutter 
scatterers.  Although  a  video  simulation  would  be  more  complex 
and  utilize  more  computer  time,  a  finer  sensitivity  to  the 
response  of  various  filters  could  be  obtained. 
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Radar  System 

A  S-Band,  20  meter  diameter  radar  system  was  used  for 
this  study.  Figure  A-1  illustrates  an  example  of  an  arbitrary 
space-based  radar  system. 


Fig.  A-1  Radar  System  Configuration 
The  satellite,  with  the  antenna's  boresight  pointed  toward  the 
horizon,  is  assumed  at  an  altitude  of  1667  kilometers. 

Since  this  study  is  concerned  with  clutter  cancelers, 
the  radar  system  described  is  for  the  purpose  of  clutter 
analysis  and  should  not  be  taken  as  a  point  design.  The 
following  radar  parameters  chosen  for  this  study  are,  however, 
representative  of  a  system  that  could  be  deployed  to  perform 
wide  area  surveillance. 

1.  Transmitted  frequency  -  3.0  GHz  (S-Band) 
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2.  Power  /  Burst 

3.  Gain  (n  =.85) 
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4.  System  Noise  Temperature 

5.  System  Losses 

6.  Azimuth  Beamwidth 

7.  Elevation  Beamwidth 

8.  PRF 

9.  Scan  Rate 


5.0  kW 

55  DB 
o 

415  Kelvin 
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o 
.5 
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30  kHz 
o 

5  /sec 


A  simple  rectangular  pulse  waveform,  illustrated  in  Figure  A- 
2,  is  used  for  the  analysis. 
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Fig.  2-A  Pulse  Doppler  Waveform 


Clutter  Model 

Due  to  the  uncertain  nature  of  the  statistics  of  radar 
clutter  as  seen  by  a  space-based  radar  system,  the 
backscattering  coefficient,  a4  ,  described  in  Ref  4:4-24,  is 


used  for  this  study.  Figure  A-3  illustrates  the  average 


backscatter  values  for  different  terrains 


The  values  of  or*  are  for  a  S-Band  radar. 

The  computer  simulation  of  the  clutter  cancellation 
problem  allows  the  user  to  choose  either  land, sea,  or  snow  and 
ice  terrain.  Backscatter  coefficient  values,  at  every  five 
degrees  of  grazing  angle,  are  stored  for  the  three  types  of 
terrain.  The  program  stores  the  DSA  backscatter  model  values 
for  land  and  snow/ice.  The  sea  cr**  values  are  taken  from  the 
Sea  State  5  curve  for  grazing  angles  less  than  70  degrees. 
For  grazing  angles  greater  than  70  degrees  the  Aerospace 
Program  Model  is  utilized.  The  change  at  70  degrees  is 
because  the  Aerospace  Program  Model  better  approximates  the 
available  data  values. 


For  a  pulse  doppler  radar,  the  amount  of  clutter  surface 
area  Illuminated  by  the  antenna  mainbeam  is  dependent  upon  the 
platform's  altitude  above  the  earth's  surface,  the  grazing 
angle  and  the  mainbeam* s  azimuth  and  elevation  beamwidths. 
Given  a  geometric  description  of  the  mainbeam  with  respect  to 
the  earth,  the  clutter  surface  area  can  be  computed.  Figure 
B-1  describes  the  angles  used  in  describing  the  geometry 
between  the  radar  and  the  surface  area. 


DEPRESSION 
-  ANGLE 


GRAZING 

ANGLE 


INCIDENT 

ANGLE 


Fig.  B-1  Radar  Geometry  Angles 


Depending  on  the  altitude  of  the  radar  platform  certain 
assumptions  can  be  made  in  computing  the  surface  area.  For  an 
airborne  or  extremely  low  spaceborne  platform,  the  radar  is 
close  enough  to  the  ground  such  that  the  surface  illuminated 
by  the  mainbeam  is  essentially  flat.  Depending  on  the  situa¬ 
tions,  two  different  cases  occur  when  computing  the  clutter 
area.  The  pulse  limiting  case  occurs  when  the  equivalent 
physical  length  of  the  pulse  is  less  than  the  value  of  the 

elevation  beamwidth  dimension  (dimension  D  ,  Figure  B-2). 

EL 

This  case  results  when  the  relationship  between  the  half  power 
elevation  beamwidth,  the  grazing  angle  V  ,  the  range  R,  and 
the  pulse  width  T  satisfy  tan  Y  <  tan  [2©  R  /cT].  Figure  B- 


2  illustrates  the  relationship  between  the  various  parameters 


RADAR 


✓  cf/2 / 
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a.  Elevation  view  showing  surface  clutter. 


RADAR  — —  R 


—  —  1  D 


b.  Azimuth  view  showing  surface  clutter. 


Fig.  B-2  Illuminated  Clutter  Area  for  Pulse  Limiting  Case 


Dimension  D  can  be  calculated  given  the  radar's  grazing 
EL 

angle  V  ,  and  the  pulse  width  T.  The  angle  0  is  defined  as 


0  =  180  -  (90  -  Y  )  -  90  =  V 


( B1  ) 


Therefore  D  is 
EL 


cf/2  =  cos  V 
D 

EL 


D  =  cT/2  sec  VF 
EL 


For  a  narrow  azimuth  beamwidth,  the  dimension  D  is 

AZ 

D  =  R  0 
AZ  BA 


(B2) 


( B3 ) 


(B4) 


and  the  illuminated  area  is  approximated  by  multiplying  the 

width  dimension  D  ,  with  the  elevation  length,  D 

AZ  EL 

A  =  (RG  cT)/2  sec  V  (B5) 

c  BA 

The  second  case  occurs  when  the  illumiated  area  is 
limited  by  the  beamwidth  and  not  by  the  pulse  length  as  shown 
in  Figure  B-3 

RADAR 


(Ref  41 : 65 ) 

Fig.  B-3  Illuminated  Clutter  Area  for  Beamwidth  Limiting  Case 

The  beamwidth  limiting  case  occurs  when  tan  V  >  tan [ 2^R/cT] . 
The  mainbeam  clutter  is  thus 

2 

A  =  .2577R  6  G  esc  V  (B6) 

C  BA  BE 

For  equation  B5  and  B6  to  apply,  the  beamwidths  are  assumed  to 
be  less  than  ten  degrees. 

For  a  high  altitude  radar  system,  the  area  illuminated 
is  shown  in  Figure  B-4  where  h  is  the  platform  altitude,  R  is 


the  slant  range  to  the  toe  of  the  footprint,  the  angles  R  , 

n  ^ 

01  ,  i>  ,  and  LS  ,  ,  4>  respectively  correspond  to  the  heel  and 

h  h  t  t  t 

toe  of  the  footprint. 
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Fig.  B-4  Radar  -  Earth  Geometry 


To  calculate  the  area  covered  by  the  radiated  energy  the 
toe  grazing  angle  must  be  calculated.  From  the  law  of  sines 


sin  +  90)  sin« 

t  t  . 

H  +  R  R 


thus 


<X  =  sin  (R  /{H  +  R  )  cos  ) 
t  e  e  t 


( B7 ) 


Again  using  the  law  of  sines,  the  heel  grazing  angle  is 
computed 


i.  .  -.  v*  •/A-  -.K  ■ 


sin  (4  4  90) 

h 


sin  <x 


h 


H  ♦  R 


where  ot  =  a  -  9  .  Thus 

h  t  BE 
-1 

4  =  cos  ( ( H  +  R  ) /R  sin  ©<  ) 

h  e  e  h 


The 


toe  earth  central  angle  Q  ,  is 


R  =  90  -  4  -  ex 

't  t  t 


and  the  heel  earth  central  angle,  ,  is 

R  =90-4  -a*. 

h  h 
-1 

*  90  -  cos  ( (H  +  R  )/R  sin  <*  )  -  *  +0 

e  e  h  t  BE 


Thus  the  length  is 


D  *  <8  -/?>* 

EL  Mt  'h  € 


-1 

=  [cos  ( ( H  +  R  )/R  sin(<x  -  O  )  -  4  -  O 

e  e  t  BE  t  BE 


To  calcultate  the  width  of  the  footprint,  D  , 

AZ 

range  must  first  be  obtained.  Again  using  the  law  of 


( B8 ) 


( B9 ) 


(BIO) 


]R  (B11) 
e 

the  slant 
sines 


thus 


R  =  R  sin  O  /  sin  cX  (B12) 

s  e  *  t  t 

The  footprint  width  is 


D  =  R  ©  { B1 3 ) 

AZ  s  BA 

The  clutter  surface  area  is  then  the  length  of  the  footprint 
times  the  width,  or 

-1 

A  =  [cos  (  (H  +  R  )/R  sin  (c<  -  0  ))  -  <t>  -  0  ]0  RR  (B14) 

c  e  e  t  BE  t  BE  BA  e  s 

(Note:  The  area  calculated  by  equation  B14  is  slightly  larger 

than  the  actual  area  covered  by  the  mainbeam.  This  is  because 
in  calculating  the  width,  the  heel  of  the  footprint  was 
assumed  to  have  the  same  width  as  the  toe  of  the  footprint.) 


APPENDIX  C 


USEFUL  FORMULAS 


Trigonometric  Identities 


sin  ( x  ±  >■)  =  sin  ( x )  cos  (.»')  ±  cos  (x)  sin  O') 
cos  (x  ±  y)  =  cos  (x)  cos  O')  +  (*)  s'n  O') 


sin  (x)  ±  sin  (>’)  =  2  sin  [(x  +  y)! 2]  cos  [(x  +  >')/2] 
cos  (x)  +  cos  O')  =  2  cos  [(x  +  >0/2]  cos  [(x  -  >0/2] 
cos  (x)  -  cos  O0  =  2  sin  [(x  +  >0/2]  sin  [O'  -  *)/2] 

2  cos  (x)  cos  O’)  =  [cos  (*  +  >0  +  cos  (x  -  >0] 

2  sin  (x)  sin  (>')  =  [cos  (x  -  >0  -  cos  (x  +  >)] 

2  sin  (x)  cos  O')  =  [s‘n  (x  +  >)  +  sin  (x  -  >■)] 


sin  (2x)  =  2  sin  (x)  cos  (x) 
cos  (2x)  =  cos2(x)  -  sin2(x) 


sin  (x/2)  —  x  \/[i  —  cos  (x)]/2 
cos  (x/2)  =  ±  VC1  +  cos  0)3/2. 


Select  sign  for  radical 
corisistcnt  with  left  side. 


2  sin2  (x)  =  1  -  cos  (2x) 
2  cos2(x)  =  1  +  cos  (2x) 


where 


4  cos3(x)  =  3  cos  (x)  +  cos  (3x) 

A  cos  (x)  -  B  sir.  (x)  =  R  cos  (x  +  8) 

r  =  J'a^Tb1 

8-,aD  uJ  - 


+je 

e  =  cos  9  +  sin  0 

j0  -j© 
cos  9  =  e  +  e 
2 


j©  -j® 

sin  ©  =  e  +  e 

2j 

j©  -j© 
tan  ©  =  -  j  (e  -  e  ) 

]©  -d© 

e  +  e 


2je 

-i  (e  -  1) 

2D© 

e  +1 


sin  A 


sin  B 


sin  C 


2  2  2 

a  =  b  +  c  -  2bc  cos  A 


C-2 
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Fourier  Transform  Theorems 


A.  Superposition 

5>,0.  (0  - 

A  *  1 

A  *  1 

H.  Symmetry 

C(/) 

</(-/) 

C.  Time  Scaling 

£t(fl0  — 

^C7(//o) 

‘0| 

D.  Time  Shifting 

9(1  -  >o)  “* 

C(/)exp  [ -J2~ft0 ] 

E.  Frequency 
Shifting 

</(r)  exp  [-r  j2-f^n  — 

GLf-h) 

F.  Time 

Differentiation 

<r 

~ 

U2-/TG(f) 

G.  Frequency 
Differentiation 

{-jlniTffO)  — 

7-'Mm 

11.  Integration 

J 

j  g(x)  fix  — 

J  -  IC 

£^J(/)  +  J_C(/, 

I  Convolution 

<j{:)  •  l‘{>i  — 

C(/)  *(/) 

J.  Multiplication 

C(/)  •  «i/) 

o-3 


Summary  of  statistical  averages* 


Definition 


Indicated  operation 


Discrete  case 


Continuous  case 


Mean  ■»  m.  Also 


called  average, 
expectation, 
ensemble  average 

E{X] 

2  x,F(x,) 

1  •  1 

J  ^xp(x)dx 

Variance  «»  cr1 

K 

(standard 

E\{X-m,Y) 

2  ix,  -  m,YP(x,) 

I  ( x-mM)ap(x)dx 

deviation  —  0 ) 

1  -  1 

J  -« 

nth  moment 

E(X') 

2  x‘P(xt) 

1  •  1 

J  ^  x*/»(x)  dx 

nth  central  moment 

E{(X-m,Y } 

(— Oforn—  1) 

or  nth  moment 
about  the  mean 

2  (Xi  -  m.YP(x,) 

1  •  1 

J  ^  (x  -  m.Ypix)  dx 

Mean  ol  a  function 
Six) 

E{g{X)) 

i  *(*.)/»(*.) 

!•  1 

J>X*X)* 

(n  +  A)th  joint 
moment 

£{*■}'*} 

2  2  x,‘y/P(x, ,  y,) 

x'y*pix'  x'i  dz  dy 

(n  -r  A)th  joint 
central  moment 

E«X  -  m.nY -  m,Y) 

2  2  (x,  -  m.TOv  -  m,)*F(x, ,  y,) 
#•1^-1 

_(x  -  m,)*(y  -  m,Yp[x,  y)  dx  dy 

Covariance 

EUX-mMY-m,)) 

2  2  (x,  -  m.Xj';  ~  m,)F(x, ,  y,) 

1- 1  /-i 

JT.J 

(x  -  m.Xr  -  m,)/Kx,  y)  dx  dy 

-  • 

Mean  ol  a  function 
Six,  y) 

E{0(X,  Y)) 

A'  M 

2  2f(x,,Xi)F(x,,Xi) 

f- 1 

f  -  -  J—  ^ 

•  For  the  discrete  cases,  outcomes  are  x»(f  —  1 . N)  and  yt  {j  —  1, ... ,  M)  with  probability  /*(x,)  and  P(jj)  respectively.  For  the  con¬ 

tinuous  cases  X  and  Y  are  defined  over  the  interval  (—00  <x,y<  00)  and  have  continuous  density  functions  p(x)  and  p{y)  respectively. 
The  means  are  indicated  by  m,  and  m, . 


B,  W.  LimJgren,  Statist  nut  Throry,  The  Macmillan  Company,  I960. 


Z  -  Transform  Theorems 


Sequence  x-Transform 


1.  x(n) 

XU) 

A,.  <  |z|  <  A,. 

2.  y(ri) 

YU) 

A,.  <  |z|  <  A,. 

3.  ax(n)  t  by(n) 

aXU)  +  b  YU) 

max  IA,_,A,_)  <  |z|  <  mm  IA.„,  A„) 

4.  x(n  +  n,) 

fXU) 

A,.  <  |z|  <  A,. 

5.  a'x{n) 

X(a-U) 

dXU) 

|o|A,.  <  |z]  <  |o|A,» 

6.  nx(n) 

2  dz 

A,.  <  |z|  <  A,. 

7.  x‘(n) 

x*u%) 

A,.  <  |z|  <  A,* 

8.  x(-n) 

X{\  lx) 

1/A,.  <  |z|  <  1/A,. 

9.  Re  [*(*)) 

i[XU)  +  X'U‘)l 

A,.  <  |  z  |  <  A,. 

10.  Im  {jr(/»)] 

i  [XU)  -  X*Um)  1 

2J 

A_  <  |z|  <  A,. 

11.  x(n)  •  y(n) 

XU)  YU) 

m«  [A,.,  A,.J  <  |z|  <  min  [A.,,  A^.J 

12.  x(n)y(n) 

A_A,.  <  |z|  <  A..A,. 

DFT  Transform  Theorems 


Finite-Length  Sequence  ( length  N) 

1.  x(n) 

2.  y{n) 

3.  ax(n)  +  by(n) 

4.  x ((«  +  m))yS.K(n) 

5.  W';x(n) 

6.  ^  -  m))ArJ*A*W 

7.  Jr(n)jW 

8.  **(«) 

9. 

10.  Re  (*(»)] 

11.  y'lm  (x(«)] 

12.  *,,(") 

13.  •*„(«) 


DFT 

X(k)  ' 

Y{k) 

cX(k)  +  bY(k) 

wf-xw 

X«k  -  l 0A*  v(*) 

AW1W 

^  [  i/((/)).v  y((*  -  o>a-J*x<*> 

**((-*))A*v(*) 

*•(*) 

AT.,(*)  -  *IW*)).V  4  A'*((-*))a.J51.v<*) 

*.,(*)  -  J[^((*)).v  -  X*((-*)).v]*.v<*> 


Re  [AW) 
j  Im  [AW) 

The  following  properties  apply  only  wheo  x(n)  is  real: 


14.  Any  real  x(n) 


1J.  x,,(n) 
16.  x „(n) 


I  AW  -  X'((-k))y*v(lc) 

Re  l  AW)  -  Re  [  X((-k))y]X  y(k) 
lm  |AW)  -  -lm  [^((-*)).vJ^.vW 

\X(k)\  -  \X((-k))y\Zy(k) 
arg  [ X{k))  -  -arg  [X((-k))x)*y(k) 
Re  \X(k)) 

J  lm  [AW) 


Table  of  Elementary 
Z  -  Transform  Pairs 


_2_ 

3 


4 

5 


fOO 


1  '  5(k) 


o(k-i) 


u(k) 


u(k-i) 


ku(k) 


6  k  u(k) 


7  e 


-akT 


8  ke 


•akT 


9  1  -  e 


-akT 


LO  a' 


11  ka1 


12  (k  1  )a 


13  (k+l)(k+2)ak/2 


14  sin(wkT) 


15  cos(uikT) 


16  e-akTsin(ukT) 


F(z) 


.-i 


z  -  1 


(z  -  1>‘ 


z(  z  +  1  ) 
(2  -  l)3 


z  -  exp(-aT) 


zexp( -aT) 


{z  -  sxp(-aT) 


z{l  -  exp(-aT)} 
(z-l)iz-exp(-aT)} 


z  -  a 


az 


(z  -  a)' 


(z  -  a)4 


(z  -  a)- 


zsin(uT) 


z  -  2zcos(u>T)  +  1 


z{z  -  COS(U)T)} 

-  2zcos(uT)  +  1 


zexp(-aT)sinCu)T) 


z  -  2zexp(-aT)cos(uT)  +  exp(-2aT) 


17  e"akTcos(i_kT) 


zexp(-aT)cos((iiT) 


z  -  2zexp(-aT)cos(uT)  +  exp(-2aT) 


ml  Relationships 


-•'ROGRAM  CLTCAN 

¥  %  %  >K  #  #  >K  $  *  $  #  *  *.  *  Xofc##  >!■■:  *  #  #  *  :*:  *  * ».  %j|o|()|c)|<#  * 

A  PROGRAM  CL. UT CAN  * 

*  WRITTEN  BYt  JOE  DEVENUTQ  * 

*  15  SEPTEMBER  1963  * 

*  *  ************  ********  IK****##  ##*•>!■:* 

THE  PROGRAM  CLUTCAN  IS  AN  INITIALIZATION  AND 
CALLING  PROGRAM  THAT  SIMULATES  THE  CLUTTER  CAN¬ 
CELLATION  PROBLEM,  CLUTCAN  CALLS  THE  SUBPROGRAM  GLUT 
TO  CALCUALTE  THE  CLUTTER  POWER  SPECTRAL  DENSITY  AS  SEEN  BY 
A  SPACEBASED  DOPPLER  RADAR  SYSTEM.  CLUTCAN  THEM  CALLS  THE 
SUBPROGRAM  CANCEL,  A  PROGRAM  WHICH  SIMULATES  THE  VARIOUS 
CANCELLATION  SCHEMES,  TO  FILTER  THE  CLUTTER  PSD  THAT  WAS 
CALCUALTED  BY  CLUT. CLUTCAN  THEN  CALLS  THE  SUBPROGRAM  CLTATT 
TO  CALCULATE  THE  CLUTTER  ATTENUATION  FACTOR  AS  A  RESULT  OF 
THE  CANCELLATION  SCHEME  SIMULATED  BY  CANCEL. 

CLUTCAN  CALLS  THE  FOLLOWING  SUBPROGRAMS 
CLUT 
CANCEL 
CLTATT 


THE  FOLLOWING  VARIABLES  ARE  USED  IN  ’CLUTCAN’ 

INPUT  RADAR  PARAMETERS 

FREG-RADAR  OPERATING  FREQUENCY  (HZ) 

PRF-PULSE  REPETITION  FREQUENCY  (HZ) 

TAU-PULSE  WIDTH  (SEC)  i 

PT-PEAK  POWER  TRANSMITTED  PER  BURST  (WATTS) 
BURST-BURST  WIDTH  (SEC) 

BURPRF-BURST  PRF  (HZ) 

LS-SYSTEM  LOSSES 

NT-SYSTEM  NOISE  TEMPERATURE  (DEGREES  KELVIN) 
CANT-INPUT  TO  CONTROL  ANTENNA  FUNCTION 

1- COMPUTE  ANTENNA  GAIN 

2- FIXED  ANTENNA  GAIN 

GMAX-MAXIMUM  1-WAY  GAIN  OF  ANTENNA  MAINBEAM  (DB) 
GMIN-MINIMUM  1-UAY  GAIN  OF  ANTENNA  SIDELOBES  (DB) 
GCK1 -CONSTANT  TO  COMPUTE  1-WAY  ANTENNA  GAIN  BETWEEN 
ANGMB  AND  ANGSL 

GCK2-C0NSTANT  TO  CONFUTE  ANTENNA  1-WAY  GAIN  BETWEEN 
ANGMB  AND  ANGSL 

GCK3-ANGLE  SCALE  FACTOR  TO  COMPUTE  1-WAY  ANT.  GAIN 
CN-EXPONENT  TO  COMPUTE  ANTENNA  1-WAY  GAIN  BETWEEN 
ANGMB  AND  ANGSL 

ANGSL-RADIAL  ANGLE  (DEG)  BEYOND  WHICH  MIN  SIDE- 
LOBE  GAIN  IS  COMPUTED 

CMB-ANTENNA  ANGLE  SCALE  FACTOR  FOR  COMPUTING  MAIN 
BEAM  GAIN 

BEL-ANTENNA  ELEVATION  BEAMWIDTH  (DEG) 

BAZ-ANTENNA  AZIMUTH  BEAMWIDTH  (DEG) 

D-ANTENMA  DIAMETER  (M) 

ALPHA-SCAN  RATE  (DEG/SEC) 

OUTPUT  VARIABLES 


KRAI)  --RADAR  POWER  CONSTANT  GIVEN  PT , PRF . BAZ . ALPHA  <W> 

I N  P  U  T  G  E  0  M  E  T  R  Y  P  .A  R  A  M  E  T  E  R  S 
H - S  A T  E L L I T  E  ALTITUDE  (KM  > 

A NT A Z~ AZIMUTH  ANGLE  OF  ANT,  MAJNBEAM  FORESIGHT  (DEG) 
DEPANG- DEPRESSION  ANGLE  OF  Mi  A  INBEAM  FORESIGHT  (DEG) 

INPUT  TARGET  PARAMETERS 

VT-TARGET  VELOCITY  (M/SEC) 

R T - T  A R G ET  R A N G E  (KM ) 

RCS-TARGET  RADAR  CROSS  SECTION  <M*M) 

CLUTTER  PARAMETERS 
INPUT  VARIABLES 

CCLUT-I NPUT  TO  CONTROL  CLUTTER  OPTION 

1 - ASSUME  GAUSSIAN  CLUTTER  SPECTRUM 

2- COMPUTE  CLUTTER  SPECTRUM 
CTER-INPUT  TO  CONTROL  TERRAIN  TYPE 

1- LAND 

2- SEA 

3- SNOW  AND  ICE 

CSGPDF-INPUT  TO  CONTROL  PDF  OF  SIGMA-0 

1- CONSTANT 

2- GAUSSI AN 

3- LOG-NORMAL 
OUTPUT  VARIABLES 

CLTPSD(N) -ARRAY  OF  CLUTTER  PSD  VALUES 

DELFQ -FREQUENCY  INTERVAL  BETWEEN  FREQ  SPECTRAL 


a 


CANCELLATION  PARAMETERS 
INPUT  VARIABLES 

CCNSCH-INPUT  TO  CONTROL  CANCELLATION  SCHEME 

1- BARTLETT  WEIGHTS 

2- HANNING  WEIGHTS 

3- HAMMING  WEIGHTS 

4- BLACKMAN  WEIGHTS 

5- BINOMIAL  WEIGHTS 

6-  IMPULSE  RESPONSE 

7- UNIFORM  WEIGHTS 

8- BUTTERWORTH  WEIGHTS 
N-NUMBER  OF  POINTS 
T-SAMPLING  INTERVAL 
TAP-NUMBER  OF  FILiER  TAPS 

OUTPUT  VARIABLES 

FLCPSD ( N ) -ARRAY  OF  FILTERED  PSD  COMPONENTS 
FLMAG ( N ) -ARRAY  OF  FILTER  MAG  COEFFICIENTS 


ATTENUATION  FACTOR  PARAMETERS 
CTAT-CLUTTER  ATTENUATION  FACTOR 


1 


CLUTCAN  PROGRAM  VARIABLES 

CTEST-INPUT  TO  CONTROL  TEST  VARIATION  OF  CLUTTER 
CANCELLATION  PROBLEM 

1 - NO  MORE  TESTING 

2- TEST  WITH  DIFFERENT  CANCELLATION  SCHEME 


o  n  n  n  o  o 


C 

c 

c 

c 

c 


c 

c 

c 


3- TEST  WITH  DIFFERENT  CLUTTER  PARAMETERS 

4 - TEST  WITH  DIFFERENT  CLUTTER  PARAMETERS  AND 
C A N CELL A T  T  0 N  S C H E M E 

DATANM-NAME  OF  FILE  RESULTS  WRITTEN  TO 


THE  FOLLOWING  CONSTANTS  ARE  USED  BY  THE  SUBPR3GR 
SUBROUTINES  CLUTCAN  CALLS  DIRECTLY  OR  INDIRECTLY 
R  E  -  E  A  R  T  H  R  A  I>  I U  S  <  K  M ) 

G R A  V - E A R T H  G R A V ITATI 0 N A L  C 0 N S T A N T  ( M / S E C * S E C  > 
CLGHT-VEL 0 C I T Y  OF  LIGHT  CM/ SEC) 

BK-BOLTZMAN ' S  CONSTANT 

DTR-DEGREES  TO  RADIANS  CONVERSION  FACTOR 
RTD-RADIANS  TO  DEGREES  CONVERSION  FACTOR 


AND 


PI-PI 

KMTNM- KILOMETER  TO  NAUTICAL  MILES  CONVERSION  FACTOR 
NMTKM-NAUTICAL  MILES  TO  KILOMETERS  CONVERSION  FACTOR 


♦DECLARE  AND  INITIALIZE  VARIABLES* 

DIMENSION  FLCPSDC 150) , CLTPSDC 150) ,FLMAG<150) ,CMAG< 150) , 
+  FMAGC150) 

REAL  FREQ ,  F'T  ,  LS  ,  NT  ,  GMAX  ,  GMIN ,  GCK1 ,  GCK2  r  GCK3  ,  BEL  ,  BAZ  ,  D  , 

+  ALPHA , H , ANTAZ , DEPANG , VT , RT  , PCS , T , CT AT , KRAD  , RE , 

+  GRAV  ,  CLGHT  ,  BK ,  DTR ,  RTD ,  PI ,  KMTNM  ,  NMTKM ,  FLCF’SD , 

+  CLTPSD,FLMAG,DELFQ,CN,ANGSL,CMB,TAU,PRF, BURST 
INTEGER  CANT  ,  CCLUT ,  CTER  ,  CSGF'DF  ,  CCNSCH, N  ,  CTEST , I , J , 

+  K ,  L  ,  M ,  BURF'RF  ,  TAP 
CHARACTER  DAT ANM*6 

COMMON  RE  , GRAV , CLGHT , BK , DTR , RTD  , PI , KMTNM , NMTKM 
DATA  FREO , PT , LS , NT , GMAX , GMIN , GCK 1 , GCK2 , GCK3 , BEL  F BAZ , D , 
ALPHA , H , ANTAZ , DEPANG , VT , RT , PCS , T  T  CTAT , KRAD , PRF  f TAU , 
DELFO , BURST/26*0 . 0/ , I , J , K , L  , M , CANT , CCLUT , CTER , 

CSGF'DF  ,  CCN3CH  ,  N  ,  CTEST  ,  BURF'RF  ,  TAP/ 14*0/ 

DO  10  1=1,150,1 
FLCPSDC I )=0,0 
CLTPSDC I )=0.0 
FLMAG  C I ) =0  *  0 
CONTINUE 

*SET  VALUES  OF  CONSTANTS* 

PI =3. 14159265 
DTR=PI/180 
RTD  =  180/F'I 
CLGHT =2 » 997925E08 
RE=6378  *  388 
GRAV=9. 80665 
BK=1 . 38E-23 
NMTKM=1 .853 
KMTNM= , 5396 

♦INPUT  DATA  VALUES* 

PRINT*, ' INPUT  RADAR  OPERATING  FREQUENCY  (HZ)' 

READ*, FREQ 

PRINT*, ' INPUT  PULSE  REPETITION  FREQUENCY  (HZ)' 

READ*, PRF 
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I  riD-A128  46?  CLUTTER  CANCELLATION  TECHNIQUES  FOR  USE  IN  A  3/3 

SPACE-BASED  RADAR  SVSTEMCU)  AIR  FORC/i  INST  OF  TECH  ** 
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PRINT* * ' INPUT  PULSE  WIDTH' 

READ*  r T AU 

PRINT*,  '  INPUT  PEAK  BURST  TRANSMITTED  POWER  twFTS!' 

READ*, FT 

PRINT*, 'INPUT  BURST  WIDTH  (SEC)  ' 

READ*. BURST 

PRINT*, ' INPUT  BURST  FRF  (HZ)' 

READ* , BURPRF 

PRINT*, ' INPUT  SYSTEM  LOSSES  (DB) ' 

READ* , LS 

PRINT*, 'INPUT  SYSTEM  NOISE  TEMPERATURE  <  DEGREES  KELVIN ) ' 

READ*, NT 

PRINT*, 'INPUT  ANTENNA  CONTROL  FUNCTION  (1-COMPUTE  ANTENNA 
+GAIN ,  2-FIXED  ANTENNA  GAIN)' 

READ*, CANT 

PRINT*, ' INPUT  MAX  1-WAY  GAIN  OF  ANTENNA  MAINBEAM  (DB)' 
READ* , GMAX 

PRINT*, 'INPUT  MIN  1-WAY  GAIN  OF  ANTENNA  SIDELOBES  (DB)' 
READ* , GMIN 

PRINT*, ' INPUT  CONSTANT  TO  COMPUTE  ANTENNA  1-WAY  GAIN 
+( RATIO)  BETWEEN  MBANG  AND  SLANG' 

READ* , GCK1 

PRINT*, ' INPUT  CONSTANT  TO  COMPUTE  ANTENNA  1-WAY  GAIN 
+( RATIO)  BETWEEN  MBANG  AND  SLANG' 

READ* , GCK2 

PRINT*, ' INPUT  ANGLE  SCALE  FACTOR  TO  COMPUTE  1-WAY 
+GAIN  (DB)  BETWEEN  MBANG  AND  SLANG' 

READ* , GCK3 

PRINT*, 'INPUT  ANT. ANGLE  SCALE  FACTOR  TO  COMPUTE  MB  GAIN' 
READ* , CMB 

PRINT*, 'INPUT  EXPONENT  TO  COMPUTE  ANT  1-WAY  GAIN  BETWEEN 
+ANGMB  AND  ANGSL ' 

READ* , CN 

PRINT*, ' INPUT  RADIAL  ANGLE  BEYOND  WHICH  MIN  SIDELOBE  GAIN 
+IS  COMPUTED' 

READ*, ANGSL 

PRINT*,' INPUT  ANTENNA  ELEVATION  BEAMWIDTH  (DEG)' 

READ*, BEL 

PRINT*, 'INPUT  ANTENNA  AZIMUTH  BEAMWIDTH  (DEG)' 

READ* , BAZ 

PRINT*, 'INPUT  ANTENNA  DIAMETER  <M>' 

READ* ,  D 

PRINT*, 'INPUT  ANTENNA  SCAN  RATE  (DEG/SEC)' 

READ*, ALPHA 

PRINT*, 'INPUT  TARGET  VELOCITY  (M/SEC)' 

READ* , VT 

PRINT*, 'INPUT  RANGE  TO  TARGET  (KM)' 

READ* ,RT 

PRINT*, ' INPUT  TARGET  RADAR  CROSS  SECTION (DB)' 

READ* , RCS 

PRINT*, 'INPUT  SATELLITE  ALTITUDE  (KM)' 

READ* ,  H 

PRINT*, ' INPUT  AZ  ANGLE  OF  ANTENNA  BORESTGHT  (DEG)' 

READ* , ANT  AZ 

PFrt NT*, 'INPUT  DEPRESSION  ANGLE  OF  ANTENNA  BORESIGHT  (DEO 


GAUSSIAN 


►  * 


40 


c 

c 


50 


READ*, DEPANG 

=RINT*, 'INPUT  CLUTTER  CALCULATION  METHOD  <1 -ASSUME 
+ CL UTTER  SPECTRUM,  2-COMPUTE  CLUTTER  SPECTRUM) ' 
READ* ,  CCL.UT 


PRINT*, 'INPUT  TERRAIN  TYPE  (1-LAND,  2-SEA 
READ* , CTER 

PRINT*, ' INPUT  PDF  OF  SIGMA-0  ( 1 -CONSTANT , 
+3-L0G-N0RMAL )  ' 


3-SNOW  &  ICE)' 
2-GAUSS I  AN , 


READ* ,  CSGF'DF 

PRINT*, ' INPUT  CANCELLATION  SCHEME  ( 1 -BARTLETT  WEIGHTS, 
t2-HANNING  WEIGHTS,  3-HAMMING  WEIGHTS,  4-BLACKMAN 
+WEIGHTS ,  5-BINOMIAL  WEIGHTS,  6-IMPULSE  FUNCTION, 
+7-UNIF0RM  WEIGHTS,  8 -BUTTER WORTH  WEIGHTS)' 

READ* , CCNSCH 

PRINT*, ' INPUT  NUMBER  OF  SAMPLE  POINTS' 

READ* ,  N 

PRINT*, 'INPUT  SAMPLING  INTERVAL  (SEC)' 

READ* , T 

PRINT*, 'INPUT  NUMBER  OF  FILTER  TAPS' 

READ*, TAP 

PRINT*, ' INPUT  FILE  NAME  RESULTS  ARE  WRITTEN  TO' 

READ  ( * , 35 ) DAT  ANM 
FORMAT (AA) 


OPEN <  5 , FILE =D AT ANM) 


♦COMPUTE  CLUTTER  POWER  SPECTRAL  DENSITY* 

CALL  CLUT ( F RED , PRF , TAU , PT , LS , CANT , GMAX , GMIN , CN , ANGSL , CMB , 
+  GCK1 ,GCK2,GCK3, BEL, BAZ,D, ALPHA, H, ANTAZv DEPANG, CCLUT , 

+  CTER ,  N ,  T  ,  CSGF'DF  ,  CLTPSD ,  KRAD ',  DELFR  ,  VT  ,  RT  ,  BURST  ,  BURF'RF  ) 


CALL  CANCEL ( CCNSCH , N , T , BURST , CLTPSD , FLCPSD , FLMAG , TAP ) 


♦CALCULATE  ATTENUATION  FACTOR* 

CALL  CLTATT ( CLTPSD , FLCPSD , N , CTAT) 

♦OUTPUT  STATEMENTS* 

WRITE(*, 100) 

DO  50  1=1, N,1 

WRITE (*, 150) CLTPSD ( I ) , FLCPSD ( I ) , FLMAG ( I ) 
CONTINUE 
SCL=CLTPSD ( 1 ) 

WRITE (*, 120 )DELFQ, CTAT 

PRINT*, 'INPUT  FILTER  SCALE  FACTOR' 

READ* , FSCL 

PRINT*, 'NORMALIZED  VALUES' 

WRITE (*,100) 

DO  60  1=1, N,1 


I 

\ 

i 


h 


60 


FLMAG ( I )=FLMAG( I >/FSCL 
CLTPSD  ( I )  =CLTF’SD  ( I )  /SCL 
FLCPSD ( I ) =FLCPSD < I ) / ( SCL*FSCL ) 

WRITE <*, 150) CLTPSD (I ) , FLCPSD ( I ), FLMAG (I ) 
WRITE (5,150) CLTPSD ( I > , FLCPSD ( I ) , FLMAG ( I ) 
CONTINUE 
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.  a  A  -  ■  . 


FORMAT  <  5X//1 1HCLUTTEF:  PSD  ,  1  OX  ,  20HF JETER? 

3X.26HFILTER  TRANSFORM  MAGNITUDE/) 
FORMAT  (5X//21H FREQUENCY  INTERVAL  *  ,E17 
+  26HATTENUATI0N  FACTOR  (DEO  =  ,£17.10) 
FORMAT  (E17. 1 0 , 4X ,  E17 .  .1 0 , 4X  ,  E17 . 1 0  > 


)  L  U  T  T  E 


*TEST  DIFFERENT  SENERIO* 

PRINTS,  ■'  (  "TO  TEST  A  DIFFERENT  CASE  It 
+  "2-TEST  DIFFERENT  CANCELLATION  SCf 

t  DIFFERENT  CLUTTER  STATISTICS  V 4-' 

+  CLUTTER  AND  CANCELLATION  SCHEME") 
READ*  r CTEST 
IF  (CTEST. ER.l)  THEN 
GO  TO  200 
ELSE 

IF  (CTEST. EQ. 2)  THEN 
GO  TO  40 
ELSE 

IF  (CTEST. EQ. 3)  THEN 
GO  TO  30 
ELSE 

GO  TO  20 
END  IF 
END  IF 
END  IF 


JT  "  / "  1  -  END 


TEREN 


200  STOP 


l»! 


FLOWCHART  FOR  SUBPROGRAM  CLUT 


CALCULATION) 


RETURN 


cr* 

=  «hj2  t r  © 

s 

AZ 

a* 

=  ( . 6/D)RV  COStf  COS© 

P 

s  d 

/  IF  \ 

CTER  =  2  > - T- 


*BARLOW'S  CONSTANT* 
A  =  1 .5E16 


a*  *  thJlh 

i 

_ 

2  2 

2  2 

a*  =  cr*  + 

O'*  +  a* 

c  i 

P  s 

• 

CALL  SIGMAO 

TO  DETERMINE  cr* 

CALL  CLTARA 

TO  COMPUTE  CLUTTER  AREA 

2  . 

K  =  K 

G  cr*A 

CLTPW  RAD  c 

4 

R 

s 

I  <  N-1 


2  2 

C  =  K  EXP[  -  ( Itf)  /2a*  ) 

PSD  CLTPW  C 


RETURN 


FLOWCHART  FOR  METHOD  2  CLUTTER  PSD  CALCULATION 


♦SEQUENCE  THROUGH  ALL  ISODOPS* 


CV  =  (((KD-I)PRF  +  .  5BW)  X. )  /2V 
R  S 


CV  =  (  ( (KD-1  jPRF)  X)/2V 
M  S 

CV  =  (({KD-DPRF  -  .5BW)\)/2V 
L  S 


♦FIND  MINIMUM  RANGE  TO  ISODOP  CENTER* 

CALL  TGECSR  TO  COMPUTE  GRAZING  ANGLE, 

EARTH  CENTRAL  ANGLE,  AND  SLANT  RANGE 


♦FIND  MINIMUM  RANGE  TO  ISODOP  MINIMUM* 

CALL  TGECSR  TO  COMPUTE  GRAZING  ANGLE, 

EARTH  CENTRAL  ANGLE,  AND  SLANT  RANGE 


♦SEQUENCE  THROUGH  ALL  ISORANGE  FOR  ISODOP* 


4- 

,1 


5 U B P G UTINE  C L U '  ( F F; EOS  ,  P F: F 9  v  T AUS  ,  ? T 5 .  L. S S  C A rJ 7  5  . 


O' AXS  ,  GMINS  ,  CN S  ,  A N G S L S ,  C M B S  ,  G C K 1 S  ,  G C h 
DS ,  AL  PHAS  ,  HS  ,  ANTAZS  ,  DF  ANGS  ,  CCL-JTE  ,  CTE 
CSGPDS  v  CL7PSS  ,  KRAPS  ,  DELFCS  ,  VIS  ,  R'TS ,  PL 


2$ , GCK3S , PELS . BAZS» 
RS.NSf 7St 
RSS , BURPRS > 


*  SUBPROGRAM  CLUT  * 

*  WRITTEN  BY! JOE  DEUENIP'n  * 

*  15  SEPTEMBER  1983  * 

*  *  ****  *  ********  *  *#*##*#  *  *****  * 


THE  SUBPROGRAM  CLUT  CALCULATES  THE  POWER  SPEC¬ 
TRAL  DENSITY  AS  SEEN  BY  A  SPACE-BASED  DOPPLER  RADAR 
SYSTEM  GIVEN  THE  RADAR  SYSTEM  PARAMETERS  AND  THE 
VIEWING  GEOMETRY  PARAMETERS.  THIS  SUBPROGRAM 
CALCULATES  THE  CLUTTER  PSD  BY  EITHER  ASSUMING  A 
SPECIFIED  CLUTTER  PSD  PDF  OR  BY  SUMMING  THE  DIS¬ 
TRIBUTED  CLUTTER  RETURNS  FORM  AMBIGUOUS  RANGE  BINS 
ALONG  A  GIVEN  UNAMBIGUOUS  DOPPLER  (ISODOP).  SUB¬ 
PROGRAM  CLUT  RETURNS  TO  CLTCAN  THE  CLUTTER  POWER 
FOR  THE  ANALOG  FREQUENCY  AXIS, 


CLUT  CALLS  THE  FOLLOWING  SUGROUTINES 

CLTARA-CALCULATES  THE  AREA  OF  CLUTTER 
ANTGN-CALCULATES  THE  ANTENNA  GAIN 
SIGMAO-CALCULATES  THE  SIGMA-0  VALUE 
TGECSR-CALCULATES  THE  GRAZING  ANGLE, EARTH 
CENTRAL  ANGLE  AND  SLANT  RANGE  ! 


THE  FOLLOWING  VARIABLES  ARE  USED  IN  'CLUT’ 

INPUT  VARIABLES 

FREQS-RADAR  OPERATING  FREQUENCY  (HZ) 

PRFS-PULSE  REPETITION  FREQUENCY  (HZ) 

TAUS-PULSE  WIDTH  (SEC) 

PTS-PEAK  POWER  TRANSMITTED  PER  PULSE  (WATTS) 
LSS-SYSTEM  LOSSES 

CANTS~INPL.IT  TO  CONTROL  ANTENNA  FUNCTION 
GMAXS-MAX  1-WAY  GAIN  OF  ANTENNA  MAINBEAM  (DEG) 
GMINS-MIN  1-WAY  GAIN  OF  ANTENNA  SIDELOBES  (DEG) 
GCK1S , GCK2S, GCK3S- ANTENNA  CONSTANTS 
CMBS-ANTENNA  ANGLE  SCALE  FACTOR  FOR  COMPUTING  MAIN 


BEAM  GAIN 

CNS-EXPONENT  TO  COMPUTE  ANTENNA  1-WAY  GAIN  BETWEEN 
ANGMB  AND  ANGSL 

ANGSLS-RADI AL  ANGLE  (DEG)  BEYOND  WHICH  MIN  SIDE- 
LOBE  GAIN  IS  COMPUTED 
BELS-ANTENNA  ELEVATION  BEAMWIDTH  (DEG) 

BAZS-ANTENNA  AZIMUTH  BEAMWIDTH  (DEG) 

DS-ANTENNA  DIAMETER 
ALPHAS-SCAN  RATE  (DEG/SEC) 


.V, 


-V- 


*5* 


L, 

C 

C 

r 


r 


t- 

C 

C 

C 

c 

c 

c 

c 


c 

C- 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DEG) 

:GH"  (DEG) 
OPTION 


RTS-RANGE  TO  TARGET  (KM) 

HS--SATELL.ITE  ALTITUDE  (KM) 

ANT A2S-AZ I MUTH  ANGLE  OF  ANTENNA  MA INBEAM  > 
DPANGS-DEF'RESSION  ANGLE  OF  MA  INBEAM  BORES) 

CCLUTS- INPUT  TO  CONTROL  CLUTTER  COMPUTATIC 
CTERS- INPUT  TO  CONTROL  TERRAIN  TYPE 
CSGPDS- INPUT  TO  CONTROL  PDF  OF  SIGMA-0 
BURSS-BURST  WIDTH  (SEC) 

BURPRS-BURST  PRF  (HZ) 

OUTPUT  VARIABLES 

CLTF'SS  ( N )  -ARRAY  OF  CLUTTER  PSD  VALUES 
KRADS-RADAR  POWER  CONSTANT  GIVEN  PT , PRF r BAZ , ALPHA ( W ) 
DELFQS -FREQUENCY  INTERVAL  BETWEEN  FREG  SPECTRA 
SUBPROGRAM  VARIABLES 
DW-DWELL  TIME 
V5-VEL0CITY  OF  SATELLITE 
ARACLT-AREA  OF  CLUTTER 
THETAG-GRAZING  ANGLE 


c 

GN1 -ANTENNA 

GAIN 

WRT 

CLUTTER 

PATCH 

IN 

1ST 

QUAD 

(DB) 

c 

GN2-ANTENNA 

GAIN 

WRT 

CLUTTER 

PATCH 

IN 

2ND. 

QUAD 

(DB) 

c 

GN3-ANTENNA 

GAIN 

WRT 

CLUTTER 

PATCH 

IN 

3RD 

QUAD 

(DB) 

Ur 

GN4-ANTENNA 

GAIN 

WRT 

CLUTTER 

PATCH 

IN 

4TK 

QUAD 

( DB ) 

DB 

DB 


RPV-RADIAL  PLATFORM  VELOCITY 
LAMDA-WAVELENGTH 
ECANG-EARTH  CENTRAL  ANGLE 
RS-SLANT  RANGE 

LSNDB-SYSTEM  LOSSES  NOT  IN  DB 
RCSNDB-RCS  NOT  IN  DB 

GMXNDB-l-WAY  MAX  ANTENNA  GAIN  NOT  IN 
GMNNDB-1 -WAY  MIN  ANTENNA  GAIN' NOT  IN 
STDSCN-CLUTTER  SPREAD  DUE  TO  SCANNING 
STDPTM-CL. UTTER  SPREAD  DUE  TO  PLATFORM  MOTION 
STDIM-CLUTTER  SPREAD  DUE  TO  INTERNAL  MOTION 
KCLT VAR- TOTAL  CLUTTER  SPREAD 
KCLTPW-CLUTTER  POWER  CGNSTANT 
A-BARLOW'S  INTERNAL  MOTION  CONSTANT 
SIGMA-CLUTTER  SIGMA-0  VALUE 
ANG-ANGLE  (DEG) 

DELR-DELTA  RANGE 

FT-TARGET  DOPPLER  FREQUENCY  (HZ) 

BWF-FILTER  BANDWIDTH  (HZ) 

NAHOR-NADIR  ANGLE  TO  HORIZON  (DEG) 

DEF’HOR-DEPRESSION*  ANGLE  TO  HORIZON  (DEG) 

ECAHOR-EARTH  CENTRAL  ANGLE  TO  HORIZON  (DEG) 

RSHOR-SLANT  RANGE  TO  HORIZON  (KM) 

KD-LOOP  VARIABLE  FOR  ISODOPS 
KR-LOOP  VARIABLE  FOR  ISORANGE 
NRM AX -CONTROL  FOR  MAXIMUM  ISORANGES 
KDMAX-CONTROL  FOR  MAXIMUM  ISODOPS 
CVR-RIGHT  CLUTTER  DOPPLER  LINE 
CVM-MIDDLE  CLUTTER  DOPPLER  LINE 
CVl.-LEFT  CLUTTER  DOPPLER  LINE 

CAANG-CLUTTER  AZIMUTH  ANGLE  TO  CENTER  OF  CLUTTER  PATCH 
CNANG-NADIF:  ANGLE  TO  CENTER  OF  CLUTTER  PATCH  (DEG) 

KRM IN-CONTROL  FOR  MINIMUM  ISORANGE 
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*  W  V  *: 


•  •  ■  'j*  •  .  ■ 


«V.\. «v 


h 


9 


a 


RMIND-MINIM'JM  RANGE  TO  ISODOP  (KM) 

RMINDN-MINIMUM  RANGE  TO  MINIMUM  ISODOP  (KM) 
CDEPA-DEPRESSION  ANGLE  TO  CENTER  OF  CLUTTER  PATCH 
RG -GROUND  RANGE  OF  CLUTTER  PATCH  (KM) 

RSF-SLANT  RANGE  TO  FAR  EDGE  OF  CLT  PATCH  (KM) 

RSN- SLANT  RANGE  TO  NEAR  EDGE  OF  CLT  PATCH  KM) 

NAF -NADIR  ANGLE  TO  FAR  EDGE  OF  CLT  PATCH  (DEG) 

NAN-NADIR  ANGLE  TO  NEAR  EDGE  OF  CLT  PATCH  (DEG) 
DEPRF-DEPRESSION  ANGLE  TO  FAR  EDGE  OF  CLT  PATCH  (DEG) 
DEF'RN-DEPRESSION  ANGLE  TO  NEAR  EDGE  OF  CLT  PATCH  (DEG) 
ECANGF-EARTH  CENTRAL  ANG  TO  FAR  EDGE  OF  CLT  PATCH (DEG) 
ECANGN-EARTH  CENTRAL  ANG  TO  NEAR  EDGE  OF  CLT  PATCH (DEG) 
A ANGFL- AZIMUTH  ANG  OF  FAR  LEFT  EDGE  OF  CLT  PATCH  (DEG) 

A ANGFR- AZIMUTH  ANG  OF  FAR  RIGHT  EDGE  OF  CLT  PATCH (DEG) 

A ANGNL- AZIMUTH  ANG  OF  NEAR  LEFT  EDGE  OF  CLT  PATCH(DEG) 

A ANGNR- AZIMUTH  ANG  OF  NEAR  RIGHT  EDGE  OF  CLT  PATCH (DEG) 
BF-LENGTH  OF  FAR  EDGE  OF  CLUTTER  PATCH  (KM) 

BN-LENGTH  OF  NEAR  EDGE  OF  CLUTTER  PATCH  (KM) 


♦DECLARE  AND  INITIALIZE  VARIABLES* 

DIMENSION  CLTPSS ( 150 ) 

REAL  FREDS , PTS , LSS , GMAXS , GMINS , BELS , BAZ5 , DS , ALPHAS, 

+  HS , ANTAZS , DPANGS , CLTPSS , KRADS , DELFQS , GCK1 S , 

+  GCK2S , GCK3S , DU , VS , ARACLT , THETAG , GN1 , GN2 , GN3 , GN4 , 

+  RPV , LAMDA , ECANG , RS , GMXNDB , GMNNDB , LSNDB , RCSNDB , 

+  STDSCN ,  STDPTM ,  STDIM ,  KCLTF’W ,  A ,  CLTVAR , SIGMA ,  TAUS ,  PRFS , 

+  BURSS  - 

REAL  ANG , DELR , FT , BUF , NAHOR , DEPHOR ,  ECAHOR , RSHOR , 

+  CVR , CVM , CVL , CAANG , CNANG , CDEPA , RMIND , RMINDN , 

+  RG , RSF , RSN , NAF , NAN , DEPRF , DEPRN , ECANGF , ECANGN , 

+  AANGFL , AANGFR , AANGNL , AANGNR , BF , BN , CNS , ANGSLS , CMBS 

INTEGER  CANTS,  CCLUTS , CTERS , CSGPDS ,  If  J,L., 

+  BURPRS , KD , KR , KRMAX , KDMAX , KRMIN 

COMMON  RE , GRAV , CLGHT , BK , DTR , RTD , PI , KMTNM , NMTKM 
1=0 
J=0 
L=0 

DO  205  1=1,150,1 
CLTPSS ( I ) =0 . 0 
CONTINUE 


♦INITIAL  CALCULATIONS* 

♦COMPUTE  PLATFORM  RADIAL  VELOCITY* 
VS=SQRT(3.986012E05/(HS+RE)  ) 

RPV=VS* ( COS ( ANTAZS*DTR ) ) * ( COS ( DPANGS*DTR ) ) 
♦  CO*  F’UTE  FREQ  SAMPLING  RATE* 

P'  .r  QS=PRFS/NS 

♦COMPUTE  RADAR  POUER  CONSTANT* 

LSNDB=10 ♦ ** ( LSS* ♦ 1 ) 

DW=BAZS/ ALPHAS 

LAMDA=CLGHT/FREQS 

KRADS= ( PTS*DW*BURPRS*LAMDA*LAMDA ) 

+  /( ( (4*PI)**3)*LSNDB) 
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* COMPUTE  SLANT  RANGE* 

CALl.  TGECSR  <  HS  -  DPANGS  ,  THETAG  ,  ECANG ,  RS ) 

* DETERMINE  METHOD  OF  CLUTTER  PSD  CALCULATION* 

IF  < CCLL'TS  ♦  EQ »  2)  GO  TO  250 

♦COMPUTE  CLUTTER  PSD  VIA  GAUSSIAN  SPECTRUM* 

♦COMPUTE  CLUTTER  SPREAD* 

STDSCN= ALPHAS/ < ( 2** . 5 ) *PI*BAZS ) 

STDPTM= ( *  6/DS ) *VS*COS ( DPANGS*DTR > *  1 000 
r  *SIN( ANT AZS*PTR ) 

IF  (CTERS  .EQ.  2)  THEN 
A=1 . 5E16 
ELSE 

A=2 . 3E17 
END  IF 

STDIM=FREQS/< (2*A)**.5) 

CLTVAR= ( STDSCN**2 )  +  <  STDPTM**2 )  +  <  STDIM**2 ) 

♦COMPUTE  CLUTTER  PSD* 

GMXNDB= 1 0 . ** ( GMAXS* . 1 ) 

CALL  SI GMAO( THETAG f CTERS, CSGPDS, SIGMA) 
SIGMA=10.**(SIGMA*.l) 

CALL  CLTARA ( HS , BAZS , BELS , THETAG , RS , ARACLT ) 

KCLTPU= ( KRADS*GMXNDB*GMXNDB*ARACLT*SIGMA*1E06 ) / 

+  < <RS*1000)**4) 

DO  210  1=0 , NS~1 , 1 

ARG= < ( I *DELFQS ) **2 ) / ( 2*CLT VAR ) 

IF  (ARG.GT.250, )  THEN 

CLTPSS < 1+1 )  =0.0  ' 

GO  TO  210 
END  IF 

CLTPSS ( 1+1 ) =KCLTPW* ( EXP ( -ARG ) ) 

CONTINUE 
GO  TO  290 

♦INITIAL  CALCULATIONS  FOR  CLUTTER  PSD  METHOD  2* 
UELR=CLGHT*TAUS/< 1000*2. ) 

BWF=1/(2.*RTS*1000/CLGHT> 

*F JND  GEOMETRY  AT  HORIZON* 

NAHOR=RTD*ASIN< RE/ <  RE+HS ) ) 

ECAH0R--90 » -NAHOR 

RSHOR=  ( RE+HS )  *SIN  ( DTR*ECAHOR ) 

KRMAX=RSHOR/DELR 

KDMAX=( ( < 2 . *VS*1000/LAMDA ) *COS ( ECAHOR*DTR ) ) )/PRFS  +  1 
1  =  1 

♦SEQUENCE  THRU  ALL  ISODOPS* 

DO  270  KD=1 , KDMAX , 1 

CVR= ( < ( KD-1 ) *PRFS+ ♦ 5*BWF )*LAMDA)/(2.*VS*1000) 

CVM=< < (KD-l )*PRFS)*LAMDA)/(2.*VS*1000) 

CVL=  < ( ( KD-1 ) *PRFS- . 5*BWF ) *LAMDA ) / <  2 . *05*1000 ) 


U  U  O  U  uu  u  u  u 


C  .  # FIN'D  MIN  RANGE:  OF  1 3 ODOR  CENTER* 

I D  E  P  A = R  T  D  *  A  C  0  S  (  C  V  M ) 

C  A  L.  L  T  G  E  C  S  R  (  H  S  ,  C  D  E  P  A ,  T  H  £  T  A  G  .EGA  N  G  .  R  S  ) 

:  •;  R  M I N  =  (  R  S  —  H  S )  /  D  E  i .  R  +  1 
R M I ND-HS + ( K R M I N - 1 ) * D E  L R 

C  *FIND  MIN  RANGE  OF  ISODOF  MIN* 

CDEPA-RTD*ACOS ( COL ) 

I F  ( CDEPA ♦ GT  *  90 . )  CDEPA=90 . 0 
CALL  TGECSR ( HS , CDEPA , THETAG f ECANG , RMINDN ) 

IF  (RMINDN  . LT .  HS)  RMINDN=HS 

SEQUENCE  THRU  ALL  I 30RANGES  FOR  PARTICULAR  ISODOF* 

DO  260  KR=KRM.I N ,  KRMAX  ,  1 
RS--HS+  (  KR-1  )  *DELR+RTS 

CNANG-RTD*ACOS ( ( ( HS+RE ) * ( HS+RE ) +RS*RS~RE*RE ) / <  2 . * ( HS  +  RE ) 
+  *RS ) ) 

ECANG=RTD*ACOS ( ( ( HS+RE ) * ( HS+RE ) +RE*RE-RS*RS ) / ( 2 . * ( HS+RE ) 
+  *RE )  ) 

THETAG=90. -CNANG-ECANG 
RG==ECANG*DTR*RE 
CDEPA=90 ♦ -CNANG 
CALL  ANGC( COM, CDEPA yCAANG) 

♦COMPUTE  ANTENNA  GAIN* 

CALL  ANTGN <  CANTS , CAANG , CDEPA , GMAXS , GMI NS , 

+  GCK1S,GCK2S,GCK3S,CNS,CMBS,BELS,ANGSLS,ANTAZS,DPANGS, 
+  GN1 f  GN2 1 GN3 » GN4 ) 

GN1=1 0 . ** ( ♦ 1*GN1 ) 

GN2=10.**( » 1*GN2) 

GN3= 1 0 . ** ( . 1*GN3> 

GN4"- 1 0  *  **  (  .  1  *GN4  ) 

♦FIND  GEOMETRY  OF  4  TRAPEZOID  CORNERS* 

♦FIND  GEOMETRY  OF  FURTHER  CORNERS* 

RSF-RS+ *  5*DELR 
IF  ( RSF , GT ♦ RSHOR )  RSF-RSHOR 
CALL  ANAD ( HS  t RSF , NAF ) 

DEPRF=90 ♦ -NAF 

CALL  TGECSR  (  HS  .  DEF'RF ,  ANG  ,  ECANGF  ,  RSF  ) 

CALL  ANGC  (  CVR  ,  DEF'RF  ,  AANGFR ) 

CALL  ANGC  (  COL  ,  DEP.RF ,  AANGFL  ) 

IF  ( AANGFR. GT. 90. )  AANGFR=90 , 0 
IF  ( AANGFR. LT.O. )  AANGFR-0 . 0 
IF  (AANGFL. GT. 90. )  AANGFL=90 . 0 
IF  (AANGFL. LT. 0.0)  AANGFL=0 . 0 

♦FIND  GEOMETRY  OF  NEAR  CORNERS* 

RSN=RS-,5*DELR 

IF(RSN.LT. RMINDN)  RSN=RMINDN 
CALL  ANAD  (HS,F'SN,  NAN) 

DEF'RN=90 .  -NAN 

CALL  TGECSR ( HS , DEPRN , ANG , ECANGN f  RSN ) 

CALL  ANGC (COR, DEPRN rAANGNR) 

CALL  ANGC (COL, DEPRN, AANGNL) 
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IF(AANGNR.GT .90.0)  AANGNR=90 . 0 

I F  <  A  A  N  G  N  R  ,  LT.0.0)  A  A  N  a  N  R =0,0 
I F  < AAN3NL . GT . 90 . 0 )  AANGNL=90 . 0 
I F  <.  AAN6NL  ,  LT  ,0.0)  AANGNL.-  0 . 0 

# COMPUTE  AREA  OF  CLUTTER  PATCH* 

BF= ( AANGFL-AANGFR > *DTR*RE*SIN ( DTR*ECANGF > 

BN= < AANGFL-AANGFR ) *BTR*RE*SI N  < DTRiECANGN ) 

ARACLT= (ECANGF-ECANGN ) *DTR*RE* , 5* ( BF+BN ) 

CALL  SIGMAO  ( THETAG ,  CTERS  ,  CSGF'DS  , SIGMA ) 
SIGMA=10.** (SIGMA*. 1  ) 

KCLTPU= ( KRADS*SIGMA*ARACLT*1 . E06 ) / ( ( RS*1000 ) **4 ) 
CLTF'SS  <  I )  =CLTPSS  (  I )  +KCLTPW*  (  GN 1  *GN1  +  GN2*GN2  + 
+  GN3*GN3  +  GN4*GN4) 

CONTINUE 
1  =  1  +  1 
CONTINUE 
RETURN 
END 
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*  »:  *  *  *■#:#***#***#####  #  #•  *  *  *  #  #*•#*  *  *  #  * 

*  SUBROUT  I  ME  SIGMA-0  * 

*  WRITTEN  BY ; JOE  DEVENUTC  * 

*  7  OCTOBER  1983  * 

********************************** 

THE  SUBROUTINE  SIGMA-0  CLCULATES  THE  5IGMA-0  VALUES 
GIVEN  THE  TERRAIN  TYPE, PDF  OF  SIGMA-O,  AND  GRAZING 
ANGLE  TO  CENTER  OF  CLUTTER  PATCH ♦ 

SIGMAO  CALLS  THE  FOLLOWING  IMSL  SUBROUTINES 
GGNML 
GGNL.G 

\ 

THE  FOLLOWING  VARIABLES  ARE  USED  IN  “SIGMAO* 

INPUT  VARIABLES 

CTERSS-INPUT  TO  CONTROL  TERRAIN  TYPE 
CSGPSS-INPUT  TO  CONTROL  SIGMA-0 ' S  PDF 
THETA-GRAZING  ANGLE  TO  CENTER  OF  CLUTTER  PATCH 
OUTPUT  VARIABLES 

SIGMAS-VALUE  OF  SIGMA-0  <DB> 

SUBROUTINE  VARIABLES 

SUBS  <  N ) -SUBROUTINE  VALUE  OF  SIGMA,  TERRAIN  DEP 

G ( N ) -ARRAY  CONTAINING  GROUND  SIGMA-0  VALUES 

S ( N >— ARRAY  CONTAINING  SEA  S1GMO  VALUES 

SAI (N) -ARRAY  CONTAINING  SNOW  AND  ICE:  SIGMA  VALUES 

ANGG ( N ) -ARRAY  OF  GRAZING  ANGLES 

NUMPT-NUMBER  OF  SIGMA-0  DATA  POINTS 

DSEED- INPUT  VARIABLE  TO  GAUSSIAN  DIVIATE  GENERATOR 

DSIG(N)  DELTA  TO  SIGMA  DUE  TO  PDF 

NR-NUMBER  OF  DIVIATES  TO  BE  GENERATED 

STD-STANDARD  DEVIATION  OF  SIGMA 

I-LOOP  VARIABLES' 

♦DECLARE  AND  INITIALIZE  VARIABLES* 

DIMENSION  G(50) ,S<50> , SAI (50) ,DSIG<5) , ANGG (50) , SUBS (50) 
REAL  THETA, SIGMAS, G, S, SAI, DSIG, STD, ANGG 
INTEGER  CTERSS , CSGPSS , NUMPT , NR , I 
DOUBLE  PRECISION  .DSEED 


♦SIGMAS  VALUES  GIVEN  THETA* 

DATA  G/-22, 5, -20. 5, -19, 5, -18. 0,-17. 5, -17. 5, -17.0, 
-16. 0,-15. 0,-14, 5, -12, 5, -10. 0,-8. 0,-6, 5, 

—4*5,— 2,5,— .5,1.5,2,5,31*0.0/ 

DATA  S/-45. 0,-41 .0,-38. 0,-35, 5, -33. 5, -32.0, 

-30 . 0 , -28 . 0 , -26 . 0 , -23 . 0 , -20 . 0 , 

— 18.0,— 15.5,  — 12.5,— 11 .0,— 3.0,  .5, 
6.5,13.0,31*0.0/ 

DATA  SAI/-31 .0,-29, 0,-28, 0,-26. 0,-25. 5, -23.0, 
-22.5,-20.0,-18.5,-17.5,-15.0,-14.5, 

-12. 5, -10. 5, -8. 5, -6. 5, -5. 0,-4. 0,-2. 5, 31*0,0/ 
DATA  ANGG/O. 0,5.0, 10.0, 15.0,20.0,25.0,30.0, 


35 •  0 • 40 ♦ 0 , 45 . 0 , 50 ♦ 0  *  55  *0-60. 0 . 65 , 0 * 
5  0  #  0  y  7  5  *  0  «*S0*0f85#  0  y  * 9  0  *0*3 1  %  0  *  0  / 

KjJJW  P"’ :  1  y 


♦  CONSTRAIN  ’'HET A  TO  LIMITS  OF  0 
IF  (THETA  . G  T .90)  THEN 
THE! A =90 . 0 
END  IF 

IF  (THETA  .  LT.  0)  THEN 
THET  A=0 . 0 
END  IF 


90  DEGRF 


♦DETERMINE  SIGMAS  FOR  TERRAIN  TYPE*' 
IF  ( CTERSS  . EQ .  1)  THEN 
DO  300  I =1, NUMPT, 1 
SUBS ( I >  =G  < I ) 

CONTINUE 

ELSE 

IF  (CTERSS  .  EP  ♦  2)  THEN 
DO  310  1=1. NUMPT , 1 
SUBS ( I ) =S ( I ) 

CONTINUE 

ELSE 

DO  320  I  =  1, NUMPT,  1 
SUBS ( I ) =3AI ( I > 

CONTINUE 
END  IF 
END  IF 


♦  LOCATE  SIGMAS  VALUE  WITHIN  ARRAY  OF  VALUES*: 
DO  330  1=2, NUMPT  , 1 

IF  (THETA  ►  LT . ANGG ( I ) )  GO  TO  340 
CONTINUE 

IF  ( ANGG  < I  —  1 >  .EQ.  ANGG(I))  THEN 
SIGMAS  =  SUBS ( I ) 

GO  TO  350 
END  IF  ' 


♦DO  2  POINT  LINEAR  INTERPOLATION  OF  SIGMAS* 
SIGMAS=SUBS ( 1-1 ) + ( SUBS ( I-l ) -SUBS ( I ) ) ♦ 

( THET A- ANGG ( 1-1 ) ) / ( ANGG ( I ) -ANGG ( I -1 > ) 

♦COMPUTE  SIGMAS- 0  DELTA  DUE  TO  PDF * 

IF  ( CSGPSS  . EQ , 1 )  THEN 
DSIG ( 1 ) =0 . 0 
GO  TO  360 
END  IF 
NR=1 

DSEED=123457 . ODO 
IF  (CSGPSS  . EQ . 2 )  THEN 

CALL  GGNML..  (  DSEED ,  NR  , DSIG) 

ELSE 

STD=1 

XM=SIGMAS 

•  CALL  GGNL.G  (DSEED,  NR,  XM.S,  DSIG) 
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*  subroutine:  toecsr  * 

*  WRITTEN  BY:  JOE  DEVENUTO  * 

*  7  OCTOBER  1983  * 

*  ♦♦♦♦♦*♦♦♦♦♦>?♦♦  ♦♦♦♦♦♦♦♦♦♦♦  ♦♦♦♦♦♦♦♦* 


THE  SUBROUTINE  TGECSR  COMPUTES  THE  GRAZING  ANGLE 
AND  THE  CORRESPONDING  EARTH  CENTRAL  ANGLE  AND  THE 
SLANT  RANGE  GIVEN  THE  ALTITUDE  (KM)  AND  THE  DEPRESSION 
ANGLE  * 

THE  FOLLOWING  VARIABLES  ARE  USED  IN  *  TGECSR* 

INPUT  VARIABLES 

HSS-ALTITUDE  OF  PLATFORM  (KM) 

DEPAG-DEPRESSION  ANGLE  (DEG) 

OUTPUT  VARIABLES 

ECANGS- EARTH  CENTRAL  ANGLE  (DEG) 

THETAS-GRAZING  ANGLE  (DEG) 

RSS-SLANT  RANGE  (KM) 

SUBROUTINE  VARIABLES 
NANG-NADIR  ANGLE 
ARG-ARGUMENT 

DECLARE  AND  INITIALIZE  VARIABLES 

REAL  HSS, DEPAG, ECANGS, THET AS, RSSt NANG 

COMMON  RE ,  G.RAV ,  CLGHT  ,  BK ,  DTR  ,  RTD ,  PI ,  KMTNM ,  NMTKM 

♦COMPUTE  THETAG* 

NANG-90. O-DEPAG 

ARG- ( SIN ( DTR*NANG ) * ( RE+HSS ) /RE ) 

IF  ( ABS ( ARG  > . GT . 1 ♦ 0 )  THEN 
ARG=  ARG- 1,0 
THET  AS-ASIN ( ARG ) ♦RTD 
GO  TO  396 
END  IF 

THETAS-  90.0  -  AS I N ( ARG ) ♦RTD 


♦  COMPUTE  EARTH  CENTRAL..  ANGLE* 
ECANGS-90 . O-NANG-THETAS 


♦COMPUTE  SLANT  RANGE* 

IF  ( NANG . EG . 0 .  )  THEN 
RSS-HSS 
GO  TO  399 
END  IF 

RSS=RE*SIN( ECANGS*DTR ) /SIN ( NANG*DTR ) 


RETURN 
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G  C  K 1 S  S  ,  G  C  K  2  S  3  ,  G  C  K  3  Jr  5  ,  C  N  3S.CMB  S  S  . 
A K A Z S  .  D  E P A S  ,  G N 1  S  .  G N2S  -  GN3S ,  G N - S ) 
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*  SUBROUTINE  ANTGN  * 

*  WRITTEN  BY:  C.R. RISER  * 

X  MODIFIED  BY?  JOE  DEVENUT0  * 

*  10  SEPTEMBER  1983  * 


X* XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

THE  SUBROUTINE  ANTGN  COMPUTES  THE  1-WAY  ANTENNA  GAIN 
TO  THE  CENTEF  OF  CLUTTER  PATCHES  IN  4  QUADRANTS  ON  THE 
EARTH  LOOKING  DOWNWARD  FROM  A  SBR  SYSTEM .  CLUTTER  PATCHES 
ARE  ASSUMED  TO  BE  SYMMETRICAL  WITH  THOSE  IN  THE  1ST  QUAD¬ 
RANT.  THE  ANTENNA  GAIN  IS  ASSUMED  TO  BE  CIRCULAR,  BEING 
FORMED  BY  A  FUNCTION  OF  REVOLUTION.  THE  ANTENNA  GAIN 
IS  COMPUTED  FROM  1  OF  4  DIFFERENT  FUNCTIONS  DEPENDING  ON 
THE  RADIAL  ANGLE  FROM  THE  MAIN  BEAM  BORESIGHT. 

THE  FOLLOWING  VARIABLES  ARE  USED  IN  'ANTGN1 
INPUT  VARIABLES 

CANTSS-INPUT  TO  CONTROL  ANTENNA  FUNCTION 
AZCTS -AZIMUTH  ANGLE  TO  CENTER  OF  CLUTTER  PATCH  (DEG) 
DEF'CTS -DEPRESS  I  ON  ANGLE  TO  CENTER  OF  CLUTTER  PATCH 
GMAXSS-MAX  1-WAY  GAIN  OF  ANTENNA  MAINBEAM  (DB) 
GMIN-SS-MIN  1-WAY  GAIN  OF  ANTENNA  SIDELOBES  (DB) 
GCK1SS-  ! 

GCK2SS -CONST ANT  TO  COMPUTE  1-WAY  GAIN  BETWEEN  MB  £  SL 
GCK3SS-ANGLE  SCALE  FACTOR  TO  COMPUTE  1-WAY  GAIN 
BETWEEN  MB  AND  SL 

BELSS-ANTENNA  ELEVATION  BEAMWIDTH  (DEG) 

CNSS-EXPONENT  TO  COMPUTE  ANTENNA  1-WAY  GAIN  BETWEEN 
ANGMB  AND  ANGSL 

C MBS S- ANTENNA  ANGLE  SCALE  FACTOR  FOR  COMPUTING 
MAINBEAM  GAIN 

ANGSLS-RADI AL  ANGLE  FROM  BORESIGHT  BEYOND  WHICH 
SIDELOBE  MAIN  IS  COMPUTED 

ANAZAS-AZIMUTH  ANGLE  OF  ANT,  MAINBEAM  BORESIGHT  (DEG) 
DEPAS -DEPRESSION  ANGLE  OF  MAINBEAM  BORESIGHT  (DEG) 
OUTPUT  VARIABLES 

GN1S-1 -WAY  GAIN  OF  ANTENNA  IN  1ST  QUADRANT 

GN2S-1-WAY  GAIN  OF  ANTENNA  IN  2ND  QUADRANT 

GN3S-1-WAY  GAIN  OF  ANTENNA  IN  3RD  QUADRANT 

GN4S-1-WAY  GAIN  OF  ANTENNA  IN  4TH  QUADRANT 

PROGRAM  VARIABLES 

ARG1 -ARGUMENT  IN  1ST  QUADRANT 
ARG2- ARGUMENT  IN  2ND  QUADRANT 
ARG3- ARGUMENT  IN  3RD  QUADRANT 
ARG4 -ARGUMENT  IN  4TH  QUADRANT 
ANG-ANGLE 

♦DECLARE  VARIABLES* 

REAL  AZCTS,  DEF'CTS,  GMAXSS , GMI NSS , GCK 1 SS , GCK2SS , GCK3SS , 
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B E L S S .  CNBS  «  CMBSS ,  AN03LS  ,  ANAZS .  D E P A 3 .  0 M 1  S  ,  G N 2 L  . 
G*»3S  « GN4S  j  ARC  1 1  ARG2  «  ARG3  »  ARG4  »  ANG 

1 1'  ■'  f  E  G  E  R  C  A  N  f  3  E> 

COMMON  RE »  GRAV  r  CLGHT » B K  ,  D T R  ,  R T D ,  P I , ! : M T N M .  N -IT K ■" 


* CHECK  GAIN  CALCULATION  METHOD* 
IF  ( CANTSS . EQ *  2  >  THEN 
GN1S«6MINSS 
GN2S-GHINSS 
GN3S=GMINSS 
GN4S=GMINSS 
GO  TO  495 
END  IF 


♦COMPUTE  ANTENNA  GAIN* 

♦COMPUTE  ANGLE  BETWEEN  MAINBEAM  BORESIGHT  AND  1ST  QUADRANT 
CLUTTER  PATCH  CENTER* 

ARG1 *SIN < DTR*DEPCTS ) *S IN ( DTR*DEPAS ) +COS ( DTR*DEPCTS ) * 

COS  (,  DTR*DEPAS  )  *COS  (  DTR*  ( AZCTS-ANAZS  )  ) 

IF  ( ABS  <  ARG 1 ) ♦ GT . 1 ♦ 0 )  THEN 
ANG=90.0 
GO  TO  405 
END  IF 

ANG=RTD*ACOS ( ARG1 ) 


♦FIND  ANGLE  SECTOR  FOR  FIRST  RADIAL  ANGLE  AND  COMPUTE 
ANTENNA  GAIN* 

IF  (ANG.LT.1.0E-5)  GO  TO  410 
IF  (ANG.LT.BELSS)  GO  TO  415 
IF  ( ANG . LT . ANGSLS )  THEN 
GO  TO  420 
ELSE 

GO  TO  425 
END  IF 
GN1S-GMAXSS 
GO  TO  430 

GN1S=(C0S(DTR*CMBSS*ANG> >**2,5 
GN1S-10. *AL0G10  <  GN1S  > 

GN1S=GN1S+GMAXSS 
GO  TO  430 

GN.1  S=GCK  1SS+GCK2SS/ (  ( SIN  (  DTR*GCK3SS*ANG  )  )  **CNSS  > 

GN1 S= 10 . *AL06 1 0 < GN 1 S ) 

GO  TO  430 
GN1S=GMINSS 


♦COMPUTE  ANGE  BETWEEN  MAINBEAM  BORESIGHT  AND 
CLUTTER  PATCH  CENTER  IN  2ND  QUADRANT* 

ARG2=SIN ( DTR*DEPCTS ) *SIN ( DTR*DEPAS > +COS ( DTR*DEPCTS ) * 
COS ( DTR*DEPAS ) *COS < DTR* (180. -AZCTS-ANAZS ) ) 

IF  (ABS  (ARG?.)  .GT,  1 .0)  THEN 
ANG=180 . 

GO  TO  435 
END  IF 

ANG=RTD*ACOS :ARG2) 


4  4  0 

445 


450 


C 

C 

C 


460 


470 

475 


*FIND  ANGLE  SECTOR  FOR  2ND  RADIAL  ANGLE  AND  COMPUTE 
A  N  T  E  N  N  A  G  A I N  * 

I F  (  A  N  G  ,  L  T  ,  1  .  E  -  5  )  G  0  T  C.:  4  4  0 
IF  ' ANG.LT.BELSS)  GO  TO  445 
I F  (  A N 3.LT *  A N G S L S )  T H E N 
GO  TO  450 
ELSE 

GO  TO  455 
END  IF 
GN2S=GMAXSS 
GO  TO  460 

GN2S= ( COS  < DTR*CMBSS*ANG ) ) **2 . 5 
GN2S= j 0 . *AL0G1 0 ( GN2S ) 

GN2S=GN2S+GMAXSS 
GO  TO  460 

GN2S=GCKlSS+GCK2SS/( ( S IN ( DTR*GCK3SS*ANG ) )**CNSS) 
GN2S=10.*AL0G10(GN2S> 

GO  TO  460 
GN2S-GMINSS 


♦COMPUTE  ANGLE  BETWEEN  MAINBEAM  BORESIGHT  AND 
CLUTTER  PATCH  CENTER  IN  3RD  QUADRANT* 

ARG'3~SIN ( DTR*DEPCTS >  *SIN ( DTR*DEPAS ) +COS ( DTR# DEPOTS ) * 
+  COS ( DTR*DEPAS ) *COS  <  DTR* (ISO. +AZCTS-ANAZS ) ) 

IF  ( ABS ( ARG3) .  GT ♦ 1 ♦ 0)  THEN 
ANG=90 ♦ 0 
GO  TO  465 
END  IF  - 

ANG=RTD*AC0S(ARG3>  .  1 

*FIND  ANGLE  SECTOR  FOR  3RD  RADIAL  ANTENNA  ANGLE  AND 
COMPUTE  GAIN* 

IF  (ANG.LT.l .E-5)  GO  TO  470 
IF  (ANG.LT.BELSS)  GO  TO  475 
IF  ( ANG . LT ♦ ANGSLS )  THEN 
GO  TO  480 
ELSE 

GO  TO  482 
END  IF 
GN3S=GMAXSS 
GO  TO  485 

GN3S=  <  COS ( DTR*CMBSS*ANG ) >**2.5 
GN3S~1 0 . *AL0G10 ( GN3S ) 

GN3S=GN3S+GMAXSS 
GO  TO  485 

GN3S=GCKlSS+GCK2SS/( <SIN(DTR*GCK3SS*ANG) )**CNSS) 
GN3S=10.*AL0G10(GN3S) 

GO  TO  485 
GN3S=GM JNSS 


♦COMPUTE  ANGLE  BETWEEN  MAINBEAM  BORESIGHT  AND  CLUTTER 
PATCH  CENTER  IN  4TH  QUADRANT* 

ARG4=S IN ( DTR*DEPCTS ) *SIN ( DTR*DEPAS ) +COS ( DTR*DEPCTS ) * 

+  COS ( DTR*DEPAS ) *COS ( DTR* ( 360 . -AZCTS-ANAZS ) ) 

IF  ( ABS ( ARG4 ) . GT . 1 . 0 )  THEN 
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ANG=90.0 
GO  TO  486 

TND  IF 

ANG=RTD*AC0S<ARG4> 

♦FIND  ANGLE  SECTOR  FOR  4TH  RADIAL  ANTENNA  AN! 
COMPUTE  GAIN* 

IF  ( ANG . LT ♦ 1 . E-5)  GO  TO  488 
IF  (ANG.LT.PELSS)  GO  TO  490 
IF  ( ANG . LT ♦ ANOSLS )  THEN 
GO  TO  492 
ELSE 

GO  TO  494 
END  IF 
GN4S=GMAXSS 
GO  TO  495 

GN4S= < COS ( DTR*CMBSS*ANG ) )**2.5 
GN4S=10 ♦ *ALQG1 0 ( GN4S ) 

GN4S=GMAXSS+GN4S 
GO  TO  495 

GN4S=GCK1SS+GCK2SS/ ( <  SIN ( DTR*GCK3SS*ANG ) ) **CNSS ) 
GN4S= 10 . *AL0G10  <  GN4S ) 

GO  TO  495 
GN4S=GMINSS 


RETURN 


46 


MSS  T SS  , 


REAL  TSS 

INTEGER  NSS  ,  I » J 

C  r  ■  ~i  P  L  E  X  3  5  c  0  E  S  <  2  0  0 ) 


DO  460  1-1,30,2 

SGCOES  < I  }-•;  1 .0,0.0) 
SGCOES ( I + 1 > «CMPLX ( 0 . , 0 . ) 
460  CONTINUE 
J=NSS-30 
DO  465  I«31,J,1 

SGCOES  ( I )  =CMF’LX  ( 0 .  ,  0  .  ) 
CONTINUE 
RETURN 
END 
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THE  SUBROUTINE  CLTARA  CALCULATES  THE  AREA  3F 
CLUTTER  SEEN  BE  A  SPACE-BASED  RADAR  SYSTEM  GIVEN  “HE 
PLATFORM  ALTITUDE, RADAR  AZIMUTH  AND  ELEVATION  BEAM- 
WIDTHS,  THE  GRAZING  ANGLE  AND  THE  SLANT  RANGE. 

THE  FOLLOWING  VARIABLES  ARE  USED  IN  *  CLTARA  * 

INPUT  VARIABLES 

HSS-SATELLITE  ALTI TUBE ( KM ) 

BAZSS-ANTENNA  AZIMUTH  BEAMWIDTH  (DEG) 
BELSS-ANTENNA  ELEVATION  BEAMWIDTH  (DEG) 
THETSS-GRAZING  ANGLE  (DEG) 

RSS-SLANT  RANGE  (KM) 

OUTPUT  VARIABLES 

ARACLS-CL UTTER  AREA  (KMHcKM) 

PROGRAM  VARIABLES 
ANG-ANGLE  (DEG) 

♦DECLARE  VARIABLES* 

REAL  HSS , BAZSS , BELSS , THETSS , RSS , ARACLS  ,  ANG 
COMMON  RE , GRAV , CLGHT , BK , DTR , RTD , PI , KMTNM , NMTKM 


♦COMPUTE  CLUTTER  AREA* 

ANG=ASIN( ( RE/ ( RE+HSS ) ) *COS ( DTR*THETSS ) ) 
ARG-  (  ( HSS+RE  )  /RE  )  *SIN  (  ANG  ) 

IF  (ABS(ARG) .GT. 1 .0)  THEN 
ARG=2 . O-ARG 
END  IF 

ARACLS= ( ACGS ( ARG ) - ( THETSS+BELSS ) *DTR ) * 
BAZSS*DTR*RE*RSS 


( BELSS*DTR ) 


RETURN 
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THE  SUBROUTINE  A NAD  COMPUTES  THE  NADAR  ANGLE  USING 
THE  LAW  OF  COSINES  GIVEN  THE  SLANT  RANGE  AND  SATELLITE 
ALTITUDE. 

THE  FOLLOWING  VARIABLES  ARE  USED  IN  'ANAD* 

INPUT  VARIABLES 

HSS-SATEL.LITE  ALTITUDE 
RSS-SLANT  RANGE 
OUTPUT  VARIABLES 
NAS -NADIR  ANGLE 
ARG-ARGUMENT 

♦DECLARE  VARIABLES* 

REAL  HSS . RSS  r  NAS  r ARG 

COMMON  RE , GRAV , CLGHT , BK , DTR » RTD , PI , KMTNM , NMTKM 
♦COMPUTE  NADIR  ANGLE* 

ARG=  (  (  (  ( HSS+RE  >  **2 )  +  ( RSSYY2 )  - <  REYY2  >  )  / 

( 2* ( HSS+RE ) YRSS ) ) 

IF  (ABS(ARG) .GT.l.O)  THEN 
ARG=2 . O-ARG 
END  IF 

NAS=RTD*ACOS< ARG) 

RETURN 

END 
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THE  SUBBROUTINE  ANGC  COMPUTES  THE  AZIMUTH  ANGLE 
BETWEEN  THE  DEPRESSION  ANGLE  VECTOR  AND  THE  CVX  VECTOF 


THE  POLL WO I NG  VARIABLES  ARE  USED  IN  *  ANGC ’ 
INPUT  VARIABLES 

CVX= ( FT  + <  KD-1 ) APRF+NABWF ) ALAMDA ) / <  2 . * VS ) 
DEPAS-DEPRESSION  ANGLE 
OUTPUT  VARIABLE 

AZANG-AZIMUTI-I  ANGLE 
PROGRAM  VARIABLES 
DENOM - DENOMINATOR 
ARG-AEGUMENT 

♦DECLARE  VARIABLES 
REAL  CVX , DEPAS  r AZANG 

COMMON  RE  ,  GRAV ,  CLGHT  ,  BK ,  DTR ,  RTD ,  PI ,  K'MTNM  ,  NMTK'M 

PENOM=COS  < DTR*DEPAS ) 

IF ( ABS ( DENOM ) .LT. IE-148)  THEN 
DEN0M=SIGN(1E-1 48, DENOM) 

END  IF 

ARG=CVX/DEN3M 
IF  <  ABS ( ARG ) . GT . 1 . 0 )  THEN 
ARG=2 . -ARG 
END  IF 

20  AZANG=RTD* ACOS  (  ARG  ' 


RETURN 


r  IF  \ 
CCNSHS 
«  7  > 


* COMPUTE  UNIFORM  WEIGHTS* 
K  *  ( (-1 )**I) 


c 

c 

c 

c: 
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s  u  b  r  o  u  t  :::  n  e  c  a  n  c  e  l  :  c  c  n  s  h  s  ,  n  s  ,  t  s ,  b  u  r  s  l 

t  FLMAGS, TAPS) 
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THE  SUBPROGRAM  CANCEL  SIMULATES  A  GIVEN 
CLUTTER  FIlTER  AND  THEN  PERFORMS  THE  FILTERING 

ON  THE  CLUTTER  PSD  CALCULATED  IN  SUBPROGRAM 
GLUT. 


CANCEL  CALLS  THE  FOLLOWING  SUBROUTINES 
FFT2C 
SIGNAL 


THE  FOLLOWING  VARIABLES  ARE  USED  IN  •CANCEL" 

INPUT  VARIABLES 

CCNSHS- INPUT  TO  CONTROL  CANCELLATION  SCHEME 
NS-MUMBER  OF  POINTS 
TS-SAMPLING  INTERVAL 

CLTPSS <N) -ARRAY  OF  CLUTTER  PSD  VALUES 
BURSS-BURST  TIME  (SEC) 

TAPS-NUMBER  OF  FILTER  TAPS 
OUTPUT  VARIABLES 

FLCPSS(N) -ARRAY  OF  FILTERED  CLUTTER  PSD  VALUES 
FLMAGS ( N  ) -ARRAY  OF  FILTER' MAGNITUDE  COEF 
SUBPROGRAM  VARIABLES 

FLCOEF ( N ) -ARRAY  OF  COMPLEX  FILTER  COEFFICIENTS 
REVAL-REAL  PART  OF  COMPLEX  VARIABLE 
IMVAL- IMA JINARY  PART  OF  COMPLEX  VARIABLE 
I f  J  >  K  t L-LOOP  VARIABLES 
I WK ( N ) -WORK  ARRAY  FOR  FFT2C  SUBROUTINE 
FLTAMG ( N ) -ARRAY  OF  H(W>  FILTER  COEF . 

SGCOEF(N) -ARRAY  OF  COMPLEX  SIGNAL  COEF » 

DELT-DELTA  TIME 

DENQM-DENQMINATOR 


♦DECLARE  AND  INITIALIZE  VALUES* 

DIMENSION  CLTPSS ( 200 ) , FLCPSS ( 200 ) , IWK(200) , 
FLMAGS <  200 ) , FLTAMG  <  200 ) 

REAL  TS , CLTPSS , FLCPSS , REVAL , IMVAL , FLMAGS , 
FLTAMG  , K r  BURSS  , BELT , NUM , DENOM 
INTEGER  CCNSHS, NS, I,L, J,IWK,TAPS 
COMPLEX  FLCOEF  <200),SGC0EF(200) 

COMMON  RE , GRAV , CLGHT , BK , DTR , RTD , PI , KMTNM , NMTKM 
DO  605  1=1,150,1 
FLTAMG ( I ) =0 . 0 
CONTINUE 


*CHOOSE  WEIGHTING  SCHEME  AND  DETERMINE  TAP  WEIGHTS* 
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DELT=TS /(TAPS+l) 

P R 3: ; n*,  •'  D E LT -  •'  ,  DEL  T  ,  ' C Z  ■  3 H S=  '  ,  Z C N S H S 
♦BARTLETT  WEIGHTS* 

IF  i CCNSHS. EQ.l ) THEN 
DO  600  1=0 , TAPS  , 1 

K= ( 1 - I *DELT/BURSS ) * ( ( - i ) ** I ) 

FLCOEF (  1  +  1 >=CMPLX(N,0.0) 

600  CONTINUE 
GO  TO  650 
END  IF 

♦HANNING  WEIGHTS* 

IF  ( CCNSHS ♦ EG ♦ 2 )  THEN 
DO  610  1=0, TAPS, 1 

K=(  » 5+ » 5* ( COS ( I *DELT*F‘I/BURSS )  )  )*<  (-!)**!  ) 
FLCOEF (  I  + 1 ) =CNPLX ( K , 0 ♦ 0 ) 

610  CONTINUE 
GO  TO  650 

END  IF 

♦HAMMING  WEIGHTS* 

IF  ( CCNSHS. EQ. 3)  THEN 
DO  620  1=0, TAPS, 1 

K=  ( ♦ 54+ . 46* ( COS  < I *DELT*PI/BURSS ) ) )*( (-1)**I) 
FLCOEF < I + 1 ) =CMPLX ( K  ,  0 . 0  ) 

620  CONTINUE 
GO  TO  650 

END  IF 


♦BLACKMAN  WEIGHTS* 

IF  ( CCNSHS, EQ. 4)  THEN 

DO  630  1=0 , TAPS , 1  : 

K= ( . 42+ . 5* ( COS ( I*DELT*PI/BURSS ) )  + 

+  ♦  08*  ( COS ( 2*I*DELT*PI/BURSS ) ) >*( <-l)**I) 

FLCOEF ( 1  + 1 ) =CMPLX <  K ,0 ♦ 0 ) 

630  CONTINUE 
GO  TO  650 

END  IF 

♦BINOMIAL  WEIGHTS* 

IF  (CCNSHS. EQ. 5)  THEN 
DO  635  1=1, TAPS+l, 1 
NUM=T  APS 
DEN0M=1 
J=T  APS-1 

636  IF  (J.GT.O)  THEN 

NUM=NUM* J 
J=J-1 
GO  TO  636 

END  IF 

J=TAPS-I+1 

637  IF  (J.GT.O)  THEN 

DENOM=DENOM*J 

J=J-1 

GO  TO  637 

END  IF 

J=I-1 

638  IF  (J.GT.O)  THEN 

DENOM=DENOM# J 


>vv 


J=J-1 
GO  TO  638 

ENT  IF 

K= ( NUH/DENOM >*<(-l> *# (1-3 ) > 

F L C 0 E F ( I ) -  C MR LX ( K , 0 .  > 

635  CONTINUE 
GO  TO  650 

END  IF 

C  * IMPULSE  RESPONSE* 

IF  C  CCNSMS . EG . 6 )  THEN 
CALL  S IGNAL  ( NS  ,  TS ,  SGCOEF  ) 

DO  640  1  =  1 , NS  , 1 

FLCOEF ( I )  =  ( SGCOEF  ( I ) ) 

640  CONTINUE 
GO  TO  660 

END  IF 

C  *UNIFGRM  WEIGHTS* 

IF  ( CCNSHS . EQ ♦ 7 )  THEN 
DO  645  1=0, TAPS, 1 
K  =  (  - 1  )  *  *  I 

FLCOEF ( 1  +  1 )  =  CMPLX ( K , 0  » 0 ) 

645  CONTINUE 
END  IF 

C  *BUTTER WORTH  WEIGHTS* 

IF  ( CCNSHS, EQ. 8)  THEN 
FLCOEF (1>  =  < 5. 1258,0.) 

FLCOEF  <2)  =  (-13. 1371,0. ) 

,  FLCOEF (3)=(21 .8462,0*) 

FLC0EF(4)=(-25, 6884,0. ) 
FLCOEF  ( 5  >  =  <  21  ♦  8462*,  0  .  ) 
FLCOEF (6)=(~13.1371,0. ) 
FLCOEF <7 >=(5.1258,0. > 

FLCOEF ( 8  >  =  ( -1 , ,0.  ) 

END  IF 

C  *ZERO  FILL  FILTER* 

650  L=NS- <  T  APS+1 > 

DO  655  I = 1 , L , 1 

FLCOEF ( I+(TAPS+1 > )=CMPLX ( 0 . 0 , 0 . 0 > 

655  CONTINUE 

J= ( AL0G1 0 ( REAL ( NS  > > ) / ( AL0G10 ( 2  .  >  > 

660  CALL  FFT2C ( FLCOEF , J, I WK > 

DO  670  1=1 , NS , 1 

REVAL=REAL ( FLCOEF ( I ) > 

IMVAl..=AIMAG  (  FLCOEF  ( I  )  > 

FLMAGS( I >=( (REVAL**2>  +  ( IMVAL**2 > > ** . 5 
670  CONTINUE 
C 

C  *COMPUTE  FILTER  CLUTTER  PSD  VALUES* 

DO  690  1=1, NS, 1 

FLCF’SS  ( I  )  =  FLMAGS  ( I )  *FLMAGS  <  I  )  *CLTF'SS  ( I ) 
690  CONTINUE 
RETURN 
END 
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SUE  :'PG  CRAM  CLTATT  * 
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THE  SUBPROGRAM  CLTATT  COMPUTES  THE  CLUTTER  ATTENUATION 
FACTOR  GIVEN  THE  CLUTTER  PSD  AS  SEEN  BY  A  SPACE-BASED  FiAEi 
AR  AND  THE  CLUTTER  CANCELLATION  SCHEME. 

THE  FOLLOWING  VARIABLES  ARE  USED  IN  ‘CLTATT" 

INPUT  VARIABLES 

CLTPSS(N) -ARRAY  OF  CLUTTER  PSD  VALUES 
FLCPSS CN> -ARRAY  OF  FILTERED  CLUTTER  PSD 
FCOFMG(N) -ARRAY  OF  FILTER  COEFFICIENTS 
NS-NUMBER  OF  DATA  POINTS 
OUTPUT  VARIABLES 

CTATS-CLUTTEF:  ATTENUATION  FACTOR 
S U B P R 0 G R A H  VARIABLES 

F CL F’W -FILTERED  CLUTTER  POWER 
CL TP W- RECEIVED  CLUTTER  POWER 
I,J,L-LOOP  CONTROL 

DIMENSION  CLTF'SS ( 200 ) , FLCPSS < 200 )  , FSGPSD < 200 > 

DIMENSION  IWK ( 200) , FFTSIG < 200 ) ,FC0FMG(200> 

REAL  VTS ,  RTS  ,  RCSS  ,  KRADS  ,  CLTF'SS ,  FLCPSS  ,  FCLTF'W  ,  CL..TPW , 

+  FSIGF'W  ,  SIGF'W  ,  GAIN ,  GMAXS  , CTATS »  FSGPSD , 

+  PCSNDB ,  REVAL. ,  IMVAL ,  FFTS IG ,  FCOFMG 
COMPLEX  SGCOES ( 200 ) 

INTEGER  NS,I ,K, J,L,IWK 

*ZERO  POWER  VARIABLES* 

FCLTF'W-0 . 0 
CLTF’W-0 . 0 

* COMPUTE  CLUTTER  POWER  IN  %  OUT  OF  FILTER* 

DO  800  1=1 , NS , 1 

FCLTPW=FCL.TPW  +  FLCF'SSd) 

CLTPW=CLTPW  +  CLTF'SS  ( I ) 

CONTINUE 


♦COMPUTE  ATTENUATION  FACTOR* 
CTATS=1 0*AL0G1 0 ( CLTPW/FCLTPW ) 


RETURN 
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